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Introduction 


The  First  Conference  on  Defense  Meteorological  Satellite  Program  (DMSP) 
Retrieval  Products  was  planned  by  scientists  at  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)'s  National  Meteorological  Center  (NMC),  NOAA’s  National 
Environmental  Satellite,  Data,  and  Information  Service  (NESDIS),  and  the  Department  of 
Defense  (DoD).  The  conference  was  held  14-15  April  1992  at  the  NOAA  Science  Center  in 
Washington,  DC.  The  objective  of  the  conference  was  to  bring  together  individuals  from 
both  the  defense  and  civilian  communities  with  a  common  goal  of  understanding  the  use  of 
DMSP  retrieval  products.  The  format  consisted  of  four  sessions  of  contributed  papers  and 
a  discussion  session  resulting  in  a  series  of  recommendations.  The  conference  Steering 
Committee  consisted  of  E.  Kalnay  (NOAA/NMC),  P.  K.  Rao  (NOAA/NESDIS),  and  R. 
McClatehey  of  the  Air  Force  Phillips  Laboratory  Geophysics  Directorate  (PL/GP).  The 
meeting  was  organized  by  Atmospheric  and  Environmental  Research,  Inc.  (AER). 

The  conference  was  attended  by  over  90  scientists  interested  in  the  generation  and 
application  of  DMSP  retrieval  products.  Among  those  represented  were  the  government 
and  academic  research  community,  the  major  operational  numerical  weather  prediction 
centers,  and  industry.  The  research  community  was  represented  by  the  Air  Force  Phillips 
Laboratory  (PL),  the  Naval  Research  Laboratory  (NRL),  NOAA/NESDIS  and  the  National 
Aeronautics  and  Space  Administration  (NASA).  Operational  centers  included  the  NMC, 
the  Navy's  Fleet  Numerical  Oceanography  Center  (FNOC),  Air  Force  Global  Weather 
Central  (AFGWC),  the  European  Centre  for  Medium  Range  Weather  Forecasts  (ECMWF), 
and  Canada's  Atmospheric  Environment  Service  (AES).  DoD  was  also  represented  by 
Headquarters,  Air  Force  Space  Systems  Division  and  the  Air  Weather  Service.  Industry 
representatives  included  those  from  GE  Astro-Space,  Lockheed,  Hughes,  and  ITT. 

This  proceedings  volume  summarizes  the  technical  presentations  at  the  conference. 
The  main  part  of  the  report  consists  of  the  abstracts  and  copies  of  the  viewgraphs  or  slides 
and  other  material  as  provided  by  the  authors  for  the  presentations.  Preceding  the 
viewgraphs  for  each  session  is  a  commentary  prepared  by  the  co-chairpersons  on  each 
presentation  within  that  session.  The  Appendix  includes  the  original  call  for  papers  and 
invitation  to  the  meeting,  a  copy  of  the  Agenda  for  the  meeting,  a  list  of  the  attendees,  and 
an  author  index. 

We  acknowledge  the  important  contributions  of  the  session  chairmen  E.  Kalnay 
(NOAA/NMC),  G.  Ohring  (NOAA/NESDIS),  and  R.  McClatehey  (PL/GP).  Special 
thanks  are  also  extended  to  E.  Stenhouse  of  AER  and  Pam  Taylor  of  NESDIS. 

Received  for  publication  21  July  1992 
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Session  1:  Overview  of  Current  and  Future  DMSP  Sensors  and  Products 

(Isaacs,  chair) 

The  first  session  provided  background  information  on  current  operational  DMSP 
sensors  and  near  term  and  far  term  plans.  James  Hollinger  of  the  Naval  Research 
Laboratory  provided  a  status  report  on  the  performance,  calibration,  geolocation,  and 
algorithm  validation  of  the  on  orbit  Special  Sensor  Microwave  Imagers  (SSM/I).  Of 
particular  note  was  data  on  the  intercalibration  of  the  three  SSM/I's  aboard  F-08,  F-10,  and 
F- 1 1 .  F-08  and  F-10  compare  well,  while  this  is  not  so  true  for  F-10  and  F- 1 1 . 
Implications  were  described  for  the  application  of  the  calibration/validation  retrieval 
algorithms. 

Mike  Griffin  of  the  Phillips  Laboratory  presented  an  overview  of  the  calibration/ 
validation  program  for  the  recently  launched  DMSP  microwave  humidity  profiler 
(SSM/T-2).  The  SSM/T-2  is  a  cross-track  scanning  microwave  radiometer  operating  in  the 
frequency  range  from  90  to  183  GHz.  The  mission  of  the  SSM/T-2  is  to  provide  high 
quality  vertical  moisture  profiles  that  can  be  made  available  to  a  variety  of  models  running 
operationally  at  the  Air  Force  Global  Weather  Central  (AFGWC).  The  SSM/T-2  sounder 
was  launched  late  in  1991.  Calibrations  efforts  will  be  performed  by  the  Phillips 
Laboratory  and  validation  efforts  by  the  Aerospace  Corporation.  Calibration  efforts  for  the 
SSM/T-2  microwave  moisture  sensor  will  be  based  on  aircraft  underflights  and  the  use  of 
radiative  transfer  models.  Imagery  from  the  DMSP  Operational  Linescan  System  (OLS) 
will  be  used  to  help  identify  heavy  cloud  areas  which  may  impact  on  water  vapor  sounding 
accuracies. 

Bruce  Thomas  of  the  Aerospace  Corporation  presented  a  companion  paper  on  the 
SSM/T-2  moisture  sounding  algorithm  validation  effort.  The  algorithm  currently  used  at 
AFGWC  to  extract  the  moisture  profiles  from  the  radiances  is  multiple  linear  regression. 
Historically,  regression  coefficients  have  been  organized  using  simple  static  discriminates 
(i.e.  latitude  and  longitude  as  for  the  SSM/T-1).  The  SSM/T-2  breaks  with  this  transition 
and  attempts  to  build  a  "dynamically  based"  set  of  piece  wise  linear  recession  coefficients 
stratified  by  air  mass.  A  brief  discussion  on  the  choice  of  discriminants  for  air  mass  typing 
was  presented.  Early  analysis  of  SSM/T-2  data  shows  the  instrument  to  be  quite  stable 
with  sensor  noise  below  that  of  the  original  specifications.  Aerospace  is  in  the  process  of 
validating  the  SSM/T-2  derived  moisture  profiles  (against  synoptic  time  and  special 
radiosondes). 

Capt.  Mark  Borden  of  Headquarters,  Air  Force  Space  Systems  Division  described 
plans  for  the  next  generation  of  DMSP  weather  satellite:  the  DMSP  Block  6.  The  Block  6 
program  features  an  aggressive  and  innovative  data  driven  strategy  for  the  acquisition  of  the 
next  generation  DoD  weather  satellite.  As  a  baseline  it  maintains  the  operational  utility  of 
the  DoD  weather  satellite  system  while  seeking  to  optimize  operational  efficiency  within  a 
specified  life  cycle  cost  cap.  It  features  the  development  of  a  complete  meteorological 
satellite  system,  embracing  both  space  based  and  ground  based  assets.  Particular  attention 
is  focused  on  the  insertion  of  maturing  technologies  into  the  developing  Block  6  effort. 

The  system  design  concept  responds  to  the  data  needs  of  the  operational  DoD  user.  System 
efficiencies  are  engendered  through  the  use  of  data  fusion,  improved  mass  memory, 
sensors,  power  concepts  and  the  application  of  Knowledge  Based  Systems  (KBS)  to 
optimize  user  efficiency. 
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SESSION  1:  OVERVIEW  OF  CURRENT  AND  FUTURE  DMSP  SENSORS 
AND  PRODUCTS 


SSM/I:  STATUS  REPORT 


James  P.  Hollinger 


Radio,  Infrared  and  Optical  Sensing  Branch 
Center  for  Advanced  Space  Systems 
Naval  Research  Laboratory 
Washington,  DC  20375 


There  are  currently  two  SSM/I's  in  orbit,  the  Defense  Meteorological  Satellite 
Project  Block  5D2  F-08  and  F-10  satellites.  A  third  is  scheduled  for  launch  in  November 
1991.  An  additional  four  SSM/I’s  along  with  five  SSM/IS's,  a  follow  on  version  of  the 
SSM/I  which  incorporates  atmospheric  temperature  and  humidity  sounding  channels,  will 
provide  data  into  the  twenty  first  century. 

A  brief  description  of  the  SSM/1  will  be  presented  followed  by  a  status  report  on  the 
performance,  calibration,  geolocation  and  algorithm  validation  of  the  SSM/I's  in  orbit. 
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SSM/I  CALIBRATION  AND  VALIDATION 
OF  SN004  ON  DMSP  F-ll 


•  SENSOR  HEALTH 

-  SCAN  PERIOD  STABILITY 

-  TEMPERATURE  STABILITY 

-  SYSTEM  NOISE 

-  GAIN  STABILITY 

-  REFERENCE  LOAD  STABILITY 

•  REGION  COMPARISONS 

-  AMAZON  JUNGLE 

-  ARABIAN  DESERT 

-  GREENLAND 

-  CLEAR  CALM  OCEAN 

-  SIMULTANEOUS  DATA 

•  GEOLOCATION 


F-ll 


•  NEAR  CIRCULAR  ORBIT 

•  STABLE  -  NO  GAIN  STEPPING 

•  PERIOD 

-  SHORT  TERM  1.899  +/-  0.002  sec 

-  LONG  TERM  1.89900  sec 

•  IIOT/COLD  LOAD  SAMPLES  UNIFORM 

•  HOT  LOAD  TEMPERATURE  UNIFORM  +/-  0.1  K 

•  SYSTEM  NOISE  EXCELLENT 

•  GEOLOCATION 

-  WITHIN  SPEC  (±7  Km) 

-  BEAMS  COLLOCATED 

•  CALIBRATION 

-  22V  and  85V  &  H  APPEAR  -  1  TO  2  DEG  HIGH 

-  OTHER  CHANNELS  <  0.5  DEG 
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F-8  F-10  F-ll 


REV  18532 

REV  719 

REV  1276 

ASC  EQUATOR  X-ING  TIME 

06:35 

19:42 

17:04 

PERIOD 

101.8  MIN 

100.7  MIN 

101.9  MIN 

ALTITUDE  MAX 

882  Km 

853  Km 

878  Km 

ALTITUDE  MIN 

838  Km 

740  Km 

841  Km 

ECCENTRICITY 

0.00145 

0.00814 

0.00129 

SEMI  MAJOR  -  SEMI  MINOR 

8  M 

238  M 

6  M 

FOCAL  DISTANCE 

21  Km 

117  Km 

19  Km 

SWATH  MAX 

1480  Km 

1427  Km 

1483  Km 

SWATH  MIN 

1400  Km 

1226  Km 

1414  Km 

INCIDENCE  ANGLE  MAX 

53.60  DEG 

53.29  DKG 

53.56  DEG 

INCIDENCE  ANGLE  MIN 

53.13  DEG 

52.10  DEG 

53.16  DEG 

MAX  INCIDENCE  ANGLE 
CHANGE  FOR  ALL  ORBITS 

0.9  DEG 

1.4  DEG 

0.5  DEG 

RADIOMETER  SENSITIVITIES 
HOT  LOAD  TARGET 

19V 

1911 

22V 

37V 

37H 

85V 

85H 

SPEC 

0.8 

0.8 

0.8 

0.6 

0.6 

1.1 

1.1 

F-8 

(Aug  1987) 

0.37 

0.37 

0.58 

0.30 

0.33 

0.69 

0.59 

F-10 

(REV  6632) 

0.50 

0.48 

0.54 

0.37 

0.37 

0.53 

0.57 

F-ll 

(REV  1764) 

0.46 

0.39 

0.55 

0.34 

0.35 

0.58 

0.44 
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AMAZON 


F-ll  AND  F-1C  NKAK  COINCIDENT 
DATA  FOR  18  PASSES 

AVERAGE  TIME  DIFFERENCE  OF  187  MINUTES 


CHANNEL 

F-ll 

MEAN 

F-10 

MEAN 

MEAN 

DELTA 

(F11-F10) 

RMS 

19V 

285.65 

286.06 

-0.65 

1.58 

1911 

284.10 

284.26 

-0.43 

1.69 

22V 

286.03 

284.36 

1.48 

1.29 

37V 

281,17 

282.29 

-1.31 

1.60 

37H 

280.14 

280.94 

-1.03 

1.69 

85V 

283.04 

283.56 

-0.58 

2.16 

8511 

282.96 

282.99 

-0.09 

2.40 

ARABIAN  DESERT 

F-ll  AND  F-10  NEAR  COINCIDENT 
DATA  FOR  18  PASSES 

AVERAGE  TIME  DIFFERENCE  OF  183  MINUTES 


CHANNEL 

F-ll 

MEAN 

F-10 

MEAN 

MEAN 

DELTA 

(F11-F10) 

RMS 

19V 

287.53 

286.24 

1.41 

1.76 

19H 

246.79 

245.52 

1.26 

1.95 

22V 

286.61 

283.44 

3.30 

1.76 

37V 

282.47 

281.95 

0.71 

2.40 

37H 

249.60 

248.60 

1.14 

2.32 

85V 

281.70 

279.80 

2.05 

3.26 

8511 

266.63 

264.18 

2.61 

3.09 

5 


GREENLAND 


F-li  AND  F-iO  NEAR  COINCIDENT 
DATA  FOR  21  PASSES 

AVERAGE  TIME  DIFFERENCE  OF  132  MINUTES 


CHANNEL 

F-ll 

MEAN 

F-10 

MEAN 

MEAN 

DELTA 

(F11-F10) 

RMS 

19V 

192.28 

192.71 

-0.57 

0.99 

1911 

141.54 

141.10 

-0_}8 

0.94 

22V 

186.94 

185.95 

0.87 

0.85 

37V 

175.76 

176.14 

-0.48 

0.53 

37H 

145.00 

145.08 

-0.09 

0.63 

85V 

191.03 

190.10 

0.82 

0.78 

8511 

169.14 

167.67 

1.44 

1.26 

CLEAR  CALM  OCEAN 

F-ll 

MODEL 

DELTA  (SSM/I 

-  MODEL) 

F-ll 

F-10 

F-08 

19V 

179.2  ±  0.7 

177.7 

1.5 

1.0 

1.1 

1911 

100.6  ±  1.0 

99.7 

0.9 

-0.3 

0.9 

22V 

189.3  ±  1.4 

187.1 

2.2 

1.6 

0.5 

37V 

200.9  ±  0.5 

205.2 

-4.3 

-2.3 

-2.8 

3711 

129.2  ±  0.9 

129.3 

-0.1 

-1.4 

0.3 

85V 

236.4  ±  1.1 

239.0 

-2.6 

-1.9 

-4.3 

8511 

172.5  ±  2.0 

173.7 

-1.2 

0.7 

-1.1 

0 


SIMULTANEOUS  DATA  MEAN  DIFFERENCES 


CHANNEL  F10  -  F08  Fit  -  HO 


MEAN 

RMS 

MEAN 

RMS 

I9V 

0.5 

1.2 

-0.2 

0.1 

1911 

0.3 

0.6 

0.3 

0.5 

22V 

-0.4 

1.1 

1.5 

0.3 

37V 

0.5 

1.3 

-0.5 

0.3 

37H 

0.8 

0.7 

0.0 

0.4 

85V 

_ _ 

0.8 

0.2 

85H 

1.3 

0.2 

SUMMARY  OF  MEAN  DIFFERENCES  FOR 
ARABIAN  DESERT,  AMAZON  JUNGLE,  GREENLAND. 
CLEAR  CALM  OCEAN,  AND  SIMULTANEOUS  DATA 

CHANNEL  F10-F08  Fll  -  F10 


MEAN 

RMS 

MEAN 

RMS 

19V 

0.0 

0.5 

0.0 

1.6 

1911 

-0.5 

0.7 

0.4 

1.6 

22V 

-0.3 

0.4 

1.6 

1.8 

37V 

0.4 

0.6 

-0.5 

1.7 

37H 

0.1 

0.9 

0.1 

1.7 

85V 

0.3 

1.5 

0.6 

2.3 

85H 

0.4 

1.2 

1.0 

2.6 
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ALGORITHM  RETREIVAL  ERRORS  FOR 
EQUAL  CALIBRATION  ERROR  FOR  EACH  CHANNEL 


CALIBRATION 

RETREIVAL 

ALGORITHM 

ERROR 

ERROR 

WIND  SPEED 

0.5 

1.13  m/s 

WATER  VAPOR 

0.5 

0.38  kg/m1 

CLOUD  WATER 

0.5 

0.01  kg/m1 

RAIN  (OCEAN) 

0.5 

0.07  mm/hr 

ICE  (WINTER) 

0.5 

0.99  % 

ICE  (SUMMER) 

0.5 

1.50  % 

ICE  (37  GHz) 

0.5 

0.56  % 

RAIN  (LAND) 

0.5 

0.06  mm/hr 

MAXIMUM  ALLOWABLE  ERRORS 

CALIBRATION 

RETREIVAL 

ALGORITHM 

ERROR 

SPECIFICATION 

WIND  SPEED 

0.9 

2.0  m/s 

WATER  VAPOR 

2.6 

2.0  kg/m1 

CLOUD  WATER 

4.7 

0.1  kg/m1 

RAIN  (OCEAN) 

38.1 

5.0  mm/hr 

ICE  (WINTER) 

6.0 

12.0  % 

ICE  (SUMMER) 

4.0 

12.0  % 

ICE  (37  GHz) 

10.7 

12.0  % 

RAIN  (LAND) 

38.5 

5.0  mm/hr 

GEOIXJCATION 


•  50  Cases  studied 

-  25  Ascending 

-  25  Descending 

•  19  Regions 

•  Revolutions  164  -  1437  (12/10/91  -  3/9/92) 

•  Within  Specification  (±7  Km) 
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PITCH  COPMCCTON 


KMMUnON  «MU 


ROa  CORRfCTION 
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CALIBRATION/VALIDATION  OF  THE  DMSP  MICROWAVE  WATER 

VAPOR  SOUNDER  (SSM/T-2) 

V.  Falcone,  J.  Morrissey,  M.  Griffin 

Phillips  Laboratory 
Geophysics  Directorate  (GPAS) 

Hanscom  AFB,  MA  01731 

D.J.  Boucher  and  B.H.  Thomas 

The  Aerospace  Corporation 
Defense  Meteorological  Satellite  Program 
PO  Box  92957 
Los  Angeles,  CA  90009 

R.G.  Isaacs  and  J.  Pickle 

Atmospheric  and  Environmental  Research,  Inc.  (AER) 

840  Memorial  Drive,  Cambridge,  MA  02139 


The  SSM/T-2  sounder  was  launched  late  in  1991  and  has  the  potential  to  provide 
global  water  vapor  profiles  for  a  variety  of  meteorological  applications.  We  describe  early 
results  to  calibrate  the  SSM/T-2  microwave  moisture  sensor  based  on  aircraft  underflights 
and  the  use  of  radiative  transfer  models.  Imagery  from  the  DMSP  Operational  Linescan 
System  (OLS)  is  used  to  help  identify  heavy  cloud  areas  which  may  impact  on  water  vapor 
sounding  accuracies.  Also  discussed  are  analyses  of  retrieved  SSM/T-2  moisture 
soundings  and  comparisons  to  upper  air  soundings. 
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CONFERENCE  ON  DMSP  RETRIEVAL  PRODUCTS 


Calibration/Validation  of  the 
DMSP  Microwave  Water  Vapor  Sounder  (SSM/T-2) 

Vincent  J.  Falcone,  Jim  Morrissey,  Michael  K.  Griffin 

PL/GP  Atmospheric  Sciences  Division 

John  Pickle,  Ron  Isaacs 

AER,  Inc. 

Bruce  Thomas,  Don  Boucher 

Aerospace  Corporation 


14  April  1992 


■  •ijmt.oesi' 


14-13  April  1992 


h 

DMSP  Retrieval  Products  Slide  No  I 


Calibration  and  Validation  of  the  SSM/T-2 


Objectives 

Obtain  Simultaneous,  Independent  Measurements  of 
Microwave  Radiance,  Water  Vapor  Mass  and  Humidity 
to: 


w 


1 .  Calibrate/Verify  SSM/T-2  Brightness  Temperatures  -5 
Over  Known  Atmospheric  and  Surface  Conditions  f- 

2.  Validate  the  AFGWC  EDRs  vs.  Rawinsonde  Data  r 

l 


14-15  Apnl  1992  DMSP  Retrieval  Products  Slide  No.  2 
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Calibration  and  Validation  Efforts 


•  PHILLIPS  LABORATORY  -  Calibration 

Absolute  Calibration  -  Underflights  of  Satellite 
Relative  Calibration  -  Radiance  Calculations 


•  AEROSPACE  CORPORATION  -  Validation 

Compile  and  Analyze  International  Radiosonde 
Network  Database  Plus  Supplemental 
Radiosonde  Launches 


14-1*  April  1992 


DMSP  Retrieval  Products  Slide  No.  3 


Calibration  Efforts 


Comparison  of  Satellite  Measured  Radiance  With 
Aircraft  Measured  Radiance: 

NASA  ER-2  Aircraft  With  a  SSM/T-2  Lookalike 
Will  Underfly  F-11 

Comparison  of  Satellite  Measured  Radiance  With 
Calculated  Radiance: 

Forward  Radiative  Transfer  Calculations  Using 
Radiosonde  Data 


14- IS  April  1992 


DMSP  Retrieval  Products  Slide  No  4 
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Calibration  Participants 


PL/GPA 

AER 

NASA 


-  Principal  Investigators 

-  Analysis 

-  ER-2  Flight  Measurements 

West  Coast  U.S.  May  '92 
East  Coast  U.S.  July  '92 


AEROSPACE  -  SSM/T-2  Data  Collection 


14-15  April  1992 


sSmXiiH:.'  *> 

DMSP  Retrieval  Products  Slide  No  5 


Millimeter-Wave  Imaging  Radiometer  (MIR) 


•  Channels:  89,  150,  183.31+1,  ±3,  ±7, 

220  and  325  GHz 

•  Scan:  Cross-track  (  50  samples  in  2.9  Seconds) 

•  Calibration:  External  to  Antenna 

•  Main  Beam:  3.5°  FOV  (1 .2  km  Footprint) 

•  Total  Power  Radiometer 


14-15  April  1992 


DMSP  Retrieval  Products  Slide  No.  6 
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DMSP  Retrieval  Products  Slide 


SSM/T 


PERFORMANCE  OF  DMSP  SPECIAL  SENSOR  MICROWAVE  HUMIDITY 
PROFILER  (SSM/T-2):  PRELIMINARY  RESULTS 

Donald  J.  Boucher  and  Bruce  H.  Thomas 


The  Aerospace  Corporation 
Defense  Meteorological  Satellite  Program 
PO  Box  92957 
Los  Angeles,  CA  90009 


The  SSM/T-2  is  a  cross-track  scanning  microwave  radiometer  operating  in  the 
frequency  range  from  90  to  183  GHz.  The  mission  of  the  SSM/T-2  is  to  provide  high 
quality  vertical  moisture  profiles  that  can  be  made  available  to  a  variety  of  models  running 
operationally  at  the  Air  Force  Global  Weather  Central  (AFGWC).  The  algorithm  currently 
used  at  AFGWC  to  extract  the  moisture  profiles  from  the  radiances  is  multiple  linear 
regression.  Historically,  regression  coefficients  have  been  organized  using  simple  static 
discriminates  (i.e.  latitude  and  longitude  as  for  the  SSM/T-1).  The  SSM/T-2  breaks  with 
this  transition  and  attempts  to  build  a  "dynamically  based"  set  of  piece  wise  linear 
regression  coefficients  stratified  by  air  mass.  A  brief  discussion  on  the  choice  of 
discriminates  for  air  mass  typing  will  be  presented.  In  addition,  Aerospace  is  in  the 
process  of  validating  the  SSM/T-2  derived  moisture  profiles  (against  synoptic  time  and 
special  radiosondes)  and  will  present  preliminary  retrieval  results. 
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Performance  of  the  DMSP  Special  Sensor  Microwave  Humidity  Profiler 
(SSM/T-2):  Preliminary  Results 


by 

Donald  J.  Boucher 
Bruce  H,  Thomas 

The  Aerospace  Corporation 
Defense  Meteorological  Satellite  Program 

Vincent  J.  Falcone 

Phillips  Laboratories 
Satellite  Meteorology 


^OverviewJ 

Calibration/Validation  data  collection  efforts 
SSM/T-2  hardware  performance  on  orbit 

Aerojet  Electro  Systems  (AES)  moisture  profile  retrieval  approach 

Moisture  retrieval  analysis 

Summary 


SSM/T-2  Calibration/Validation  Datasave 


1  -  —  l 

Data  collection  has  been  accomplished  by  Aerospace  Corporation  personnel 
assigned  to  HQ/AFGWC,  Offutt  AFB,  Nebraska 

Initial  version  of  the  SSM/T-2  ground  processing  software  used  to  support 
Cal/Val 

Software  errors  have  impacted  segments  of  the  SSM/T-2  data  collected  to 
date;  data  are  recoverable 

Initial  dataset 

30  Jan  -  8  Mar  1992 

Additional  datasaves  as  required;  field  experiments,  seasonal  samples,  etc. 
Data  types  collected 

DMSP  OLS  &  NOAA  AVHRR  visual  and  infrared  imagery  (smooth,  3nm); 
AFGWC'S  Satellite  Global  DataBase  (SGDB),  poiar  stereographic 

F-1 1  DMSP  special  sensor;  SSM/I,  SSM/T-1 ,  and  SSM/T-2 


Raw  sensor  data  has  been  archived  for  recall  as  necessary 


SSM/T-2^CaMbration/^ 


Sensor  Data  Records  (SDRs);  earth  located,  antenna  pattern  corrected 
brightness  temperatures 

Environmental  Data  Records  (EDRs);  temperature  and  moisture  profiles, 
surface  winds,  rain  rates,  etc. 


Conventional  observations;  WMO  radiosondes  worldwide,  synoptic  and 
asynoptic 

AFGWC  model  data;  every  3  or  6  hrs  where  available 

High  Resolution  Analysis  System  (HiRAS) 

Analysis  fields;  temperatures,  heights,  U-V  winds,  and  relative  humidity 
spanning  the  surface  to  10  mb 

Global  Spectra  Model  (GSM) 

0  hr  forecast  fields  (normal  modes  initialization);  temperatures,  heights, 
U-V  winds,  and  relative  humidity  spanning  1000-10  mb 

6  hr  forecast  fields;  temperatures,  heights,  U-V  winds,  and  relative 
humidity  spanning  1000-10  mb 


Real-Time  NEPHanalysis  (RTNEPH) 


COLO  LOAO  COUNTS 


SSM/T-2  Hardware  Performance  on  Orbit 


All  channels  meet  or  exceed  the  required  specifications  (.5  Kelvins 
NEDT) 

First  1600  satellite  revolutions  have  demonstrated  the  sensor's  stability 

Revolution  1655  summary  {units  =  Kelvins) 

Channel  1  (183  +1-  3  GHz)  NEDT  =  .34 

Channel  2  (183  +/-  1  GHz)  NEDT  =  .25 

Channel  3  (183  +/-  7  GHz)  NEDT  =  .46 

Channel  4  (91.655  GHz)  NEDT  =  .27 

Channel  5  (150.00  GHz)  NEDT  =  .33 

The  absolute  calibration  of  the  sensor  will  be  determined  by  the  Phillips 
Laboratory  component  of  the  CalA/al  team  by;  utilizing  radiative  transfer 
modeling,  aircraft  flights  with  similar  microwave  instruments,  etc. 

Details  were  presented  in  the  previous  briefing! 


15001 


1000  b 


500 


u  I 
-  1000 


CHAN5. 


CHAN2- 


CHAN3* 
CHAN  1. 


F-1  1  SSMT/2  REV 

■  r-  y— i — i — i — i — i — r— « — r  y  j 


1655 


SAMPLE  NUMBER  (A  samples/scan) 


"l — I — 1 — i — ! — r~ 


.CHAN  5 


.CHAN2 


CHAN4., 

_ _ CHAN  4 

CHAN  1  ORBIT  MEAN  - 

742.704 

ORBIT  STD  DEV  = 

7  07012 

C-'AN  2  0R3IT  MEAN  = 

1 100.04 

ORBIT  STD  DEV  - 

7.27648 

CHAN  3  OR  Bit  MEAN  = 

859.463 

ORBIT  STD  DEV  « 

22  8283 

CHAN  4  ORBIT  MEAN  = 

457.246 

ORBIT  STD  DEV  - 

23  7463 

CHAN  5  ORBIT  MEAN  = 

1249.84 

ORBIT  STD  DEV  = 

33  9327 

1 

,  .  1  ,  ,  .  ■  , 

,  ,  .  1 . 

^  A*rtM 

0 

1000 

2000 

3000 

4000 


Mor/coLD  oirrcRtNrt  (counis)  hoi  load  counts 


SSMT/2  REV  = 


CHAN  5, 
CHAN2, 


3200  b 


I.CHAN3 


CHAN  4 


i 

CHAN  1  0R8IT  MEAN  = 

3473.86 

ORBIT  STD  OEV  - 

6.80334 

CHAN  2  ORBIT  MEAN  = 

3551.76 

ORBIT  STD  DEV  = 

7.15602 

! 

CHAN  3  ORBIT  MEAN  = 

3334.95 

ORBIT  STD  OEV  > 

22.9619 

CHAN  4  ORBIT  MEAN  « 

3243.64 

ORBIT  STD  OEV  = 

24.1056 

CHAN  5  ORBIT  MEAN  = 

3529.40 

ORBIT  STD  DEV  = 

36.4263 

1000  2000 
SAMPLE  NUMBER  (4  samples/scan) 

F-1  1  SSMT2  REV  =  1655 


CHAN 4, 
CHAN  1 . 


.CHAN  4 
JOHAN  1 


2600  ,~ 


CHAN3 

CHAN2'* 


2000 


OPBIT  VARIANCE  oe  The  OIEEERENCE 
OPB'T  VARIANCE  OE  the  OIEEERENCE 
OPQiT  VARIANCE  OE  THE  OIEEERENCE 
OPBIT  VARIANCE  OE  THE  OIEEERENCE 
ORBIT  VARIANCE  OE  THE  OIEEERENCE 


3.31974  COUNTS/DEG  K  -  9.28964 

1.65445  COUNTS/DEC  K  -  8.33920 

4.38223  COUNTS/DEG  K  -  8.42001 

1  42932  COUNTS/DEG  K  -  9.47753 

3.81401  COUNTS/DEG  K  -  7.75360 


0  1000  2000  3000 

NUMBER  OE  SAMPLES  (4  somples/scar ) 
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HOT  COUNT  STD  DEV 


F-1  1  REV  =  1655 


CHAN3 


CHAN  1 


CHAN  1  ORBIT  MEAN  COUNTS/DEC  K 
CHAN  2  ORBIT  MEAN  COUNTS/DEC  K 
CHAN  3  ORBIT  MEAN  ■'  ;NTS/DEG  K 
CHAN  4  ORBIT  MEAN  COUNTS/DEG  K 
CHAN  5  ORBIT  MEAN  COUNTS/DEC  K 


11.8122 
12.0772 
1 1.3391 
1 1 .0285 
12.0000 


NUMBER  OF  DISCRETE  8  SCAN  CROUPS 


Aerojet  SSM/T-2  Retrieval  Approach  | 

Statistical  Approach  (multiple  linear  regression,  D-matrix) 

Piece-wise  linear  approximation  to  non-linear  water  vapor  distribution 
Air  Mass  Stratification  (AMS) 

Atmosphere  Index  is  determined  by: 

1)  Surface  type  from  eighth  mesh  database  (ocean,  land,  ice,  and 
coast) 

2)  Estimating  the  total  water  vapor  mass  within  the  column 

3)  Estimated  high  altitude  water  vapor  mass:  using  terrain  height 
adjusted  brightness  temperatures  if  necessary 

4)  SSM/T-1  brightness  temperatures  of  lower  troposphere 

Ocean:  SSM/T-1  53.2  GHz  brightness  temperature  threshold 

Land:  Terrain  height  adjusted  SSM/T-1  50.5  and  53.2  GHz 
brightness  temperatures  and  estimation  of  "high  altitude"  water 
vapor  mass 

Ice:  SSM/T-1  53.2  GHz  temperature  threshold 


Aerojet  SSM/T-2  Retrieval  Approach  Cont 


Coast:  Terrain  height  adjusted  SSM/T-1  53.2  GHz  temperature 
threshold  **  Note,  SSM/T-2  91  and  150  GHz  and  SSM/T-1  50.5 
GHz  not  used  due  to  high  variations  of  surface  emissivity 

Once  the  "air  mass  index"  and  terrain  height  is  known,  the  proper 
D-matrix  partition  can  be  selected 

The  air  mass  stratification  approach  removes  seasonal  and  latitude 
zone  dependance  when  compared  to  the  SSM/T-1  D-matrix  partitions 

**  Note:  An  air  mass  stratification  scheme  for  the  SSM/T-1  is 
planned  for  replacing  the  seasonal  and  latitude  dependant 
D-matrix  mechanism  at  AFGWC  (1992-1993) 


Preliminary  Moisture  Retrieval  Performance  Analysis 


Early  focus  has  been  placed  on  relative  humidity  retrieval  accuracy 

Collocations  of  SSM/T-2  EDRs  to  RAOBs  have  provided 
approximately  470  matches/orbit  within  1  hour  and  100km  (ideal  orbit, 
12z  +/-  50  minutes) 

1 47  coastal  matches  from  22  RAOBs 

1  ice  match  from  1  RAOB 

22  ocean  matches  from  4  RAOBs 

301  land  matches  from  36  RAOBs 

Retrieva'  statistics  not  yet  available 

SSM/T-2  ground  processing  software  errors  have  prevented  analysis  of 
non-ocean  atmospheres,  below  500  mb,  because  data  are  stored  as 
missing.  The  problem  is  under  investigation  and  will  be  corrected 
ASAP 

Ocean  retrievals  appear  to  be  very  sensitive  to  SSM/T-1  lower 
atmospheric  channel  data  (as  expected),  however,  resultant  EDR 
imagery  shows  the  non-Dhvsical  gradients/hatching  based  upon 
SSM/T-1  scan  geometry  and  resolution  (/. 4] 


Summary^ 


Data  collection  efforts  have  provided  an  initial  dataset  for  the  Cal/Val 
analysis  with  provisions  for  additional  samples  as  necessary 

Initial  SSM/T-2  ground  processing  software  contains  errors  which 
have  impacted  the  data  quality  supplied  to  the  Cal/Val  team; 
corrections  pending  and  data  are  recoverable 

SSM/T-2  hardware  preforming  well  on  orbit;  well  within  contract 
specifications 

Absolute  calibration,  TBD  by  Phillips  Labs,  may  impact  validation  work 
and  the  corresponding  accuracy  assessment 

Aerojet's  "piece  wise  multiple-linear  regression"  mechanism  provides 
interesting  approach  to  water  vapor  retrievals;  remains  under 
investigation  for  accuracy 
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DMSP  BLOCK  6:  A  DATA  DRIVEN  APPROACH 


M.W.  Borden,  F.  Kelly  and  G.  Mandt 

HQ,  Air  Force  Space  Systems  Division 
(SSD/MWB) 

PO  Box  92960 

Los  Angeles  AFB,  CA  90009-92960 

D.  L.  Glackin  and  J.  Bohlson 

The  Aerospace  Corporation 
Defense  Meteorological  Satellite  Program 
PO  Box  92957 
Los  Angeles,  CA  90009 


The  DMSP  Block  6  program  features  an  aggressive  and  innovative  data  driven 
strategy  for  the  acquisition  of  the  next  generation  DoD  weather  satellite.  The  Block  6 
program  maintains,  as  a  baseline,  the  operational  utility  of  the  DoD  weather  satellite  system 
while  seeking  to  optimize  operational  efficiency  within  a  specified  life  cycle  cost  cap.  It 
features  the  development  of  a  complete  meteorological  satellite  system,  embracing  both 
space  based  and  ground  based  assets.  Particular  attention  is  focused  on  the  insertion  of 
maturing  technologies  into  the  developing  Block  6  effort.  The  system  design  concept 
responds  to  the  data  needs  of  the  operational  DoD  User.  System  efficiencies  are 
engendered  through  the  use  of  data  fusion,  improved  mass  memory,  sensors,  power 
concepts  and  the  application  of  Knowledge  Based  Systems  (KBS)  to  optimize  User 
efficiency. 
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DMSP  BLOCK  6: 

A  DATA  DRIVEN  APPROACH 


PRESENTED  BY: 

OPT  MARK  W.  BORDEN 


DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM 

UNCUHnD 


UNCLASVXD 

OVERVIEW 


•  DMSP  MISSION 

•  DMSP  BLOCK  6 

•  BLOCK  6  REQUIREMENTS 

•  BLOCK  6  CAPABILITIES 

•  BLOCK  6  STATUS 


DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM 

UWCUiWO 
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UNCLAMFD 


DMSP  MISSION 


COLLECT 

I  VISIBLE  AND  IR  CLOUD  DATA 

l  OTHER  METEOROLOGICAL,  TERRESTRIAL, 
OCEANOGRAPHIC  AND  SOLAR-GEOPHYSICAL  DAT) 

DISSEMINATE 

>  GLOBAL  DATA  TO  STRATEGIC  USERS 

>  REAL-TIME  DATA  TO  TACTICAL  USERS 


DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM 

UNC IAMIFWD 


gNcusms 


DMSP  BLOCK  6 


•  EVOLUTIONARY  DEVELOPMENT  OF  DMSP  BLOCK  5D-3 

•  CONTINUE  TO  PROVIDE  50-3  CAPABILITY 

•  ATMOSPHERE,  OCEANS,  LAND,  SOLAR-GEOPHYSICAL 

•  TRI-SERVICES 

•  ADD  LOW  COST  ANO  IMPACT,  HIGH  YIELD  IMPROVEMENTS 

•  WITHIN  COST  CAP  FOR  DEVELOPMENT  +  10  YEARS 
OPERATION 

•  ENHANCED  CAPABILITIES  VIA  PRICED  OPTIONS 

•  FY  2005  INITIAL  OPERATIONAL  CAPABILITY 

•  CONCEPT  STUDY  COMPLETED  (JAN  88  -  APR  91)  -  4 
CONTRACTORS 

•  RISK  REDUCTION  PHASE  IN  PROGRESS  -  2  CONTRACTORS 

•  SYSTEMS  QUAL/INITIAL  PRODUCTION  FY  98 


DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM 
woLAnno 


UNdASSFCD 


BLOCK  6  REQUIREMENTS 


•  BASELINE  5D-3  CAPABILITY 

•  IMPROVEMENTS  FROM  JCS  ENVIRONMENTAL  DATA  LIST 

•  VISIBILITY,  SEA  SURFACE  TEMP,  IMAGER  RESOLUTION... 

•  INDIVIDUAL  AGENCY  INPUTS 

•  AIR  >  FATHER  SERVICE 

•  SPACE  COMMAND 

•  NAVY 

•  ARMY 

•  DEPARTMENT  OF  COMMERCE  (POSSIBLE  REQUIREMENTS) 


DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM 


tMOASSFCb 

BASELINE  SENSOR  COMPLEMENT 


•  ELECTRO-OPTICAL  IMAGING  SUITE 

7  BANDS,  0.4  -  12.5  um 

•  MICROWAVE  IMAGER  SOUNDER  SUITE 

•  IMAGER  £  7  CHANNELS,  4  WAVaENGTHS,  DUAL 
POLARIZATION.  19-91  GHz 

•  TEMPERATURE  SOUNDER  TBD  CHANNaS,  50  -  60  GHz 

•  MOISTURE  SOUNDER  4  CHANNaS,  150  -  183  GHz 

•  SEPARATE  OR  COMBINED  SCANS 

•  SOLAR  GEOPHYSICAL  SENSORS 
•TBD 

•  NAVY  OPTION 

•  ALTIMETER  AND  SCATTEROMETER  POSSIBLE 


DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM 


UNOAfMD 

DATA  PARAMETER  CORRELATION 

PARAMETER/PROOUCT  SENSINO  TECHNOLOGY 

CLOUO  IMAGERY 
VMWUTY 

LAND  SURFACE  TEMPERATURE 
SEA  SURFACE  TEMPERATURE 
SURFACE  ALBEDO 
PROCESSED  CLOUO  DATA 
VCOTATTON  I NOCX/ SURFACE  TYPE 
PR£OP(TATTON 
TOTAL  WATER  CONTENT 
SURFACE  MOISTURE 
SNOW  *  LANDLOCKED  ICE  COVER 
SEA  ICC 

OCEAN  SURFACE  WINO  SPEED 


DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM 

UNCLAttfYD 


-  - - — 

V1SIBLE/INFRARE0 

IMAOINO 

- 

r — £ 

■ 

- 1 

MICROWAVE 

IMAOINO 

mmrnm  — . 

il  H 

** 

UMCL WfD 

DATA  PARAMETER  CORRELATION 


PARAMETER/PROOUCT  SENSINC  TECHNOLOGY 


DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM 

wcli— ci 
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TECHNOLOGY  DEVELOPMENT  IMPACTS 


DATA  FUSION 

•  SYNERGISM  OF  MICROWAVE,  IR  AND  VISIBLE  DATA 

•  ONGOING  PROJECTS  THROUGH  AF  SBIR 
MASS  MEMORY 

•  TAPE  RECORDERS  LIMITING  TODAY  -  MORE  SO  IN  BLOCK  6 


•  POTENTIAL  MASS  MEMORY  UPGRADES 

•  SOLID  STATE  MEMORY  DEVICES 

•  REWRITABLE  OPTICAL  DISK  DEVICES 

•  SENSORS 

•  ELLSE 

•  LIDAR 

•  HIGH  ANGLE  MICROWAVE 

•  POWER 

•  PAS P  PLUS 

•  NIH  BATTERIES 

•  KNOWLEDGE  BASED  SYSTEMS 

•  SAGE 

•  ANOMALY  AND  STATION  KEEPING  ANALYSIS 


DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM 


^  UNCLASSMS 

i  KEY  SENSOR  RELATED 
^  DEMONSTRATIONS/EXPERIMENTS 


•  OMIS  ENGINEERING  DEVELOPMENT  UNIT 

•  BENCH  TEST  UNIT 

•  INCORPORATE  CHANNELS  1  -  6 

•  EXPERIMENTAL  LOW  LIGHT  SENSOR  (ELLSE) 

•  OMIS  CHANNEL  6  (NIGHT  CLOUD  *  AURORAL 
IMAGERY) 

•  IMPROVED  RESOLUTION  *  LOW  RADIANCE  GOALS 

•  OPTIONAL  SPACE  FLIGHT  TEST 

•  HIGH  EARTH  INCIDENCE  ANGLE  MICROWAVE  SENSING 

•  NORTH  SEA  TOWER  TESTS  (NRL) 

•  SURFACE  EMISSIVITY  STUDIES/MODELING  (PL) 

•  MAY  GIVE  CONTIGUOUS  MISS  DATA 
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CONTRACT  STRUCTURE 


DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM 

UNCLASSns 
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Session  2:  Operational  and  Quasi-operational  Retrieval  Approaches  and  Products 

(McC  late  hey,  chair) 

This  session  dealt  with  the  interpretation  and  dissemination  of  current  DMSP 
microwave  products,  the  algorithms  used  to  generate  them  and  the  system  used  to  distribute 
the  results.  There  also  was  considerable  focus  on  research  efforts  being  conducted  to 
determine  the  accuracy  of  current  operational  algorithms  and  to  evaluate  new  and  different 
algorithms.  There  were  two  papers  on  the  SSM/T-1,  the  DMSP  temperature  sounder,  and 
ten  on  the  SSM/I,  the  DMSP  microwave  imager.  The  paper  presented  by  Capt.  T.  Neu  of 
the  Air  Force  Global  Weather  Central  described  the  processing  of  data  from  these  two 
systems  as  well  as  the  new  mc  uure  sounder,  SSM/T-2  and  the  future  SSMIS  which  will 
integrate  the  DMSP  sounder  and  imager  into  a  single  system. 

The  two  papers  on  SSM/T-1  dealt  with  the  shared  processing  network  and 
products.  Temperature  products  are  being  made  available,  but  are  being  used  by  NMC  as 
backup  to  TOVS.  T.  Reale  of  NOAA/NESDIS  suggested  that  the  DMSP  accuracy  appears 
to  be  worse  than  TOVS,  although  he  indicated  that  the  change  from  seasonal  to  weekly 
updating  of  the  statistical  data  base  used  in  the  retrieval  process  has  improved  the  results 
somewhat.  In  a  comment  made  from  the  audience,  H.  Woolf  of  NOAA/NESDIS 
(Madison,  Wisconsin)  expressed  the  view  that  the  DMSP  results  are  not  all  that  much 
different  than  TOVS.  Perhaps  we  need  to  determine  a  more  satisfactory  basis  of 
comparison  for  satellite  retrieved  products  which  is  more  consistent  with  the  nature  of  the 
sensing  technique. 

The  future  availability  of  data  from  the  Microwave  Water  Vapor  Profiler  (SSM/T-2) 
discussed  by  Capt.  Neu  was  met  with  substanti  J  interest.  Of  course,  the  determination  of 
water  vapor  profiles  will  depend  strongly  on  the  temperature  profiles  determined  from  the 
SSM/T-1. 

The  remaining  ten  papers  all  related  to  SSM/I  data  processing,  data  distribution  and 
alternative  algorithm  design  and  testing.  M.  Colton  described  how  the  Navy's  Fleet 
Numerical  Oceanographic  Center  (FNOC)  applies  the  Hughs  algorithm  and  distributes  the 
results  over  the  Shared  Processing  Network  (SPN).  J.  Fiore  of  SM  Systems  and  Research 
Corp.  described  a  testbed  that  has  been  developed  for  the  same  data  base.  He  emphasized 
the  comparison  of  total  precipitable  water  algorithms.  J.  Alishouse  described  alternative 
algorithms  for  Total  Precipitable  Water  and  die  need  for  further  validation  and  quality 
control.  N.  Grodv  and  R.  Ferraro  are  evaluating  precipitation  and  surface  products. 

W.  Pichel  is  evaluating  ice  and  wind  products  by  comparison  with  ice  analyses  produced 
using  visible  and  infrared  imagery.  And  finally,  the  NASA-sponsored  project,  WETNET, 
was  described  by  Jim  Dodge  of  NASA  Headquarters,  standing  in  for  M.  Goodman  of 
NASA  Marshall  Space  Flight  Center  (MSFC).  WETNET  is  a  data  distribution  system 
offering  people  in  the  R&D  community  an  opportunity  to  work  with  the  necessary  data  sets 
for  algorithm  improvement  with  respect  to  SSM/I  and  eventually  also  SSM/T-1  and 
SSM/T-2. 


SESSION  2:  OPERATIONAL  AND  QUASI-OPERATIONAL  RETRIEVAL 
APPROACHES  AND  PRODUCTS 


NOAA/NESDIS  OPERATIONAL  PROCESSING  OF  SSM/T  DATA 

F.llen  L.  Burdsall 

NOAA/NESDIS 
Washington,  DC  20233 


Since  1987,  the  National  Oceanic  and  Atmospheric  Administration  (NOAA)/ 
National  Environmental  Satellite  Data  and  Information  Service  (NESDIS)  has  been 
operationally  producing  vertical  temperature  profiles  (soundings)  from  the  Special  Sensor 
Microwave/Temperature  (SSM/T)  instruments  on  the  Defense  Meteorological  Satellite 
Program  (DMSP)  satellites.  In  the  years  since  1987,  the  NOAA/NESDIS  SSM/T 
processing  system  has  been  enhanced  and  modified  to  improve  the  quality  of  soundings 
from  the  SSM/T  instrument.  A  brief  history  of  the  SSM/T  processing  system  will  be 
presented.  An  overview  of  the  current  system  will  be  discussed  and  data  distribution  and 
archive  will  be  addressed  in  detail. 
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BRIEF  HISTORY 
OF  NESDIS 

SSM/T  SOUNDING  PROCESSING 


1982  NESDIS  RECEIVES  SSM/T  SOFTWARE 
DEVELOPED  BY  AEROJET 

1985  AIR  FORCE  SENDS  SSM/T  DATA  VIA 
LAND  LINE  (DATA  ARCHIVED) 

1986  SYNTHETIC  REGRESSION  RETRIEVAL 
METHOD  REPLACED  WITH  SIMPLE 
MULTIPLE  LINEAR  REGRESSION 

1986  AIR  FORCE  SENDS  SSM/T  DATA  VIA 
-1987  SHARED  PROCESSING  LINK 


1987  JUNE  -  SSM/T  OPERATIONAL 

PROCESSING  AT  NESDIS 

1989  IMPLEMENT  PRECIPITATION 

SCREENING  ALGORITHM  DEVELOPMENT 
BY  DR.  GRODY 

1 991  LIMB  CORRECTION  IMPLEMENTED 
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SSM/T  PROCESSING  AT  NESDIS 


DPSS/CEMSCS 


HDS 


NMC 


DATA  - 

LAYER  MEAN  TEMPERATURE 
BRIGHTNESS  TEMPERATURE 


DISTRIBUTION  - 

DISK  FILE  SHARED  BY  NMC 
AND  NESDIS 


USE  - 

BACKUP  TO  TOVS 

4* 

ECMWF 


DATA  - 

LAYER  MEAN  TEMPERATURE 
BRIGHTNESS  TEMPERATURE 

DISTRIBUTION  - 

RECEIVED  IN  REALTIME  FROM 
NMC 

USE  - 

ANALYSIS  AND  BACKUP  TO  TOVS 
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NAVY  -  FNOC 


DATA  - 

LAYER  TEMPERTURE 
BRIGHTNESS  TEMPERATURE 
GEOPOTENTIAL  HEIGHTS 

DISTRIBUTION  - 

TEMPERTAURE  VALUES  RECEIVED  IN 
REALTIME  FROM  NMC 

GEOPOTENTIAL  HEIGHTS  RECEIVED  IN 
REALTIME  FROM  NESDIS 


USE  - 


USED  IN  OPERATIONAL  PROCESSING 

i  miiiihf ^ 


ARCHIVE 


DATA  - 


LAYER  TEMPERTURE 
LEVEL  TEMPERATURE 
BRIGHTNESS  TEMPERATURE 


PROCESS  - 

IB  DATA  SETS  ARCHIVED  DAILY 
PROC'JCT  ARCHIVE  PERFORMED  WEEKLY 


STORAGE  - 

SATELLITE  DATA  SERVICES  DIVISION  OF  THE 
NATIONAL  CLIMATIC  DATA  CENTER  SERVICES 
ALL  USER  REQEUSTS  FOR  SSM/T  DATA 
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NOAA/NESDIS  SSM/T  ALGORITHMS  AND  EVALUATION 

Tony  L.  Reale  and  David  Donahue 

NOAA/NESDIS 
Washington,  DC  20233 


Defense  Meteorological  Satellite  Program  (DMSP)  sounding  products  are  currently 
produced  operationally  by  NESDIS  and  distributed  to  national  and  international  users. 

Two  sun  synchronous  DMSP  polar  orbiting  satellites,  each  containing  the  seven  channel 
Special  Sensor  Microwave/Temperature  (SSM/T).  sounder  are  being  flown  operationally. 
This  provides  an  average  of  four  independent  measurements  daily  for  most  locations  over 
the  earth.  The  temperature  soundings  produced  are  from  the  surface  to  about  10  mb  (30 
km)  with  a  nominal  horizontal  resolution  of  approximately  175  km.  An  overview  of  the 
current  scientific  algorithms  to  compute  atmospheric  sounding  products  from  the  SSM/T 
sounder  measurements  is  presented.  Results  showing  data  at  various  stages  of  the  radiance 
and  sounding  computation  process  are  displayed.  An  evaluation  of  the  sounding  products 
accuracy  based  on  vertical  accuracy  statistics  relative  to  collocated  radiosondes,  and 
horizontal  temperature  fields  relative  to  other  satellite  and  conventional  meteorological  data 
is  presented.  Concluding  remarks  on  future  plans  to  upgrade  the  DMSP  operational 
sounding  products  are  provided. 
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NOAA  NESDIS  OPERATIONAL  SSM/T 
SOUNOING  ALGORITHM  AND  METEOROLOGICAL 
EVALUATION 


Tony  Reale 
David  Donahue 
Pamela  Taylor 


NOAA/NESOIS/ORA 

Satellite  Applications  Lab 


OUTLINE 

o  Operational  System 

o  Science  Algorithms  and  Applications 

o  Sounding  Accuracy 

o  Mcteoro logical  Consistency 
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o  Preprocess 


o  Terrain/Surface  Correction 


o  Precipitation  Filter 


o  Statistical  Regression  Retrieval 
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SSM/T  PROCESSING 


DAILY 


OFF-LINE 


WEEKLY 


N 


Pi  ■  Fj  -  BT,) 


where:  I  =  Channel 

j  =  Pressure  Level 
Pj  =  Atmospheric 
Parameter 

pj  =  Mean  Atmospheric 
Parameter 


N  =  Number  of  Channels 

C 1 , j  =  Coefficients 

BTI  =  Brightness  Temperature 

BT 1  =  Mean  Brightness 
Temperature 


Present  Near-Term  Future 

OtySP  System  OMSP  System  DMSP  Sytem 

(1992)  (1993)  (Late  1990’s) 


REGRESSION 


PHYSICAL/ 


IFDMSP 


NESDIS/PSB  S3E  VERTICAL  DISPLAY 
3on  -  jos  no  lno/sla  DAr/wa 
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SMPLE  SIZE 


CONCLUSIONS 


o  OPERATIONAL  SOUNDINGS  USING  A  STATISTICAL  APPROACH 

o  ACCURACY  COMPATIBLE  TO  TOYS  CLOUDY  SOUNDINGS 

o  CONSISTENT  GLOBAL  HETOROLOGICAL  REPRESENTATION 

o  VITAL  PROOF-OF-CONCEPT  FOR  AHSU 
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DEFENSE  METEOROLOGICAL  SATELLITE  PROGRAM  MICROWAVE 
RADIOMETER  PROCESSING  AT  AIR  FORCE  GLOBAL  WEATHER 

CENTRAL 


T.J.  Neu 

HQ  AFGWC,  SYS/Satellite 
106  Peacekeeper  Drive,  STE  2N3 
Offutt  AFB,  NE  6811 3-4039 


The  Department  of  Defense  operates  polar  orbiting  meteorological  satellites  as  part 
of  the  Defense  Meteorological  Satellite  Program  (DMSP).  Along  with  the  primary  visual 
and  infrared  sensor  (OLS),  the  DMSP  satellites  carry  three  passive  microwave  radiometer 
sensors.  These  special  sensors  consist  of  the  Microwave  Imager  (SSM/I),  Microwave 
Temperature  Sounder  (SSM/T),  and  Microwave  Water  Vapor  Profiler  (SSM/T-2).  Air 
Force  Global  Weather  Central  (AFGWC)  at  Offutt  AFB  is  the  sole  Department  of  Defense 
location  where  all  of  the  ground  processing  for  these  special  sensors  is  performed.  A  brief 
review  of  DMSP  orbital  characteristics,  data  stream  transmission,  and  data  ingest  to 
AFGWC  computer  systems  will  be  presented.  The  general  description,  function,  and  uses 
of  each  microwave  radiometer  will  then  be  discussed.  Emphasis  will  be  placed  on 
operational  retrieval  approaches  and  products  generated  at  AFGWC.  An  overview  of  the 
Microwave  Imager  Sounder  (SSM/IS)  will  present  technical  enhancements  and  future 
integration  of  the  sensor  at  AFGWC. 
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DMSP  MICROWAVE  RADIOMI  i’EH 


PROCESSING!  AT  AFQWC 


presented  hy 

Ca.pt  Tliomaa  Neu 


HQ  AFGWC ,  SYS 

109  Peacekeeper  Dr.  STE  2N3 
Offutt  AFB ,  NE  68113-4039 
(402)  294-7444 


TOPICS  OE  DISSCUSSION 


1) 

DMSP 

h 

> 

0 

R  V  I 

E  W 

2) 

SS 

M  /  I 

3) 

SS 

M  /  T 

- 1 

4) 

SS 

M  /  T 

-  2 

5) 

US 

ERS 

O  E 

DA 

T  A 

6) 

SS 

M  /  I 

S 
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DMSP  SATELLITE  SENSOR  CONFIGURATION 


BLOCK 

5D-2 

UL0< 

IK  5 D  - 3 

SENSOR 

SATELLITE 

SB 

S10 

SI  l 

S  12 

S  13 

S  1 4 

S  1 5 

S  16 

thru 

S20 

FLIGHT  ID 

F9 

F0 

F 10 

_ 

FU 

_ 

OLS  Operational 
Linescan  System 

■ 

■ 

N 

B 

N 

X 

X 

X 

SSM/I  Microwave 

I mager 

■ 

DL 

N 

B 

N 

X 

X 

X 

SSM/T1  Microwave 
Temperature  Sounder 

DL 

■ 

N 

X 

N 

X 

X 

X 

SSM/T2  Microwave 
Water  Vapor  Profiler 

N 

X 

X 

X 

SSM/IS  Microwave 
Imager  Sounder 

_ 

_ 

X 

F  =  FAILURE  D  =  DEGRADED  L  =  LIMITED  N  =  NOMINAL  X  =  PLANNED  AFGWC 

1 )  Orbital  Character lgtics : 

Orbit  is  polar  sunsynchr onous 
Altitude  of  450  Nautical  miles  (833  KM) 

Period  of  101.4  Minutes 
Inclination  of  98.7  Degrees 

Local  Time  of  Ascending  Nodal  Crossing:  (LTAN) 

F-8  =  0611,  F-9  -  0930,  F- 1 0  =  1934  to  2112  »•  ,  F-ll  =  1700 

»»  F-10  has  an  elliptical  orbit  which  causes  the  nodal  crossing  time  to 
process  50  seconds  per  week.  Limited  use  after  2112  LTAN  in  March  1993. 

2 )  Data  Stream  Transmi ss i on : 

Satollito  recorder  to  downlink  site:  Pogo,  Boston,  Fairchild.  Hula 
Downlink  site  to  DOMSAT:  Atlantic  or  Pacific 
DOMSAT  to  1000  SOG  at  Offutt  AFB .  NE 
1000  SOG  to  AFGWC  Site  III 

Site  III  to  AFGWC  mainframe  system  5  (UNISYS  1100/90) 
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SPECIAL  SEMSOR  MICROWAVE  IMAGER  (SSM/I) 


1 )  Passively  collects  microwave  energy 

a)  Seven  channels:  19V,  19H,  22V,  37V,  3711.  85V.  and  85H  GHz 

b)  Antenna  and  sensor  rotate  at,  31.6  RPM 

c)  Conical  sea  view  of  earth  at  an  angle  of  53.1  degrees 

d)  Swath  width  is  1400  km,  angular  sector  of  102.4  degrees 

e)  Scan  period  is  1.9  seconds,  sub' satel 1 i te  point  travels  12.5  km 

e)  Resolution  is  25  Km  (12.5  Km  for  85  GHz). 

f)  Samples:  128  discrete  samples  per  scan  for  both  85  GHz  channels 

64  discrete  samples  on  alternate  scans  for  others 


FIGURE  3.0-3  ACTIVE  SCAN  GEOMETRY 


1 
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SPECIAL  SENSOR  MICROWAVE  IMAGER  (SSM/I)  Cont. 


2)  Calculates  Sensor  Data  Records  (SDR'fl):  Kelvin 

a)  Earth  Located:  Geodetic  1  a t i tude/ 1 ong i tude  coordinates,  oblate  spheroid, 

LaGrange  interpolation  algorithm. 

b)  Surface  Tagged:  La t 1 tude/ 1 ong i tude  look  up  table;  ocean,  land,  coast, 

and  ice 

c)  Calibrated:  Channel  radiometric  voltages,  cold  and  hot  temperature 

voltages,  linear  voltage  to  temperature  conversion 

d)  Antenna  Pattern  Corrected:  Sidelobe  contributions,  Channel  Selection 

Table 


3 )  Calculates  Environmental  Data  Records  (EDR's) 

a)  Earth  located  values  for  specific  environmental  parameters 

b)  Dependent  upon  surface  type 

-  Ice:  Ice  age,  Ice  Edge,  Ice  Concentration 

-  Land:  Rain  Rate,  Soil  Moisture,  Surface  Type,  Surface  Temperature 

Snow  Depth 

-  Ocean:  Surface  Wind.  Rain  Contamination,  Rain  Rate,  Cloud  Water, 

Water  Vapor 

-  Coast:  None 

c)  Surface  Types:  Land,  Vegetation,  Arctic,  Ice,  Possible  ice,  Ocean, 

Coast,  Flooded  soil,  Dense  vegetation,  Land,  Dry  arabie 
soil,  Moist  arable  soil,  Desert,  Rain  over  vegetation, 
Rain  over  soil.  Composite  water  and  vegetation,  Wet 
soil,  Dry  snow,  Wet  snow,  Refrozen  snow 
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SPECIAL  SENSOR  MICROWAVE  TEMPERATURE  SOUNDER  (SSM/T-1) 


1 )  Paaaively  collects  microwave  energy 


a) 

Seven  Channels:  in  oxygen  absorption  band 

Channel  Peaking  Height  (km)  Frequency  (GHz) 

Bandwidth  (MHz) 

1 

0 

50.5 

400 

2 

2 

53.2 

400 

3 

6 

54 . 35 

400 

4 

10 

54 . 9 

400 

5 

30 

58.4 

115 

6 

16 

58.825 

400 

7 

22 

59.4 

250 

b) 

Scan  Type: 

Spin  scan,  cross-track 

nad  i  r 

c) 

Positions : 

Seven  cross-track,  two 
0,  +  /-  12.  +/-  24.  and 

cal  ibration 

+  /-  36 

(hot  and  cold) 

total  32  seconds,  2.7  second  dwell  time  cross-track 
e)  Resolution:  174  km,  185  km,  220  km,  and  304  km  (view  angle  dependent) 


Figure  D-l.  Q  -  0°,  Calm 
Sea  Background 


Figure  D-2.  0  -  4  1.67°,  Calm 

Sea  Background,  Horizon¬ 
tal  Polarization 
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SPECIAL  SEHSOR  MICROWAVE  TEMPERATURE  SOUWDER  (SSM/T-1)  Cont. 


2 )  Calculates  Sen3or  Data  Records  (SDR'fl)  :  Kelvin 

a)  Earth  Located:  Oblate  spheroid,  vector  Interpolation  to  obtain 

sub-satel 1 i te  latitud e/longitude. 

b)  Calibrated:  Dlcke  radiometer;  channel  radiometric  voltages,  cold  and 

hot  temperature  voltages  (10  scan  average),  linear  voltage 
to  temperature  conversion. 

c)  Antenna  Pattern  Corrected:  Antenna  temperatures  multiplied  by  antenna 

pattern  correction  coefficients. 

d)  Terrain  and  1000  mb  Heights:  Obtains  values  from  database. 


3)  Calculates  Environmental  Data  Records  (EDR's) 


a)  Precipitation  Check:  Discards  if  channel  1  liquid  precipitation 

threshold  is  exceeded  over  ocean. 

b)  Terrain  Height  Corrected:  Corrects  four  lower  channels  to  equivalent 

sea  level  brightness  temperatures  if  terrain 
height  exceeds  25  m. 

c)  Multiple  Regression  for  Fast  Inversion:  less  computational  intensive 

Linear  relation  between  a  parameter  vector  and  a.  set  of  measured  data, 
p  -  <p>  =  D  (  d  -  <d>  ) 

where  the  bar  is  estimate  and  the  bracket  is  expected  values 


p  =  Atmospheric  parameter  vector  (31) 
d  =  Measured  data  vector  (8) 

D  =  Rectangular  31  X  8  matrix.  Determined  by  minimizing  the  expected 
value  of  the  square  of  the  difference  between  the  predicted  vector 
and  the  true  vector.  Stratified  by  season,  zone,  and  view. 

-  Air  Temperature:  1000,  850,  700,  500,  400,  300,  250,  200,  150,  100, 

70.  50.  30.  20,  and  10  mbs 

-  Atmospheric  Thickness:  1000-050.  850-700,  700-500,  500-400,  400-300 

300-250,  250-200,  200-150.  150-100,  100-70,  70-50,  50-30.  30-20, 
20-10  mbs 

-  Tropopaufle  Temperature  and  Pressure 

d)  Atmospheric  Heights:  Stacks  the  thickness  values  on  retrieved  1000  mb 
height  i  rui.'  HIRAS.  1000,  850,  700,  500,  400,  300,  250,  200,  150, 

100,  70,  50,  30.  20,  10  mbs 
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SPECIA1.  SENSOR  MICROWAVE  TEMPERATURE  SOUNDER  (SSM/T-1)  Cent. 


4)  Verification  and  Matrix  Update 


a)  Rea i dual  orrorfl  in  the  derived  atmospheric  profiles.  Evaluate: 
effectiveness  of  inversion  matrix 


-  Bias  and  RMS  statistics  for  latitude  zone  and  pressure  level 

-  Rawinsonde  and  rocketsonde  data  used  as  truth 

-  Rawinsonde  collocates  it  within  3  hours  and  ICO  Nautical  miles  (IJm) 

-  Rocketsonde  used  if  within  9  hours  and  200  Nm  of  rawinsonde 

-  Uses  climatology  above  10  mb  ii  no  rocketsonde  is  found 


b)  Sensor  random  and  bias  errors.  Evaluates  sensor  performance  and 
determine  if  an  update  is  required. 


-  Coincident  rawinsonde/rocketsonde  and  microwave  measurements  used  to 

derive  a  data  covariance  matrix 

-  Compute  error  covariance  matrices  (sensor  random  and  bias  error) 

-  Estimate  a  and  fi  regression  coefficients  from  DM  (measured)  and  DC 

(calculated) 

-  If  statistical  tests  are  met,  update  the  D  inversion  matrices 


Level 

Rpts 

Mean 

RMS. 

Rea 

Mean 

RMS 

Rciii 

Msaa  RMS. 

10  MB 

484 

-1.9 

3.4 

591 

.3 

5.0 

134 

1.4 

3.4 

20  MB 

656 

-1.8 

3.0 

1335 

-1.3 

3.2 

384 

.2 

2.7 

30  MB 

721 

-.6 

2.0 

1619 

-.3 

2.7 

487 

-.1 

2.1 

50  MB 

811 

.5 

1.8 

1801 

-.1 

2.5 

552 

.7 

2.5 

70  MB 

863 

.5 

1.7 

1886 

-.2 

2.6 

596 

1.8 

3.4 

100  MB 

934 

.4 

1.6 

2078 

-.4 

2.0 

671 

-1.8 

3.1 

150  MB 

934 

-.2 

1.9 

2078 

-.6 

2.4 

671 

-3.5 

4.5 

200  MB 

934 

-.7 

2.7 

2078 

-.5 

2.9 

671 

-1.8 

3.3 

250  MB 

934 

.0 

2.4 

2078 

-.0 

3.1 

671 

.2 

2.2 

300  MB 

934 

.0 

2.5 

2078 

-.3 

2.6 

671 

.5 

2.1 

400  MB 

934 

-.1 

2.2 

2078 

.1 

2.3 

671 

.5 

2.1 

500  MB 

934 

-.3 

2.2 

2078 

.1 

2.3 

671 

.7 

2.0 

700  MB 

934 

-.3 

2.3 

2063 

1.5 

2.9 

671 

.9 

2.1 

650  MB 

909 

-.5 

3.3 

2013 

1.8 

3.8 

658 

1.0 

3.8 

1000  MB 

398 

2.2 

4.7 

622 

-.7 

4.2 

41 1 

2.3 

3.8 

F- 10  Error  Stats  10  Apr.  92 
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SPECIAL  SENSOR  MICROWAVE  WATER  VAPOR  PROFILER  (SSM/T-2) 


l )  Padfl i vo ly  col  tecta  microwave  energy 


a ) 

Fivn  Channels: 

:  183 

310 

GHz 

*  / 

1 

wa  t  o  r 

vapor  rofion.'iiioo 

1 1  ii>' 

183 

3  10 

GHz 

<  /  - 

3 

water 

vapor  resonance 

line 

183 

310 

GHz 

♦  /  - 

7 

water 

vapor  resonance 

1  i  ne 

150 

ooo 

GHz 

low  humidity  channel 

91 

658 

GHz 

window  channel 

b)  Requires  Four  SSM/T-1  Channela:  50.50  GHz  oxygen  absorption 

53.20  GHz 
54.35  GHz 
54.90  GHz 

c)  Scan  Type:  Spin  scan,  cross-track  nadir,  synchronized  with  SSM/T-1 

d)  Positions:  28  cross- track,  four  calibration  (hot  and  cold! 

View  angles;  0  to  40.5  in  3  degree  increments 
Total  scan  time  is  8  second 

e)  Resolution:  50  km  at  nadir  through  75  km  at  view  angle  of  40.5 


Man4« 


fwncN«M  (\<«  *} 


Mm  CMqwfU.  A.  *•«*«*  ()/l|/«t| 


53 


SPECIAL  SEKSOR  MICROWAVE  WATER  VAPOR  PROFILER  (SSM/T-2)  Cont. 


2)  Calculates  Sensor  Data  Records  (SPR's):  Kelvin 

a)  Earth  Located:  Accounts  for  earth  rotation  within  1  minute  ephemens 

block,  ellipsoidal  model  replaced  by  exact,  geometry, 
adjustable  reference  height. 

b)  Calibrated:  Total  Power  radiometer;  channel  radiometric  voltages,  cold 

and  hot  temperature  voltages  (8  scan  average),  linear 
voltage  to  temperature  conversion. 

c)  Terrain  and  1000  mb  Heights:  Obtains  values  from  datatase  and  SSM/T- 1 

3)  Calculates  Environmental  Data  Records  (EDR’s) 


a)  Terrain  Height  Corrections:  SSM/T- l  background  eifect  correction, 

five  reference  heights  (0,  1.5,  3.0,  4.5,  6.2  km),  uses  linear 
interpolation 

b)  Inversion  Method  Based  on  Dynamic  Stratification:  uses  multiple  linear 

regression,  four  background  types  (land,  ocean,  coast,  ice)  with 
stratification  based  on  air  mass  characteristics  (25),  quality  flag 
for  wet  profiles  (.15  kg/m**2) .  Water  vapor  weighting  functions  show 
wide  range  of  variability. 


Background  Index  Air  Mass 


Ocean : 

01  - 

05 

Vapor 

Mass 

,  Cold,  High  Altitude  Vapor 

Land  (0000 

m)  : 

06  - 

08 

Cold  . 

High 

Altitude  Vapor,  Not 

Land  (1500 

m)  : 

09  - 

11 

* 

Land  (3000 

m)  : 

12  - 

14 

* 

Land  (6200 

m)  : 

18  - 

20 

• 

Sea  Ice: 

21  - 

22 

Cold  , 

Not 

Coast : 

23  - 

25 

Cold  , 

High 

Altitude  Vapo",  Neither 

c) 

Relative 

Humidity : 

1000 

.  850,  700,  500, 

400, 

and 

300  mb 

levels  in  V. 

d) 

Speci f ic 

Humidity : 

1000 

,  850,  700,  500, 

400, 

and 

300  mb 

levels  in  G/KG 

e) 

Water  Vap 

nr  Mass: 

SFC  - 
400- 

1000.  1000-850.  050-700, 
300,  and  above  300  mbs  in 

700-500, 
KG/M«  *2 

500-400, 
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USERS  OF  DMSP  MICROWAVE  RADIOMETER  DATA 


1)  HI RAS :  High  Resolution  Analysis  System 

-  Provides  initial  conditions  for  GSM  and  RWM 

-  Uses  SSM/T-1  temperature  and  height  data 

-  Evaluation  to  incorporate  SSM/T-2  moisture  data 

2)  RWAM :  Relocatable  Window  Analysis  Model 

-  Planned  to  use  both  SSM/T-1  and  SSM/T-2  profiles 

3)  ASPAM :  Atmospheric  Slant  Path  Analysis  Model 

-  Uses  SSM/T-1  temperature  and  height  data 

-  Will  incorporate  SSM/T-2  moisture  data 

-  Driving  force  behind  SSM/T-2 

4)  RTNEPH :  Real  Time  Nephanaiysifl  Model 

-  Uses  SSM/ I  data;  SDR's  for  earth  skin  temperature 

5)  AGRMET :  Agricultural  Meteorological  Model 

-  Uses  SSM/I  data;  rain  rates  over  land 

-  Interested  in  soil  moisture 

6)  MI ST 1C :  Mission  Sensor  Tactical  Imaging  Computer 

-  SSM/I  imaging  and  fast.  F.PR’s 

-  Technique  development  and  validation 

7)  SDHS :  Satellite  Data  Handling  System 

-  SSM/T-1  and  SSM/I  displays  to  enhance  weather  forecasts 

8)  GEON:  Global  ETAC  OL-A  Network 

-  Satellite  transmission  of  databases  to  Asheville,  NC 

9)  SSM/T-2  CAL/VAL:  SSM/T-2  calibration  and  validation  study 

-  Process  real  time  SSM/T-2  for  data  collection 

-  Process  other  data  sets  required  for  validation 

10)  Tropical  Forecast:  Tropical  Storm  fixing 

-  Uses  SSM/I  data;  rain  rates  and  surface  wind  speed 

11)  Shared  Processing:  Near  real  time  transmission  to  ENOC  and  NOAA  NESDIS 

-  Raw  data  (prepfiles)  SSM/I,  SSM/T-1, and  SSM/T-2 
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SPECIAL  SEMSOR  MICROWAVE  IMAGER  SOUNDER  (SSMIS) 


1 )  One  aenflor  for  temperature,  thickneaa.  moisture,  and  imaging 


a)  24  Channels:  Frequencies  range  from  19  to  183  GHz. 


c) 

cl) 


e ) 


Channel 

Frequency (GHz) 

Polarization 

Footpr i nt (km) 

Comments 

1 

50.3 

H 

37.7 

X 

38.8 

Sound i ng 

LW 

2 

52.8 

H 

37.3 

X 

38.8 

Sounding 

LW 

3 

53.596 

H 

37.7 

X 

38.8 

Sounding 

LW 

4 

54 . 4 

H 

37.7 

X 

38.8 

Sounding 

LW 

5 

55.5 

H 

37.7 

X 

38.8 

Sounding 

L 

6 

57.29 

V 

37.7 

X 

38.8 

Sounding 

L 

7 

59.4 

V 

37.7 

X 

38.8 

Sounding 

L 

8 

150.0 

H 

13.2 

X 

14.7 

Imaging 

37 . 7 

X 

38.8 

Sounding 

W 

9 

183.31  +-.7 

H 

13.2 

X 

14.7 

Imaging 

37.7 

X 

38.8 

Sounding 

W 

10 

183.31  +-.3 

H 

13.2 

X 

14.7 

Imaging 

37.7 

X 

38.8 

Sounding 

W 

1 1 

183.31  +  -.1 

H 

13.2 

X 

14.7 

Imaging 

37.7 

X 

38  8 

Sounding 

W 

12 

19.35 

H 

44.7 

X 

73.6 

Imaging 

13 

19.35 

V 

44.7 

X 

73.6 

Imaging 

14 

22.235 

V 

44.7 

X 

73.6 

Imaging 

15 

37.0 

H 

31.2 

X 

45.0 

Imaging 

16 

37.0 

V 

31.2 

X 

45.0 

Imaging 

17 

91.655 

V 

13.2 

X 

14.7 

Imaging 

18 

91.655 

H 

13.2 

X 

14.7 

Imaging 

37.7 

X 

73.6 

Sounding 

W 

19 

63.283  *-.285 

H  +  V 

75.2 

X 

75.0 

Sounding 

U 

20 

60.792  +-.357 

H  +  V 

75.2 

X 

75 . 0 

Sounding 

U 

21 

60.792  +-.359 

H  +  V 

75.2 

X 

75.0 

Sounding 

U 

22 

60.792  +-.363 

H+V 

75.2 

X 

75.0 

Sounding 

u 

23 

60.792  +-.373 

H  +  V 

75.2 

X 

75.0 

Sounding 

u 

24 

60.792  +-.408 

H  +  V 

75.2 

X 

75.0 

Sounding 

u 

37.2 

X 

38.8 

Sounding 

L 

Scan  Type : 

45  degree  conical  scan  with  concurrent 

beams 

Scan  Rate: 

31.6  RPM,  one 

revolution  in  1 

.9  seconds 

Swath  Width- 

143  2  degrees 

equals  1707  km 
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SPECIAL  SENSOR  MICROWAVE  IMAGER  SOUNDER  (SSMIS)  Coni. 


2)  Summary  of  Changes  from  Existing  Sensors:  SSM/T-1,  SSM/’T'-2,  and  SSM/I 

a)  Earth  Located:  Three  reference  heights,  improved  accuracy 

-  Environmental  Data:  7  km  at  sea  level 

-  Lower  Air  Sounding:  12.5  km  at  H  =  11  km 

-  Upper  Air  Sounding:  12.5  km  at  H  r  60  km 

b)  Lower  Air:  Temperature,  thickness,  water  vapor 

-  Makes  use  of  channels  18  (W)  and  24  (T) 

-  Dynamic  stratification  uses  tropopause  pressure  to 

choose  between  regression  matrices 

-  Atmospheric  thicknesses  derived  from  hypsometric 

equation 

c)  Upper  Air:  Temperature,  thickness  (7,  5.  2.  1.  0.4,  0.2,  0.1,  0.03  mb) 

-  Magnetic  dipole  moment  of  oxygen  complicates  radiative 

transfer  (shifts  weighting  functions).  Solar  storms 

-  Doppler  shift  due  to  earth’s  rotation 

-  Instability  of  center  frequencies  (narrow  band  width) 

d)  Imaging:  Similar  environmental  parameters  and  algorithms 

-  Additional  rainfall  information  from  150  GHz  and  183  GHz 

-  Uses  91.65  GHz  instead  of  85.5  GHz 


Absorption  Coallldant  Ol  Dry  Atr 
(T  =  293  K,  So  Laval) 


Vapor  Sounding 


v  (GHr) 


Absorption  Coelllclant  Ol  Watar  Vapor  In  Air  At  Sea  Laval 
(T  =  293  K,  P  =  7.5  g/n>’) 


Vapor 


Abtorpllon  Coalllclanl  Ol  Small  Walar  Droplalf 
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DMSP  SPECIAL  SENSOR  MICROWAVE/IMAGER  PROCESSING 
AT  THE  FLEET  NUMERICAL  OCEANOGRAPHY  CENTER: 
OPERATIONAL  STATUS,  APPLICATIONS  AND  PLANS 


Marie  C.  Colton  and  C.  James  Cornelius 


Fleet  Numerical  Oceanography  Center 
Satellite  Division 
Monterey,  CA  93943 


The  most  recent  Defense  Meteorological  Satellite  Program  (DMSP)  polar  orbiting 
satellites  have  added  a  microwave  component  to  their  visible/inffared  sensor  suites  by 
including  a  seven-channel,  four  frequency  microwave  radiometer,  known  as  the  Special 
Sensor  Microwave  Imager.  Since  1987  and  the  launch  of  the  first  SSM/I  aboard  satellite 
F8,  the  Fleet  Numerical  Oceanography  Center  (FNOC)  has  had  the  responsibility  of 
processing,  assimilating  and  distributing  the  real-tme  SSM/I  data  to  the  operational  defense 
community.  In  addition,  FNOC  currently  provides  the  data  to  the  Naval  Research 
Laboratory  for  subsequent  archival  at  NOAA/National  Environmental  Satellite  Data 
Information  Service  and  use  by  the  civilian  research  community.  With  the  addition  of  the 
third  SSM/I  aboard  satellite  FI  1  in  November,  1991  and  the  public  release  of  the 
environmental  data  records  from  both  satellites  (summer,  1991),  it  is  now  possible  for 
users  to  demonstrate  the  usefulness  of  global  microwave  data  and  assess  its  impact  on 
operational  meteorology. 

The  purpose  of  this  paper  is  to  provide  a  way  point  to  current  and  potential  users  of 
the  SSM/I  data  via  an  outline  of  the  operational  status  of  SSM/I  data  processing  at  FNOC. 
The  computer  platforms  on  which  the  data  are  received,  processed  and  displayed  are  briefly 
described,  as  well  as  statistics  and  relevant  information  regarding  the  acquisition  of  the  data 
thus  far.  Basic  characteristics  of  the  SSM/I  software,  delivered  by  Hughes  Aircraft 
Company,  including  the  algorithms  by  which  microwave  brightness  temperatures  are 
converted  to  geophysical  parameters  are  listed.  Standard  daily  products  and  applications 
which  use  the  environmental  parameters  in  both  image  and  digital  formats  and  the  means  by 
which  these  products  are  currently  being  communicated  to  regional  sites  are  described  The 
development  of  new  applications  of  SSM/I  data  is  predicated  on  the  amount  of  data  sharing 
that  occurs;  therefore,  near-term  and  future  methods  of  data  communications  are  discussed. 
The  paper  concludes  with  a  discussion  of  planned  improvements  and  future  goals  of  the 
SSM/I  processing  at  FNOC. 
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FNOC  SSM/I  Processing  : 
Status,  Applications,  and  Plans 


Mission  Objectives 
Operations  Focus 
Processing  and  Display  Hardware 
Ground  Processing  Software 
toducts- 


Mission  Objectives 


•  Provide  operational  Navy  users  with  global, 
environmental  products  from  the  DMSP  SSM/I 

•  Provide  SSM/I  brightness  temperatures  and 
environmental  data  records  to  the  Shared 
Processing  Network  (SPN)  partners,  as  set 
forth  in  the  SPN  MOA  of  18  Jun  1984 
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•  EARLY  ORBIT  SUPPORT  (as  required) 

•  INPUT 

•  Complete  acquisition  of  FI  0  &  F1 1  DMSP  SSM  / 1  &  OLS 

•  PROCESS 

•  Benchmark  SSM/I  software  and  algorithms 
•Accurate  geolocation:  onboard  ephemerls  (7km) 

•  SSM/1  Software  Configuration  Management 

•  OUTPUT 

•  Dally  tape  archival  to  NESDIS  via  NRL  (5/day) 

•  Oigltal  data  to  CNODDS  (NOAA,  Monterey) 

•  Ice  concenvra'.lon  product  to  JIC/NPOC 

•Wind  speea  to  Navy  Op.  Global  Atm.  Pred.  Sys.  (NOGAPS) 

•  Integrated  water  vapor  to  moisture  analysis 

•  Imagery  extractions  from  “Quadsphere"  data  base 

•  Shared  Processing  Network  transmissions 
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Operational  Software  Status 


•  Software  updated  to  include  final  changes  recommended  by 
DMSP  SSM/I  Cal/Val  team  in  their  Final  Report  Vol  II,  May 
1991 

•  Changes  Incorporated  into  FNOC  and  AFGWC  software  by 
Hughes  programmers  via  direction  from  System  Program 
Office,  after  formal  recommendation  from  NRL. 

•  Installed  on  3  Dec  9,  effective  with  F10  #5251  and  F11  SSM  / 1 
power-up 

•  Copied  to  SUN  for  Source  Code  Configuration  System 

•  Updates  pending:  SPN  file  selection  and  Small  Area 
transmits  to  regional  centers 
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Data  Communication  Networks 


•  Navy  Environmental  Data  Network  (NEDN) 

•  Four-bit  images  at  9.6Kbaud 

•  Transmitted  to  Regional  Centers  at  Guam, 
Pearl  Harbor,  Norfolk,  Suitland  and  Rota 
(Spain) 

•  Shared  Processing  Network  for  digital  data 
sharing  among  FNOC,  NOAA/NESDIS, 
AFGWC  and  NAVOCEANO  (7  Jan  92) 


Digital  files  to  Guam,  Rota,  NPOC  via  Defense 
Data  Network 


t.-CC  14  flp'r  92 


Processing  Plans 


•  Improve  reliability  of  raw  data  acquisition 
phase 

•  Exploit  move  to  Open  System  architecture 
hardware 

•  Transfer  data  from  Concurrent  to  Sun  via 
LAN  (TCP/IP) 

•  Migrate  SSM  / 1  processing  to  Sun  Galaxy 

•  Store  archival  data  sets  on  optical  discs 


•  Turnover  operational  data  archival 
responsibility  to  NOAA 


RB3@^ 


Product  Development 


Subarea  digital  files  in  standard  format  (eg.,  BUFR) 

Graphic  plots  :  ice  concentration  values 

Contour  plots  :  30  kt  wind  radius 

Multiplexed  images  and  conditional  extractions 

Simplified  SDR  and  EDR  files 

Expanded  data  base  :  land  parms  to  model  input  files 

Rain  "Index"  field 


Moisture  input  to  NOGAPS  (NRL) 


. 

ii/  :  1  *4 ♦  ? ! f * 'V.* i* f * '* 

t  *  ' 

/m. 
:ft»:  /  A 

/  M 

~  •  _ 

SUN  Hosted  Applications 


Mission  Sensor 
Tactic*!  Imaging 
Computer  (Hughes 
Aircraft  Co.) 


Navy 

Environmental 


Operational 
Nowcasting  System 
(developed  at 
NOARL 


Atmospheric  Du.) 


Image  Analysis  Precision  Visual  Co. 

software  (developed 
at  Univ.  of 
Colorado) 
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NODDS  vugraph 


ISSUES 


•  SSM  / 1  Configuration  Management 

•  Mechanism  for  evolution  of  processing 
software  and  algorithms 

•  NESDIS/NGDC  Archival 

•  SPN  automation  enhancements  (polling  v. 
windowing) 
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Summary 


•  FNOC  processes  all  SSM  / 1  data  from  two  DMSP  Satellites 
(F10  and  F11) 

•  The  operational  software  uses  the  most  current  algorithms  as 
specified  in  the  DMSP  Cal  /  Val  Report,  Vol  II  (1991) 

•  The  benchmark  software  is  under  version  control  at  FNOC 

•  Standard  products  include  wind  speed  assimiliation  by 
NOGAPS 

•  SPN  transmits  occur  24  hr/day 

•  Data  archival  to  tape  continues 

•  Current  emphasis  at  FNOC  is  on  improvement  of  graphic 
displays  and  SSM  / 1  data  communications,  including  SPN 


AN  OVERVIEW  OF  THE  EXPERIMENTAL  SSM/I  ORBIT-BY-ORBIT 
PRODUCTS  SYSTEM  (SSMIPROD)  AT  NOAA/NESDIS 

Joseph  V.  Fiore,  Jr. 

S  M  Systems  and  Research  Corporation 
8401  Corporate  Drive  (Suite  510) 

Landover,  MI  20745-2224 


A  prototype  products  system  for  deriving  experimental  SSM/I  orbit-by-orbit 
products  (SSMIPROD)  has  been  developed  at  NOAA/NESDIS.  SSMIPROD  was 
developed  to  obtain  a  better  understanding  of  the  accuracy  of  the  SSM/I  geophysical 
algorithms,  to  enhance  the  reliability  of  pre-launch  algorithms  being  developed  for  the 
Advanced  Microwave  Sounding  Unit  (AMSU),  and  to  gain  experience  in  delivering 
microwave  products  to  NOAA  users  such  as  National  Weather  Service/National 
Meteorological  Center/Development  Division  (NWS/NMC/DD). 

SSMIPROD  has  been  running  in  a  routine  "experimentally  operational" 
environment  since  July  26,  1990.  The  SSMIPROD  output  file  contains  two  days  of  SSM/I 
orbit-by-orbit  products  stored  in  a  rotating  direct  access  file  on  the  NOAA  HDS  computer. 
Currently,  there  are  1 1  products  (six  ocean,  three  sea  ice,  and  two  land)  being  generated  by 
SSMIPROD.  The  ocean  products  include  the  output  of  two  precipitation  algorithms,  two 
near  surface  wind  speed  algorithms,  and  two  total  precipitable  water  (TPW)  algorithms. 
The  sea  ice  products  include  the  output  of  total  ice  concentration  from  three  sea  ice 
algorithms.  The  land  products  include  the  output  from  two  precipitation  algorithms. 

Extensive  analysis  and  evaluation  of  these  products  has  been  performed  at 
NOAA/NESDIS.  A  plan  for  the  validation  of  the  products  was  developed  to  test  and 
evaluate  the  quality  of  the  orbit-by-orbit  products  more  completely.  This  plan  involves 
using  many  independent  data  sets  as  "ground  truth"  for  validating  the  products.  In  this 
paper  the  validation  of  the  TPW  algorithms  using  matched  radiosonde  data  will  be 
discussed.  The  quality  of  the  precipitation  algorithms  will  also  be  discussed. 
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OUTLINE  OF  TALK 


1.  Overview/Mistory  of  SSM/I  Processing  at  NOAA/NESDIS 

2.  Data  Flow:  SSM/I  Products 

3.  Current  SSM/I  Orbit-by-Orbit  (SSMIPROD)  Products 

4.  SSMIPROD  Products/Images 

5.  Summary/Plans/Related  Talks 


History  of  SSM/I  processing  at  NOAA/NESDIS 

NESDIS  -  has  a  shared  processing  agreement  Air  Force,  Navy 
NESDIS  COE  for  Polar/DMSP  Soundings 


NESDIS  Office  of  Research  and  Applications  (ORA)  receives 
raw  SSM/I  data  near  real  time  from  AFGWC  (Orbit-by-Orbit)  for 

Sounding  Improvement 


NESDIS  creates  an  SSM/I-lb  data  set  with  earth  located 
radiances  (since  1989) 


SSMIPROD  (Orbit-by-Orbit)  SSM/I  Products  -  Running 

"Experimentally"  since  July  26,  1990 


SSMIPROD  developed  to  evaluate  SSM/I  Cal/Val  algorithm 
products  prior  to  routine  availability  from  FNOC 


MASTERMAPS  -  mapped  SSM/I  products  (512  X  512  (l/8th  mesh)), 
(1024x1024  (l/16th  mesh))  -  polar  stereographic 
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CURRENT  "EXPERIMENTAL"  PRODUCTS  (F-10  and  F-ll 


OCEAN 

WIND  SPEED  (CAL/VAL) 

TEST  GOODBERLET  WIND  SPEED  (CAL/VAL) 

TOTAL  PRECIPITABLE  WATER  (CAL/VAL) 

TOTAL  PRECIPITABLE  WATER  (PETTY) 
RAIN  RATE  (GRODY) 

RAIN  RATE  (CAL/VAL) 


LAND 

RAIN  RATE  (CAL/VAL) 
RAIN  RATE  (GRODY) 
Not  Used 
Not  Used 
Not  Used 
Not  Used 


ICE 

CAL/VAL  ICE  CONCENTRATION 
CAL/VAL  ICE  TYPE 

GRODY  ICE  CONCENTRATION 
NASATEAM  TOTAL  ICE  CONCENTRATION 
NASATEAM  FIRST  YEAR  ICE  CONCENTRATION 
NASATEAM  MULTI  YEAR  ICE  CONCNETRATION 
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Ocean:  Wind  Speed  (in/sec): 


Cal/Val:  Linear  I9V.22V.37V.37I  I  'll)  (quality  flags) 

Tb’s 

Quality  Flags 
range:  l-25m/sec 

'lest  Cal/Val:  Similar  to  Original  Cal/Val 

Weather  Filter  (developed  for  Typhoon  Wind  Speeds) 
Quality  Flags 


Ocean:  Total  Precipitable  Water  TPW  (mm): 

Cal/Val:  Non-Linear  (exponential)  19V,  22V,  37V 

Tb’s 

range:  10mm-60mm 

TPW  =  232.89-0.148*Tbl9V-1.829*Tb22V-0.369*Tb37V-0.006193*Tb22V**2 


Genera  1  : 


SSMI PROD  Ice  Products 


All  %  <  0  set  to  0;  all  %  >  100  set  to  100 

Flags:  0  -  25%-100% 

1  -  0%-25% 

7-  <  0%  or  >  100% 

Climatological  Ice  Limit  (current  vs.  latest  JIC  (G.  Wohl)) 
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Ice  Products 


1.  Navy  Cal/Val  Ice  Concentration: 

Linear  (19V,  19H,  37V,  37H) 

Winter  and  Summer  Coefficients 

%  Ice  Concentration  (25%-100%),  Ice  Type  (Fyr,  Myr) 

2.  Grody  Ice  Concentration 

Emissivity  based  Linear  (19V,  22V) 

Variable  Emissivity  Threshold  =  0.7  (current) 

3.  NASATEAM  Ice  Concentration 

Linear  (19V,  19V,  37V) 

Tie  Points  (winter,  summer)  -  Open  water,  ice 
covered,  Gradient  Ratio  (weather  filter) 

Total,  Firstyear,  Multiyear  (%) 

Grody  Decision  Tree  Algorithm 
Non-Linear  (19V,  19H,  22V,  37V,  85V) 

SIS 5  and  SI37  (If  SI37  +  3  >  SI85  use  SI37) 

85GHz  Precipitation  (rain  rate  mm/hr)  Ocean  and  Land 
Snow  Cover  Index  (Land) 

Desert  Sand  (screened  out) 

Vast  improvement  over  37GHz  algorithms  (especially  over  land) 
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SUMMARY/PLANS: 


Two  "Experimental"  SSMIPROD  data  sets:  F-10  and  Fll 

FNOC  data  (EDR's)  Primary  Source  of  SSM/I  data  1992  (reformat) 

AFGWC  SSM/I  data  will  be  used  as  a  backup 

Plans  to  make  SSMI PROD/MASTERMAPS  Operational: 

Merge  code,  test  products,  Transfer  to  CEMSCS 

VALIDATION:  Q.M.  Winds  already  developed  (continue),  TPW, 
Precip,  Snow  Cover,  Sea  Ice 


Related  Talks  (Users) : 

"Real  Time  Quality  Control  of  SSM/I  Wind  Speed  Data"  M.  Waters 
W.S.  Richardson,  W.H.  Gemmill,  C.M.  Caruso  (NOAA)  (SSMIPROD) 

"Uses  of  SSM/I  Wind  Speed  data  in  Operational  Numerical 
Weather  Prediction  System  at  NMC"  T.W.  Yu,  W.H.  Gemmill,  J. 
Wollen  (NMC/DD)  (SSMIPROD) 

"Impact  of  SSM/I  based  Snow/Ice  Analyses  in  NMC's  Eta  Model  - 
K.  Mitchell  (NMC/DD)  (MASTERMAP) 

Related  Talks  (NOAA/NESDIS) : 

Sea  Ice  and  Wind  Products  Production  and  Evaluation  -  W. 
Pichel  (SSMIPROD) 

Grody/Ferraro:  Use  of  SSM/I  for  Precipitation  at  NOAA/NESDIS 
SSM/I  MASTERMAP  Products:  D.  Donahue  C.  Boettcher  (SMSRC) 
Operational  Plans  for  SSM/I  PRODUCTS:  P.  Taylor 
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SSM/I  MAPPED  EXPERIMENTAL  ENVIRONMENTAL  IMAGE 

PRODUCTS 

David  R.  Donahue  and  Corey  M.  Boettcher 

S  M  Systems  and  Research  Corporation 
8401  Corporate  Drive,  Suite  510 
Landover,  MD  20785 

Pamela  M.  Taylor  and  Norman  C.  Grody 

NOAA/NESD1S 
Washington,  DC  20233 


Globally  mapped  experimental  image  data  of  snow  cover,  sea  ice,  precipitation 
scattering,  rain  rate,  and  total  precipitable  water,  derived  from  the  Special  Sensor 
Microwave/Imager  (SSM/I),  are  now  updated  on  a  daily  basis  at  NESDIS.  NESDIS 
receives  the  SSM/I  data  in  near  real  time  via  the  Shared  Processing  Network  between 
NESDIS,  the  Navy,  and  the  Air  Force.  This  data  is  available  to  NOAA  users  on  the  HDS 
computer  facility  in  the  form  of  IB  data  sets,  which  include  earth-located  raw  counts  and 
calibration  parameters,  produced  as  soon  as  the  orbital  data  are  received  at  NESDIS.  Once 
transferred  to  the  HDS  system,  each  SSM/I  IB  data  set  is  calibrated  and  an  Intermediate 
Data  Base  (IDB)  of  antenna  temperatures  is  generated.  The  data  in  the  IDB  is  then 
processed  by  the  DMSP  SSM/I  Master  Maps  software,  which  provides  gridpoint  analyzed 
fields  of  SSM/I  channel  measurements  and  derived  products  for  access  by  users. 
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SSM/I  MAPPED 
EXPERIMENTAL 
ENVIRONMENTAL 
IMAGE  PRODUCTS 


Corey  M.  Boettcher  and  David  R.  Donahue 

S  M  Systems  and  Research  Corporation 
Landover,  Maryland 


Pamela  M.  Taylor 
Norman  C.  Grody 

NOAA  /  NESDIS 
Washington,  D.C. 


OVERVIEW 


•  NOAA/NESDIS  receives  SSM/I  data  in  near  real-time,  currently  over  the 
Shared  Processing  Network  from  Air  Force  Global  Weather  Central. 


•  The  Sounding  Implementation  Branch  (SIB)  has  developed  the 
Experimental  DMSP  SSM/I  Master  Maps  Image  Product  System  as  part 
of  an  elTort  to  use  SSM/I  data  to  improve  the  DMSP  SSM/T  temperature 
soundings.  This  work  also  enhances  the  ability  of  NESDIS  to  process  and 
distribute  the  Fleet  Numerical  Oceanography  Center  SSM/I  products. 


0  As  each  SSM/I  data  set  is  received  and  processed,  SIB  has  automated 
software  which  computes  and  maps  environmental  products.  These 
products,  along  with  the  antenna  temperatures  and  lime-identifying 
parameters,  are  written  as  image  pixel  values  onto  individual  gridded 
polar-stercographic  fields. 


•  The  image  files  containing  these  gridded  fields  are  accessible  on  the 
NESDIS  mainframe  to  all  interested  users.  Testing  of  these  experimental 
products  has  shown  their  usefulness  in  improving  SSM/T  soundings  and 
in  many  other  weather  and  climate  studies. 


•  As  a  result  of  the  interest  in  these  experimental  SSM/I  products,  SIB  has 
received  formal  requests  from  the  National  Meterological  Center,  the 
National  Hurricane  Center  and  the  Navy/NOAA  Joint  Ice  Center  for 
operational  access  to  these  products.  Work  is  underway  to  make  these 
experimental  SSM/I  Master  Maps  products  operational. 
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Diagram  of 
SSM/I  Data 
As  Recorded  By 
DMSP  Satellite 


Resulting  DMSP 
SSM/I  Master  Maps 
Polar-Stereographic 
Image 


l)MSI*  SSM/I  MASTER  MAPS  I  MAC  1C  FILES 


One  Image  File  Per  Hemisphere 
Per  Resolution 
Per  Satellite 


MAPPEQ  RESOLUTIONS 


1  /8th  Mesh 

512  X  512  Pixel 

47.6  km  @  ±60°  latitude 

Satellites:  F-8  (no  85  GHz) 

F-10 
F-ll 


1/ J 6th  Mesh 

1024  A  1024  Pixel 

23.8  km  @  ±60°  latitude 

Satellites:  F-10 

F-ll 
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MAITKI)  SSM/I  I'AKAMKTKKS  IN  KACII  IMA<;K  KILL 


ALGORITHM 

CURRKNT;  Channel  Antenna  Temperatures 

Latitude  &  Longitude  of  source  data 
for  each  pixel 


Time-Identifying  Information 
fi.e.  Julian  Day  &  GMT) 


Solar  Zenith  Angle 

Snow  Cover  Scattering  Index 

N.  Grody 

Sea  Ice  Concentration 

Navy  Cal/Val 

Precipitation  Scattering  Index; 

Land  &  Ocean 

N.  Grody 

Precipitation  Rain  Rate; 

Ocean-Only 

N.  Grody 

Total  Prccipitable  Water 

Navy  Cal/Val 

KLOlLRLj  Sea  Ice  Concentrations:  NASA  team 

Total,  Multi-year  &  First-year 

Precipitation  Rain  Rate;  N.  Grody 

Land  &  Ocean 

Surface  Wind  Speeds  Navy  Cal/Val 
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SUMMARY 


The  development  of  the  DMSP  SSM/I  Master  Maps  products  has 
proven  to  he  excellent  preparation  for  the  NESDIS  processing  and 
distribution  of  FNOC  SSM/I  products. 


The  experimental  SSM/I  Master  Maps  products  provide  a  unique 
troubleshooting  ability  for  both  science  and  software  development.  As 
a  result,  NESDIS  could  potentially  provide  science  support  to  FNOC 
products. 


•  Progress  continues  on  using  the  mapped  SSM/I  products  to  improve 
SSM/T  soundings. 

•  Testing  will  soon  begin  using  the  SSM/I  Master  Maps  sea 
ice  product  as  a  dynamic  ice  field  within  the  DMSP  SSM/T 
sounding  system,  thus  replacing  a  static  ice  field. 

•  Work  is  ongoing  to  use  the  mapped  SSM/I  channel  antenna 
temperatures  to  aid  in  applying  surface  corrections  to  the 
SSM/T  channels. 


For  details  on  how  to  access  these  experimental  SSM/I  Master  Maps  products 
or  for  further  information,  please  feel  free  to  contact: 

Corey  M.  Boettcher 
(301)  763-4384 

Pamela  M.  Taylor 

(301)  763-4380  /  FAX:  (301)  420-0932 

NOAA/NESDIS 
Satellite  Applications  Lab 
Sounding  Implementation  Branch 
E/RA23 

Federal  Building  4 
Washington,  D.C.  20233 
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SSM/I  TOTAL  PRECIPITABLE  WATER  VAPOR  ALGORITHMS: 
A  REPRISE  AND  UPDATE 

J.C.  Alishouse  and  R.R.  Ferraro 

NOAA/NESDIS 
Satellite  Research  Laboratory 
Washington,  DC  20233 


This  presentation  will  provide  a  brief  review  of  the  development  of  the  present 
operational  total  precipitable  water  (TPW)  algorithm  and  the  compilation  of  the  data  base 
used  to  develop  and  validate  it  Other  algorithms  for  the  SSM/I  TPW  have  been  developed 
and  published.  Brief  comparisons  between  the  various  algorithms  and  the  validation  data 
base  will  be  made.  A  continuing  and  real  time  validation  and  quality  control  effort  for 
water  vapor  algorithms  is  being  implemented  as  part  of  the  experimental  products 
development.  Factors  of  relevance  to  future  algorithm  and  instrument  development  will  be 
presented. 
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RETRIEVED  IKC/M.«2J 


SSM/I  TOTAL  PRECIPITABLE  WATER  ALGORITHMS: 
A  REPRISE  AND  UPDATE 

JOHN  C.  ALISHOUSE 

AND 

RALPH  R.  FERRARO 
NOAA/NESDIS 

SATELLITE  RESEARCH  LABORATORY 
PRESENTATION  AT: 

DMSP  RETRIEVAL  PRODUCTS  WORKSHOP 
APRIL  14,  1992 


SSM/I  TPW  vs.  RAOB 


NESD1S  Non-Line«r 
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INDEPENDENT  DATA  SET 


0  30  40  SO 

OBSDCVEZXKO/MmCS 


SSM/I  TB(22V)  VS.  RAOB  TPW 


587  DATA  OBSERVATIONS 


9 


Log(290-TB22VI 


SET 

N 

MEAN 

(KG/M2) 

RMSD 

(KG/M2) 

PCTD 

BIAS 

(KG/M1) 

!  T  SQUARED  ALGC 

DEP 

586 

26.3 

3.0 

11.4 

0 

/A ID 

171 

38.9 

- - -  ... 

_ 

10.8 

-0.7 

1  T™  ONLY  | 

|  DEP 

587 

26.3 

4.3 

16.3 

0 

|  LOG  (290  -  T„22V)  | 

DEP 

587 

26.3 

3.7 

14.1 

Vo 

|  4  LOG  VARIABLES  \ 

|  DEP 

586 

26.3 

3.0 

11.4 

0 

|  IND 

171 

38.9 

4.3 

11.2 

-0.4 

TPW=A0+A1  *  T™  V+A2  *  Tf  V*  A3  *  r”  V*  A,  *(T2B2V)2 


TPW=B0+Y,  b i* LOG (29 0-T/) 

ISSUES  FOR  FURTHER  STUDY 

o  FREQUENCY  OF  WATER  VAPOR  CHANNEL 

REDUCED  SENSITIVITY  AT  LARGE  TPW 
AMOUNTS 

CHANGE?  IN  CONE  ANGLE  FOR  SSM/IS-2 

PRESSURE  BROADENING  EFFECTS  AT  LINE 
CENTER 

SMMR  AND  AMSU-A  WORK  OFF  LINE  CENTER 
S/N  CONSIDERATIONS 

o  VALIDATION  AND  QC  OF  ALGO  RESULTS 

DOES  F8  ALGO  WORK  FOR  F10  AND  F11? 

QC  FOR  RESULTS 

SUMMARY 

o  NEED  WATER  VAPOR  CHANNEL 

o  RETRIEVALS  OF  LARGE  WV  AMOUNTS  (SATURA¬ 
TION) 

o  VALIDATION  AND  QC  FOR  NEW  INSTRUMENTS 


81 


THE  USE  OF  THE  DMSP  SSM/I  FOR  THE  GENERATION  OF 
PRECIPITATION  PRODUCTS  AT  NOAA/NESDIS: 

PART  I:  A  STATUS  REPORT 

PART  II:  SCIENTIFIC  RESULTS 

Norman  C.  Grody  and  Ralph  R.  Ferraro 
NOAA/NESDIS 

Satellite  Research  Laboratory/Land  Sciences  Branch 
Washington,  DC  20233 


The  demand  for  precipitation  measurements  from  satellites  ranges  from 
instantaneous  rain  rates  in  support  of  operational  weather  and  hydrological  forecasts  to 
weekly  or  longer  rain  accumulations  for  agricultural  and  climate  monitoring  needs.  At 
NOAA,  the  demand  for  these  measurements  comes  from  several  divisions  of  the  National 
Weather  Service  and  the  National  Environmental  Satellite,  Data,  and  Information  Service. 
To  support  this  demand,  there  are  presently  several  activities  ongoing  at  the  Satellite 
Research  Laboratory  using  SSM/I  measurements  to  generate  experimental  precipitation 
products. 

The  activities  can  be  categorized  as  the  following: 

(1)  SSM/I  rain  retrieval  algorithm  development  and  evaluation 

(2)  Algorithm  calibration  and  verification 

(3)  Near  real  time  rain  product  generation 

(4)  Applications  for  synoptic  scale  processes 

(5)  Applications  for  climate  scale  processes 

(6)  Applications  for  numerical  weather  prediction 

It  is  the  purpose  of  the  first  presentation  to  discuss  the  objectives,  current  status, 
and  future  plans  of  these  activities.  The  second  presentation  will  discuss  some  of  the  key 
scientific  issues  and  current  results  of  these  projects. 
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The  Use  of  the  DMSP  SSM/I  for  the  Generation 
of  Precipitation  Products  at  NOAA/NESDIS: 
Part  I  -  A  Status  Report 


Ralph  Ferraro 
Norman  Grody 


Summary  of  SSM/I  Rain  Activities 

1.  "Experimental  Operational"  Activities 

o  Orbit  by  Orbit  Products 
o  Mapped  Products 
o  VDUC  Products 
o  Level  1b  Case  Studies 

2.  Research  Activities  (in  support  of  Climate  and  Global.  Change) 

o  Algorithm  Development/Improvements 
o  Algorithm  Calibration 
o  Algorithm  Intercomparisons 
o  Climate  Products 
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VAS  Data  Utilization  Center  (VDUC)  Products 


1.  Composite  product  files  (F-ll)  generated  every  6  hours 

o  Whenever  the  computers  co-operate! 

o  Ocean  -  Winds,  Rain  Rate  (37  GHz),  Sea-ice 

o  Land  -  Rain  Rate  (37  GHz),  Snow  Cover,  Surface  "roughness" 

0  24  HOUR  ROTATING  FILE 

2.  Used  by: 

o  NESDIS/SAB  -  support  operational  rain  fall  forecasts 
o  ORA  -  in-house  R&D 

O  NHC  -  WIND  FIELDS  AND  STORM  CENTERS 

3.  Future  upgrades: 

o  Dual-satellite  system 
o  85  GHz  algorithms 


NESDIS  Levf.l  1b  Data  Sets 


1.  Used  for  in-house,  near-real  time  algorithm  evaluations 

o  "Field  Experiments" 
o  "Independent"  algorithm  performance 
o  SAB  interactions 

o  Extreme  events  investigations 

2.  Sensor  status/intercomparisons  (limited) 


w  /till,  iv/l  KAtJAK  1401  19AU91 


Algorithm  Development/Improvements 


1.  Grody  Global  85  GHz  Scattering  Index  (SI) 

2.  Scattering  vs.  Emission  over  ocean 

3.  Frequency  vs.  Polarization  (AMSU) 

4.  IFOV  VARIABILITY  IMPACT  STUDIES 

5.  SI  VARIABILITY  WITH  SENSOR  AND  GEOGRAPHY 


Algorithm  Calibration 

1.  Co- INCIDENT  SSM/I  AND  RADAR  DATA  SETS 

o  RADAP-II  (NWS) 
o  AMeDAS  (JMA) 
o  CTRMM/Darwin  (NASA)] 
o  [FRONTIERS  (UK) 3 

2.  Pixel  by  Pixel  vs.  Binning 
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CAL/VAL  IMM/HRI  grody  imm/hri 


SSM/I  VS  AMeDAS  Radar 


(iRODV 


0  2  4  6  8  10  12  U  16 


AMEDAS  RADAR  IMM/HR) 

SSM/I  VS  AMeDAS  Radar 


CAL/VAL 


0  2  4  e  S  10  12  H  If. 


AMEDAS  RADAR  IMM/HRI 
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AVERAGE  PREDICTORS  VS.  RR  (mm/hr) 


AMeDAS  -  OCEAN  ONI  V 


Algorithm  Intercomparisons  for  Climate 


1.  CREATF  3  YEAR,  MONTHLY,  CLIMATE  SCALE  (100  KM)  RAIN  AVERAGES 

l.  Task  supported  by  Phil  Arkin  for  Remote  Sensing  Systems  (Wentz) 

3.  Algorithms  to  be  used: 
o  Grody  85  GHz  SI 
o  Barrett/Kidd  (Univ.  Bristol) 
o  Petty/Olson  (Purdue/Wisconsin) 
o  Adler  (NASA)  Climate  Products 

1.  Generate  from  retrospective  data;  to  be  placed  on  MO  disks 

2.  Weekly,  monthly,  etc.  for  >5  years  F-8  SSM/I 

3.  Prototypes  from  Grody  SI  now  being  generated  from  tape 
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PRODUCTION  AND  EVALUATION  OF  EXPERIMENTAL  SSM/I 
ICE  AND  WIND  PRODUCTS  AT  NOAA/NESDIS 

William  Pichel 


NOAA/NESDIS  E/RA 13 
Satellite  Research  Laboratory 
NOAA  Science  Center  Room  102 
Washington,  DC  20233 


Ice  and  wind  products  are  among  the  suite  of  products  being  routinely  produced  on 
an  experimental  basis  from  DMSP  SSM/I  data  at  NOAA/NESDIS.  Each  product  is  being 
produced  from  two  or  more  algorithms  to  evaluate  available  algorithms  as  to  their  utility  to 
operational  ocean  and  atmospheric  analysis  and  forecasting  activities  of  the  Navy/NOAA 
Joint  Ice  Center  (JIC),  the  National  Meteorological  Center,  and  the  NOAA  Ocean  Products 
Center. 

Maps  of  ice  concentration  produced  by  three  algorithms  (Navy  Calibration/Vali¬ 
dation,  NASA  Sea  Ice  Working  Group,  and  NESDIS/Norman  Grody)  are  compared  with 
operational  JIC  ice  analyses  produced  using  visible  and  infrared  satellite  imagery.  Wind 
speed  observations  generated  with  two  algorithms  (Navy  Calibration/V alidation,  and  an 
improvement  on  that  algorithm  obtained  from  Mark  Goodberlet)  are  evaluated  using  closely 
matched  NOAA  moored  buoy  wind  reports. 
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PRODUCTION  AND  EVALUATION  OF 


EXPERIMENTAL  SSM/I  ICE  AND  WIND  PRODUCTS 
AT  NOAA/NESDIS 

o  SSM/I  ICE  AND  WIND  PROCESSING  SYSTEM 
o  WIND  PRODUCTS 
o  ICE  PRODUCTS 


DATA  FLOW:  SSM/I  EXPERIMENTAL  PRODUCTS 


SSM/I 

RAW 

DATA-1 

FROM 

AFGWC 


TARGET 

MATCH 

FILE 


8U0Y 

REPORTS 


DAILY 
CHANNEL/ 
PRODUCT  ' 
MAPS 


ORBITAL 

PRODUCTS 


USERS 

NMC 

CAC 

JIC 

OPC 


SAT/BUOY 

MATCH 

FILE 


GENERATE 
WIND 
[I  STATISTICS 


MONTHLY 

STATISTICS 


SSM/I  MARINE  WIND  SPEED  ALGORITHM 
(SSM/I  CAL/VAL  TEAM) 


=  '  ^7.90  +  1.0969  (Teigv)  -  0.4555  (TB22v) 


-  1  .7600  (Tb37V)  +  0.7860  (TB37H) 


no  1 7ow  i60M  now  now  now  now  now  ioow  oow  tow  row  60w 


NOAA  BUOYS  USED  TO  MONITOR  SSM/I  WIND  SPEED  ACCURACY 

SSM/I  CAL/VAL  WIND  SPEED  VS.  BUOY 


DECEMBER  1991 


BUOY  WIND  SPEED  at  19.5  m  (m/sec) 
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SSM/I 


BUOY  WIND  MATCH  STATISTICS 


# 

CAL/VAL 

TEST 

ALG. 

TIME 

PERIOD 

MATCHES 

R 

RMSD 

R 

RMSD 

AUG 

1991 

258 

0 . 68 

2.04 

- 

- 

NOV 

13-30,  1991 

168 

0 . 65 

2.34 

0.68 

2.37 

DEC 

1991 

189 

0.71 

2.64 

0.73 

2.70 

JAN 

1992 

104 

0.79 

2.33 

0.82 

2.30 

FEB 

1992 

17  4 

0.87 

2 . 15 

0.87 

2.34 

MAR 

1992* 

284 

0.85 

2.17 

0.85 

2.24 

March  1- 

5  (F-8) 

March 

6-31 

(F-10) 

R  =  Correlation  Coefficient 
RMSD  =  Root  Mean  Square  SSM/I  -  Buoy  Difference 


NAVY  CAL/UAL 
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SSM/I  ICE  AND  WIND  PRODUCTS 
FUTURE  PLANS 


SYSTEM 


o  TRANSITION  TO  OPERATIONS 

o  SINGLE  SYSTEM  FOR  ORBIT-BY-ORBIT  AND  MASTER  MAPS 
o  INCORPORATE  FNOC  PRODUCTS 

WINDS 


o  CONTINUE  ALGORITHM  ASSESSMENT  (RAIN  FLAGS) 
o  WIND  MATCH  TWO  SATELLITES 
O  USE  BUOY  AND  WIND  QC  FROM  OPC 
o  EVALUATE  USE  OF  BUOY  CONTINUOUS  WINDS 

ICE 


o  CONTINUE  ALGORITHM  ASSESSMENT 
o  IMPLEMENT  LATEST  NASA  ALGORITHM  UPGRADES 
o  FORMAT  ICE  PRODUCTS  FOR  JIC  DIFAS 
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REAL  TIME  QUALITY  CONTROL  OF  SSM/I  WIND  SPEED  DATA 
Marshall  Waters,  W.S.  Richardson,  W.H.  Gemmill  and  C.M.  Caruso 


NOAA  Ocean  Products  Center 
5200  Auth  Road 
Camp  Springs,  MD  20746 


There  is  a  critical  need  to  increase  the  availability  of  accurate  marine  observations/ 
measurements.  A  new  set  of  satellite  derived  (SSM/I)  ocean  surface  wind  speed  data  is 
now  available  in  real  time  through  the  NOAA  central  computer  facility.  These 
measurements  will  allow  NOAA  to  improve  the  accuracy  of  operational  atmospheric  and 
oceanographic  forecast  models  and  existing  marine  warning  and  forecast  systems,  and  to 
study  and  quantify  the  long-term  implications  of  climate  and  global  change  programs. 

More  accurate  marine  measurements  are  being  made  available  by  QUIPS  (Quality 
Improvement  Processing  System).  This  system,  which  is  operational  at  the  NOAA  Ocean 
Products  Center,  is  described  and  discussed,  as  it  relates  to  the  quality  control  (QC)  of 
SSM/I  wind  speeds.  Discussion  includes  pre-processing  and  comparisons  with  an  NMC 
first  guess  field,  interactive  QC  on  QUIPS  and  QC  feedback  for  improved  operational 
forecasts  and  SSM/I  wind  speed  algorithm  assessment/improvement. 
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U15AL  TIME  QUALITY  CONTROL 


or 

SSM/I  DERIVED  WIND  SPEED  DATA 

o  OVERVIEW 

Data  Flow 

-  QC  of  In-Situ  Data 

o  QC  OF  SSM/I  WIND  DATA 

-  "Super  Obu" 

Comparison  to  a  First  Guess 

o  FUTURE  PLANS 

Interactive  QC 

-  Data  Monitoring  and  QC  Reports 
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SURFACE  OBSERVATIONS/MEASUREMENTS  QC'd 

Four  Synoptic  Periods  Each  Day 

Sea  level  pressure,  wind  speed,  and 
and  air  and  sea  surface  temperatures 

600  -  Ships 
350  -  Buoys 
90  -  CMANs 

SUBSURFACE  MEASUREMENTS  QC'd 

Once  Each  Day 

100  -  Temperature  Profiles 
10  -  Salinity  Profifles 


OPC  REAL-TIME  QUALITY  CONTROL 
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POP  BOX  SETS  AND  CLEARS  ALL  FLAGS] 


REAL  TIME  QUALITY  CONTROL 


OF 

SSM/I  DERIVED  WIND  SPEED  DATA 


o  OVERVIEW 

Data  Flow 

-  QC  of  In-Situ  Data 


o  QC  OF  SSM/I  WIND  DATA 

-  "Super  Obs" 

-  Comparison  to  a  First  Guess 

o  FUTURE  PLANS 

Interactive  QC 

-  Data  Monitoring  and  QC  Reports 


SSM/I  "SU^ER  OBS” 


o  ALL  "0"  level  (no  precip  contamination,  no  winds 

GE  25  m/s)  measurements  are  placed  in  a  1  deg  by 
1  deg  lat/long  bin 

-  Bins  are  determined  by  truncating  latlong 
to  whole  degrees 

The  Super  Ob  is  the  linear  average  of  all 
measurements  in  the  1  deg  by  1  deg  bin 

The  location  of  the  Super  Ob  is  the  linear 
average  of  lat/long's  in  the  bin 

The  time  of  the  Super  Ob  is  the  linear  average 
of  the  times  of  the  Obs  in  the  bin 


o  5,000  -  6,000  "Super  Obs"  a  Synoptic  Period 


o  Data,  within  a  3  hr  window,  are  collected  3  hrs 
into  Synoptic  Period 
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COMPARISONS  TO  A  FIRST  GUESS 


o  First  Guess 

06  Hr  Forecast  of  wind  Speed  45  a  above  the  Sfc 

Global  Data  Assimilation  (GDAS)  run  at  06Z  and  182 
—  Verifies  at  12Z  and  00Z 

Aviation  (AVN)  run  at  00Z  and  12Z 

—  Verifies  at  06Z  and  18Z 

o  Approximately  1.5  %  (75  Super  Obs)  of  the  nearly 
5,000  Super  Obs  per  Synoptic  Period,  DEPART 
by  +/-  15  kts  or  more  from  the  First  Guess 


Comparison  of  SSM/I  Derived  Winds 
To  NvfC  First  Guess 
November  5—19.  1991 


ABS  (SSMTI  -  First  Guess)  >=  15  Kta 
(Approximately  1.5%  of  data) 
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REAL  TIME  QUALITY  CONTROL 


OF 

GGH/L  DERIVED  WIND  SPEED  DATA 


o  OVERVIEW 

Data  Flow 

QC  oC  In-situ  Data 

o  QC  OF  SSM/I  WIND  DATA 
“Super  ObB" 

-  Comparison  to  a  Firet  Cuecs 


o  FUTURE  PLANO 

Interactive  QC 

-  Data  Monitoriny  and  QC  Reports 


OPC  REAL-TIME  QUALITY  CONTROL 


OBSERVATIONS/ 

MEASUREMENTS 


FIRST  GUESSES/ 
CLIMATOLOGIES 


OMNET 


□ 


jojaaitl 


DATA 

QC 

MONITORING/ 

REPORTS 

- 

Mot 

- 

Buoy 

- 

SSM/I 

- 

Batliy 

- 

Sa  1  lnlty 

QUIPS/ 

SMARTQC 


PLATFORM  STATS 
HISTORIES 

NEIGHBORING  PLATFORMS] 
QC  RULES 


CHANGES /RECOMMEND ATI OHS/ EXPLANATIONS 


QCd  DATA 

Forecast  Modola | 
-  GTS 


WORKSTATION  REVIEW 
(Color-Coded  Locations) 


i 


SMARTQC  WINDOW 


SMARTQC  RECOMMENDATION (S) 

ELAT  TYPE  ID  LAT  LONG  COURSE  MON  DAY/TIME 

SSMI  SATSSMI  61 . OSS  137. 19E  APR  08/1200 

D/S  LOCATION 
HIST  LOCATION 


OB 

FGS 

DIF 

BIAS 

STD  #  CORK  TRNS  BCF 

/ 

FC 

39.0 

22.0 

17.0 

95% 

3 

95% 

SSMI  WINDS  ARE  IN  AGREEMENT  WITH  ALTIMETER 
AND  SCATTEROMETER  WINDS 

ACCEPT  SSMI  WINDS 


FEEDBACK  TO  ALGORITHM  DEVELOPERS 
AND  MODELERS 


Monthly  Summary  of  Departures 

-  By  Regions  (latitudinal  bounds) 

-  Purged  by  QUIPS  Operator 

Extreme  Departures  (GE  to  +/-30  kts) 

-  By  Precip  Areas 


SSM/I  DERIVED  WIND  SPEED  DEPARTURES 
FROM  FIRST  GUESS  BY  LATITUDINAL  REGIONS 
November  5-19,  1991 


REGION 

NUMBER  OF  DEPARTURES 
(15  KTS) 

PERCENT  NEGATIVE 

GE  50  N 

1,085 

9.9 

31  N  -  50  N 

1,708 

37.3 

11  N  -  30  N 

381 

1.6 

0  N  -  10  N 

388 

0.0 

0  S  -  10  S 

353 

0.0 

11  S  -  30  S 

305 

0.7 

31  S  -  50  S 

1,623 

45.7 

GE  50  S 

1,043 

10.5 

TOTAL 

6,886 

23.6 

THE  6,886  WIND  SPEED  DEPARTURES  REPRESENTS  ONLY  1.5%  OF  THE 
MORE  THAN  HALF  MILLION  SSM/I  SUPER  OBS 
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STATUS  OF  NESDIS  DISTRIBUTION  OF  SSM/I  PRODUCTS 

P.  Taylor 

NOAA/NESDIS 
Washington,  DC  20233 


In  late  1991,  limited  testing  began  for  transmission  of  SSM/I  products  from  Fleet 
Numerical  Oceanographic  Center  (FNOC)  in  Monterey,  CA  to  NOAA/NESDIS  in 
Suitland,  MD  via  the  Shared  Processing  Network  (SPN).  These  products  generated  at  the 
SSM/I  Shared  Processing  Center  of  Expertise  (COE)  -  FNOC  -  include  Antenna 
Temperatures  (TDRs),  Brightness  Temperatures  (SDRs),  and  Derived  Environmental 
Products  (EDRs).  By  mid-summer  1992  NOAA/NESDIS  plans  to  1)  have  the  processing 
system  to  unpack  and  decode  these  products  in  place,  2)  have  these  transmissions 
validated,  and  3)  increase  the  frequency  of  transmissions  to  include  all  SSM/I  data  orbits 
from  two  operational  DMSP  spacecraft.  Throughout  the  past  year  experimental  generation 
by  NESDIS  of  some  of  the  SSM/I  EDR's  has  enabled  NOAA  users  to  analyze  these 
products  prior  to  their  routine  availability  from  FNOC.  Parallel  testing  of  the  products  has 
demonstrated  their  utility  for  weather  and  climate  applications  and  has  resulted  in  formal 
requests  from  the  National  Meteorological  Center,  and  National  Hurricane  Center,  and  the 
NAVY/NOAA  Joint  Ice  Center  for  operational  access  to  SSM/I  products.  In  order  to 
respond  to  these  requests  a  system  to  format  the  FNOC  products  for  NOAA  users,  provide 
science  support  to  the  Navy  for  product  enhancement,  and  distribute  the  products  is 
required.  This  presentation  will  review  the  preliminary  design  for  the  proposed  NESDIS 
SSM/I  Operational  System  for  distribution,  science  support,  and  archive  of  the  FNOC 
products. 
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PROCESSING  AND  DISTRIBUTING  SSM/I  DATA  VIA  THE  WETNET 

PROJECT 

H.  Michael  Goodman 


Earth  Science  and  Applications  Division 
NASA/Marshall  Space  Flight  Center 
Huntsville,  Alabama 


WetNet  is  a  NASA  sponsored  research  project  designed  to  promote  inter¬ 
disciplinary  research  in  the  atmospheric,  related  oceanic,  hydrologic  and  land  surface 
sciences.  The  project  currendy  provides  SSM/I  data  and  derived  data  products  to  over 
forty  U.S.  and  international  scientists.  The  project  is  more  than  a  data  service.  The  intent 
is  to  provide  common  data  sets,  analysis  software,  and  image  display  workstation  to  a 
diverse  groups  of  researchers  in  order  to  facilitate  cooperative  research  across  the  Earth 
sciences.  The  current  project  status  and  future  developments  will  be  presented.  The 
emphasis  will  be  on  data  processed  to  date,  problems  encountered  and  solutions  devised, 
data  distribution  media,  services  provided,  and  the  future  direction  of  the  project. 
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(VASA 


The 

WetNet  Program 


Marshall  Space  Flight  Center 
March  1992 


IWVSA 


WetNet 


WetNet  is  an  interdisciplinary,  interactive  data 
analysis  and  research  project  devoted  to  the  study  of 
the  global  hydrologic  cycle 


Wet:  studies  of  the  global  hydrologic  cycle 


Net:  develop  and  test  a  remote  interactive 
computer  network  in  an  Earth  system 
science  research  environment 


Marshall  Space  Flight  Center 
March  1992 


Page  1 
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IWSA 


Science  Objective 


To  promote  Earth  system  science  research 
through  cooperative  studies  of  the  global 
hydrologic  cycle. 


The  Special  Sensor  Microwave  /  Imager  (SSM/I) 
instrument  data  sets  are  provided  as  a  common 
data  set  from  which  interdisciplinary  working 
groups  can  conduct  passive  microwave 
investigations. 


Marshall  Spaca  Flight  Canlar 
March  1992 


Network  Objective 


To  utilize  existing  computer  communication 
networks  and  optical  disk  technology  to 
distribute  data  sets  to  a  diverse  group  of  Earth 
system  science  researchers. 


Marshall  Spaca  Flight  Cantar 
March  1992 
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Earth  Observing  System 

NASA 

•  EOS  Observatories  -  coordinated  international  suite  of 

instruments  and  platforms  for  long  term  remote 
sensing. 

•  EOS  Scientific  Research  Program  -  sustained  Earth  science 

research  commitment  for  the  purpose  of  data  analysis 
and  integration  into  predictive  models. 

•  EOS  Data  and  Information  System  (EOSDIS)  -  science 

driven  data  system  with  an  open  and  distributed 
architecture.  Learn  and  evolve  from  existing  data 
systems  and  prototypes. 

Marshall  Space  Flight  Cantar 
March  1992 


NASA 


WetNet  Organizational  Chart 

NASA  Earth  Science  and  Applications  Division 


NASA  I  HQ 


EOSOIS  Program  Office 
Dr.  Dixon  Butler 


Hydrologic  System  and  Radiation 
Program  Manager 


Dr.  James  Dodge 

Project  Scientist 

1 

NASA  1 MSFC 

Mr.  Michael  Goodman 

l 

Science 

Processing 

Applications 

Programming 

User  Services  | 

Mr.  Frank  LaFontaine 
Mr.  Don  Moss 

Mr.  Matt  Smith 

Ms.  Vada  LaFontaine 
Mr.  James  Dobbs  j 

Marshall  Space  Flight  Canter 
March  1992 
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WetNet  Research 
and  Data  Locations 


WetNet  Scientists 


Bob  Adler 
John  C.  Allshousa 
Richard  l_  Armstrong 
Susan  Avsry 
Eric  Barrstt 
Francis  Brathsrton 
Robert  A.  Brown 
Donald  J.  Cavalier) 
Alfred  T.  Chang 
Alaine  Chedln 
Bhaskar  J.  Choudhury 
Robert  G.  Crane 
James  Dodge 
Robert  N.  Eli 


NASA/Goddard  Space  Flight  Center 
NOAA/NESDIS 

National  Snow  and  Ice  Data  Center 
University  of  Colorado 
University  of  Bristol 
University  of  Wlsconsln-Madison 
University  of  Washington 
NASA/Goddard  Space  Flight  Center 
NASA/Goddard  Space  Flight  Center 
Ecole  Polytechnique 
NASA/Goddard  Space  Flight  Center 
The  Pennsylvania  State  University 
NASA  Headquarters 
West  Virginia  University 


Marshall  Space  Flight  Center 
March  1992 
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WetNet  Scientists 


William  Emery 
Gerry  W.  Felde 
Kevin  P.  Gallo 
Catherine  Gautier 
Steven  Goodman 
Kenneth  Hardy 
Jim  Hollinger 
Tony  Hollingsworth 

Greg  Hunott 
Nicole  Husson 
John  E.  Janowlak 
John  F.  LeMarshall 
Timothy  Uu 
James  P.  McGuIrk 


University  of  Colorado 
U.  S.  Air  Force  Geophysics  Laboratory 
National  Climatic  Data  Center 
University  of  Californla-Santa  Barbara 
NASA/Marshall  Space  Right  Center 
Lockheed  Mlsale  &  Space  Co.  Inc. 

Naval  Research  Laboratory 
European  Centre  for  Medium-Range 
Weather  Forecasts 
NASA  Headquarters 
Ecole  Polytechnique 

NOAA/NWS/NMC  -  Climate  Analysis  Center 
Australian  Bureau  of  Meteorology 
NASA/Jet  Propulsion  Laboratory 
Texas  A&M  University 


Mirahall  Space  Flight  Carter 
March  1992 _ 


NASA 


WetNet  Scientists 


Alberto  Mugnai 
Christopher  Neale 
William  S.  Olson 
Grant  Petty 
Franklin  R.  Robertson 
Barry  Rock 
Richard  C.  Savage 
Akira  Shibata 
Eric  A.  Smith 
Roy  Spencer 
Jeff  Star 

Graeme  L.  Stephens 
Fran  Stetlna 
Frank  J.  Wentz 


Instituto  Dl  Flslca  Dell'Atmosfera/CNR 
Utah  State  University 
University  of  Wlsconsln-Madlson 
Purdue  University 
NASA/Marshall  Space  Right  Center 
University  of  New  Hampshire 
Hughes  Aircraft  Co. 

Meteorological  Research  Institute 
Rorlda  State  University 
NASA/Marshall  Space  Right  Center 
University  of  Californla-Santa  Barbara 
Colorado  State  University 
NASA/Goddard  Space  Right  Center 
Remote  Sensing  Systems 


Marshall  Space  Flight  Canter 
March  1992 


Page  5 


113 


3/13/92 


vMllll)  Examples  of  Research 

CUAS l\ 


Savage  &  Heinrichs  (Hughes  Aircraft)  -  vegetation  classification  using  SSM/I  and 
comparisons  with  OLS  derived  vegetation  classification.  Stress  Identification 
to  distinguish  healthy  from  stressed  vegetation. 

Phallpou  (ECMWF)  •  Comparisons  of  global  preclpitable  water  with  ECMWF 
model  predictions.  Comparisons  of  radiosonde  moisture  calculations  with 
SSM/I  total  preclpitable  water. 

Smith,  at  al.  (Florida  State  U.)  -  Spatial  resolution  ehancament  of  19  &  37  GHz 
SSM/1  brightness  temperatures  through  deconvolution  procedures. 

Barrett,  et.  al.  (U.  of  Bristol,  U.K.)-  Global  precipitation  Irttercomparlson  project 
(PIP-1).  Demonstrate,  test,  intercompare,  and  validate  techniques  for  global 
monthly  rainfall  monitoring  by  SSM/1  .  Determine  which  techniques  function 
best  over  land,  ocean,  Ice  and  snow  covered  surfaces. 

ZIpser  &  McGuIrk  (Texas  A&M  U.)  -  Determination  of  SSM/I  usefulness  In 
categorizing,  measuring,  and  parLmeterlzing  the  effects  of  global  rainfall  with 
emphasis  on  the  effect  of  mesoscale  convective  systems. 

Marshall  Space  Flight  Cantor 
March  1992 


IUASA 


Data  Sources 


•  Utilize  space-based  measurements  from: 
.  SSM/I  (F8.F10&  F1 1) 

.  Geo  Satellites  (GOES,  Meteosat,  GMS) 


•  Utilize  ground-based  measurements  from: 
.  Radar  (USA  coverage)  (future) 

.  Surface  observations  (future) 

.  Rawinsonde  observations  (future) 


Marshall  Space  Flight  Canlar 
March  1992 
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Satellite  &  Instrument 
Characteristics 


•  Passive  microwave  radiometric  system 

•  Four  bi-polarlzed  frequencies 

V85.V37,  V22,  VI 9 
H85,  H37,  HI  9 

•  Spatial  Resolution 

85  GHz:  15x13  km 

37  GHz:  38x30  km 

22  GHz:  60  x  40  km 

19  GHz:  70x45  km 

•  Swath  Width:  1394  km 

•  Altitude:  833  km 

•  Scan  spacing:  12.5  km  (1.90  s) 

•  Orbits /day:  14.2 


Marshall  Space  Right  Canter 
March  1992 


SSM/I  Data  Availability 


NASA 


1987 

1988 

1989 

1990 

1991 

1992 


V85  bad 


H85  bad 


H85  good  i  F10 


1185  bad 


Marshall  Spacs  Flight  Csntsr 
March  1992 


NASA 


Data  Set  Composition 
SSM/I  and  Geostationary 


•  SSM/l  Daily  Browse  (Ascending  &  Descending) 

V85,  H85,  V37,  H37,  V22,  V19,  H19 

•  SSM/l  Full  Resolution  Swaths  (Ascending  &  Descending) 

V85,  H85,  V37,  H37,  V22,  V19,  H19 

•  SSM/l  Full  Resolution  Swath  Catalog 

•  GOES  Infrared 

GOES  East  at  00Z  and  12Z, 

GOES  West  at  02Z  and  14Z 


Marshall  Spaca  Flight  Cantor 
March  1992 


NASA 


Data  Set  Composition 
SSM/I  Products 

•  SSM/l  Products  Daily  Browse 

Rain  Rata  (Olson,  at.  al.) 

Total  Praclpltable  Water  (Allshouse,  et  al.) 

Cloud  Uquld  Water  (Aliahouse,  et  al.) 

Marins  Wind  Speed  (Goodberlet  et  al.) 

Land  Classification  (Neale,  et  al.) 

Land  Surface  Temperature  (McFarland,  et  al.) 

•  SSM/I  Products  Cartridge  Composite 

Rain  Rate  (Olson,  et.  al.) 

Total  Praclpltable  Water  (Allshouse,  et  al.) 

Cloud  Uquld  Water  (Allshouse,  et  al.) 

Marine  Wind  Speed  (Goodberlet  et  al.) 

Lend  Classification  (Neele,  et  al.) 

Land  Surface  Temperature  (McFarland,  et  el.) 
Preclpltaton  Identifier  (Spencer) 

Snow  Water  Equivalent  (Cheng) 

Vegetation  Index  (v37  -  h37) 

Sea  Ice  Fraction  (Ramseler,  et  al.) 

Sea  Ice  Age  (Ramseler,  et  al.) 

Marshall  Space  Flight  Center 
March  1992 
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Data  Set  Composition 
Applications  and  Ancillary  Software 

•  Software 

MclDAS-  OS2 

WotNet  Software 

User  Contributed  Software 

•  Documentation 

Help  Files 

System  Documentation 

•  Demonstrations 

Tutorials 


Marshall  Spaca  Right  Cantar 
March  1992 


( imp 

NASA 


WetNet  Data  Flow  Chart 


Marahall  Spaca  Flight  Can  tar 
March  1992 
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Research  Results  Investigator  Workshops 


Science  Investigator 


Marshal!  Space  Flight  Center 
March  1992  _ 


NASA 


Data  Management 
&  Distribution 


•  Develop  browse  and  catalog  algorithms  for 
managing  large  image  based  datasets. 

•  Explore  the  use  of  optical  disks  as  a  reliable 
storage  media  and  as  a  "high  bandwidth”  data 
transfer  mechanism. 

•  Use  existing  computer  networks  (where  feasible) 
to  distribute  daily  browse  datasets. 


Mar&tall  Spac*  Flight  Csntsf 


March  1992 


NASA 


Hardware  and  Software 
Specifications 


3/13/92 
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NAS/V 


WetNet  Workstation  Configuration 


Sony  Magneto-Optical  IBM  PS/2  Modal  70-A16 

650  MBytes  160  MByte  Hard  Disk 

6  MBytes  RAM 
3Com  or  Western  Digital 
Ethernet  card 

Marshall  Spaca  Right  Cantar 
March  1992 


Mous# 


NASA 


WetNet  Hardware 


•  System  Unit 

•  IBM  PS/  2  Model  70  -A1 6 

•  Intel  80386  Processor/80387  Math  Co-Processor 

•  6  MBytes  RAM 

•  160  MByte  Hard  Disk 

•  three  slots  (two  16-blt,  one  8-blt) 

•  Mouse 

•  Monitor 

•  VGA 


•  Magneto-Optical  Drive 

•  Sony  Rewritable /Erasable 

•  Adaptec  SCSI  Adapter 


Ethernet  Adapter 

•  3COM  or  Western  Digital 


Marshall  Spaca  Flight  Cantar 
March  1992 


3/13/92 
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NASA 


System  Software 


•  Operating  System 

a  IBM  Operating  System/2 

Standard  Edition  (SE)  vl.3  or  v2.0 

•  MclDAS 

•  MclDAS-OS2  vS.5 

•  TCP/IP 

.  IBM  TCP/IP  for  OS/2  vl.2 


Marshall  Space  Flight  Canter 
March  1992 


NASA 


McIDAS-OS2  Programming 


What  you  need 

.  Microsoft  FORTRAN  5.1 
a  Applications  Programming  Manual 
a  Knowledge  of  FORTRAN 
a  MclDAS-OS2  source  code 


Document  and  share  your  code 


Marshal  Space  Right  Canter 
March  1992 
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Session  3:  New  Retrieval  Approaches  and  Products  and  Relevance  to  AMSU 

(Ohring,  chair) 

In  this  session,  research  results  were  presented  on  the  development  of  an  advanced 
satellite  retrieval  system,  new  retrieval  approaches,  and  the  use  of  DMSP  SSM/T  and 
SSM/I  data  to  develop  retrieval  methodologies  for  NOAA’s  Advanced  Microwave 
Sounding  Unit-A  (AMSU-A).  AMSU-A,  which  is  scheduled  for  launch  on  the  NOAA-K 
satellite  in  1995,  is  a  15-channel  temperature  sounder  with  many  channels  similar  to  those 
on  SSM/T  and  SSM/I. 

W.  Baker  of  NOAA/NMC  discussed  a  joint  NESDIS/NASA/NMC  research 
program  to  develop  an  advanced  satellite  retrieval  system  that  is  an  integral  part  of  NMCs 
numerical  weather  analysis  and  prediction  model.  The  new  system  incorporates  the  best 
components  of  the  operational  retrieval  system  at  NESDIS  and  the  experimental  approach  at 
NASA.  A  key  aspect  of  the  new  system  is  the  use  of  the  6-hour  forecast  temperature 
profiles  as  first  guesses  in  the  satellite  retrieval  process.  Test  results  are  encouraging  and 
the  system  is  scheduled  for  operational  implementation  during  1992. 

L.  McMillin  of  the  NOAA  Satellite  Research  Laboratory  reported  on  his  work  in  the 
use  of  a  classification  approach  to  provide  an  initial  temperature  profile  for  the  retrieval 
process.  Application  of  the  method  to  TO  VS  data  suggests  an  increase  in  accuracy  of  0.2- 
0.3°C  compared  to  the  use  of  a  library  search  first  juess.  Similar  results  are  obtained  when 
the  method  is  applied  to  the  microwave-only  observations  of  the  SSM/T. 

In  a  companion  pair  of  papers,  H.  Fleming  of  the  NOAA  Satellite  Research 
Laboratory  and  E.  Kratz  of  SM  Systems  and  Research  Corporation  discussed  the 
derivation  and  application  of  a  physical  retrieval  algorithm  for  SSM/T.  The  method  is  a 
two-part  retrieval  procedure;  the  first  part  produces  an  atmospheric  temperature  profile  and 
the  second  part  produces  surface  (skin)  temperature  and  surface  emissivity  at  50  GHz. 
Comparisons  of  retrieved  temperatures  with  coincident  radiosonde  temperature  profiles 
indicate  that  the  accuracy  of  this  physical  approach  is  comparable  to  the  operational, 
regression  based,  retrieval  algorithm,  but  produces  the  two  surface  parameters  which  the 
regression  approach  cannot  do. 

J.I.F.  King  of  the  USAF  Phillips  Laboratory  described  the  application  of  the  zeta- 
funcdon  transform  to  the  satellite  temperature  retrieval  problem.  In  conventional 
approaches,  the  radiance  integral  is  expressed  as  a  finite  sum  quadrature  with  the  matrix 
then  inverted  to  yield  the  temperature.  In  the  proposed  approach,  a  moment  method  is  used 
to  approximate  the  temperature  profiles  at  three  interacting  bilateral  wave  packets. 
Application  of  the  zeta-function  transform  to  the  three  packets  reduces  the  complex 
interaction  of  the  radiance  and  temperature  profiles  to  a  simple  subtractive  operation. 
Examples  illustrated  the  quantitative  relationship  between  radiance  observation  noise  and 
the  retrieval  accuracy. 

P.W.  Ros  akranz  of  M.I.T.  discussed  SSM/I  measurements  as  predictors  of  the 
response  of  AMSU-A  to  surface  and  atmospheric  phenomena.  SSM/I’s  hannels  (19, 22, 
37,  and  85  GHz;  are  similar  to  the  AMSU-A  window  channels  (24,  31, 50,  and  89  GHz). 
Differences  in  in-: :•  Vice  angle  and  polarization  between  the  two  instruments  are  reconciled 
by  reasonable  ass  ptiors  about  angle  and  polarization  dependence.  The  SSM/I  data 
show  signatures  ot  rough  seas,  sea  ice,  snow  cover,  and  precipitation  that  are  relevant  to 
the  interpretation  of  the  measurements  to  be  made  by  AMSU-A. 
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SESSION  3  -  NEW  RETRIEVAL  APPROACHES  AND  PRODUCTS; 
RELEVANCE  TO  AMSU 


JOINT  NESDIS/NASA/NMC  EFFORT  TO  DEVELOP  AN  ADVANCED 
SATELLITE  RETRIEVAL  SYSTEM 

Wayman  Baker 

National  Meteorological  Center 
Washington,  DC  20233 


1 .  Introduction 

The  accuracy  of  operational  numerical  weather  forecasts  has  steadily  increased 
during  the  past  decade  to  the  point  where  it  is  becoming  difficult  to  routinely  expect 
operational  satellite  temperature  soundings  to  have  a  positive  impact  on  the  forecasts  over 
the  Northern  Hemisphere.  The  realization  of  the  importance  of  the  satellite  retrieval 
problem  to  numerical  weather  prediction  (NWP)  culminated  in  the  convening  of  a  U.S. 
government  Interagency  Satellite  Retrieval  Workshop  in  July  1989  at  NMC  which  included 
satellite  retrieval  and  NWP  experts  from  NASA  and  NOAA.  The  Workshop  participants 
agreed  that  an  advanced  satellite  retrieval  system  for  numerical  weather  forecasting  and 
climate  monitoring  should  be  developed  as  soon  as  possible.  This  system  should  take 
advantage  of  the  best  components  of  the  operational  system  at  NESDIS  and  those  from  the 
experimental  approaches  at  NESDIS  and  at  NASA.  There  was  a  strong  consensus  among 
the  Workshop  participants,  that,  over  most  of  the  globe,  the  present-day  accuracy  of  the 
model  6-hour  forecast  during  the  data  assimilation  should  be  beneficial  to  the  retrieval 
process.  The  approach  of  producing  and  using  the  satellite  soundings  "interactively"  then 
naturally  follows,  where  the  assimilation  model  forecast  is  used  as  a  first  guess  by  the 
retrieval  scheme  (Baker  et  al.,  1984).  The  resulting  retrievals  are  then  analyzed  in 
conjunction  with  other  available  data  to  produce  initial  conditions  for  the  subsequent 
forecast  by  the  assimilating  model.  Such  an  approach,  developed  and  used  extensively  in  a 
research  mode  by  J.  Susskind  and  collaborators  at  the  NASA  Goddard  Space  Flight 
Center,  has  been  shown  to  have  a  beneficial  effect  on  both  the  accuracy  of  the  6-hour 
assimilation  model  forecast  and  that  of  the  retrievals.  More  importantly,  the  negative 
impacts  on  medium  range  forecasts  over  the  Northern  Hemisphere  are,  for  the  most  part, 
eliminated  (Susskind  and  Pfaendtner,  1989).  In  addition  to  temperature  profiles  and 
precipitable  water  estimates,  a  wide  range  of  geophysical  parameters  (e.g.,  land  and  sea 
surface  temperature,  ice  and  snow  extent,  cloud  top  pressure  and  amount,  and  estimates  of 
outgoing  long  wave  radiation,  soil  moisture,  total  ozone,  and  precipitation)  can  be 
produced.  A  particularly  strong  point  with  the  interactive  approach  is  that  the  above 
parameters  are  obtained  in  a  way  which  is  consistent  with  the  mass  and  wind  fields  in  the 
assimilation  because  the  6-hour  forecast  from  the  assimilation  is  used  as  the  first  guess  for 
the  retrievals.  Therefore,  not  only  is  an  accurate,  self-consistent  set  of  data  produced  for 
numerical  weather  forecasting,  a  powerful  climate  monitoring  system  is  also  created. 

2.  Proposed  Strategy 

There  was  a  general  agreement  among  the  Workshop  participants  that  the  various 
sophisticated  inversion  methods  (e.g.,  the  NESDIS  operational  system  (Fleming  et  al., 
1986b,  1988);  the  NASA  research  scheme  (Susskind  et  al.,  1984))  should  have  about  the 
same  accuracy,  given  the  same  radiances.  There  was  also  a  strong  consensus  that  the  most 
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accurate  first  guess  possible  should  be  used  for  the  inversion  problem.  Over  most  areas  of 
the  globe  that  is  most  likely  a  forecast  provided  by  a  general  circulation  model  during  the 
data  assimilation.  In  areas  where  the  model  forecast  might  be  poor  as  determined  by 
comparing  observed  (satellite)  and  computed  (model)  radiances  from  the  forward  problem, 
a  classification  approach  (i.e.,  McMillin,  1986b)  blended  with  the  model  first  guess,  seems 
promising. 

The  Workshop  participants  strongly  agreed  that  the  most  important  remaining 
weakness  (assuming  the  use  of  a  forecast  first  guess)  in  the  NESDIS  operational  system 
was  in  the  approach  used  for  cloud-clearing,  where  an  angle  correction  (to  nadir)  is 
performed  before  cloud-cle?  mg,  and  hence,  the  angle  correction  procedure  must  also 
simultaneously  account  foi  -louds.  All  the  participants  agreed  that  the  most  accurate  off- 
nadir  retrievals  are  produced  if  the  off-nadir  radiances  are  cloud-cleared  at  the  same  angle 
that  the  satellite  views  the  atmosphere  (as  in  the  NASA  scheme).  In  addition,  the  NASA 
scheme  also  uses  the  forecast  temperature  information  in  the  cloud-clearing,  while  the 
NESDIS  scheme  does  not  use  this  information.  On  the  other  hand,  a  limitation  identified  in 
the  NASA  retrieval  scheme  was  in  not  utilizing  the  SSU  data  for  the  stratospheric  sounding 
problem,  as  is  done  in  the  NESDIS  system. 

Various  aspects  of  the  NESDIS  and  GLA  retrieval  schemes  will  be  intercompared 
and  tested.  Retrievals  from  a  combined  system  will  be  evaluated  using  collocated 
radiosondes,  and  by  conducting  forecast  impact  tests  with  the  NMC  global  data 
assimilation  system  (GDAS).  Additional  forecast  impact  tests  will  then  be  conducted  in 
parallel  to  the  operational  system  in  real  time  with  the  combined  retrieval  scheme. 

With  the  above  considerations  in  mind,  the  Workshop  participants  recommended 
the  following  steps  in  order  to  combine  the  strengths  of  the  different  approaches: 

1 .  NESDIS  and  NASA  will  exchange  retrieval  software. 

2.  Both  retrieval  systems  will  be  implemented  on  the  NMC  CRAY  YMP  for 
experimentation. 

3 .  The  NASA  approach  for  cloud-clearing  (producing  cloud-cleared  radiances 
at  the  angles  of  observation)  will  be  tested  by  NESDIS  (but  without  iteration 
as  in  the  NASA  system). 

4.  NASA  will  provide  NESDIS  cloud-cleared  radiances;  NESDIS  will  produce 
retrievals;  NMC  will  evaluate  the  analysis/forecast  impact. 

5.  *'  Vt  \  wii!  conduct  tests  to  determine  the  sensitivity  of  its  cloud-clearing 

Ix-me  to  the  number  of  iterations.  Cloud-cleared  radiances  computed  from 
the  f,  ue  s,  after  one  iteration,  and  after  the  full  adjustment  made  in  the 
NY  \  system,  will  be  provided  to  NESDIS  for  comparison. 

6.  NESDIS  will  implement  their  bias  correction  (on  observed  radiances  during 
the  forward  problem)  in  the  NASA  code.  This  will  be  compared  with  the 
NASA  approa  ;h. 

7.  Other  aspects  of  the  retrieval  problem  will  also  be  evaluated,  such  as  the 
ability  to  retain  an  accurate  first  guess  and  improve  a  poor  one. 

The  initial  thrusts  of  the  joint  project  involve:  1)  exporting  the  NASA  approach  for 
cloud-clearing  to  NESDIS  for  evaluation,  and,  2)  developing  the  capability  to  combine  ar 
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air-mass  classification  guess  with  the  NMC  6  hour  forecast  guess  in  regions  where  the 
forward  calculation  differences  are  large. 

3 .  NESDIS/NMC  Pilot  Study 

A  joint  NESDIS/NMC  pilot  effort  was  initiated  in  early  1989  to  develop  an 
interactive  retrieval  system. 

Phase  II  of  the  pilot  effort  is  now  underway  which  involves  developing  the 
capability  to  retrieve  temperature  and  specific  humidity  at  the  NMC  model  levels.  This  has 
two  advantages:  1 )  the  error  due  to  vertical  interpolation  of  the  model  6-hour  forecast  to  the 
levels  where  the  retrievals  are  produced  (previously,  constant  pressure  levels)  is  eliminated, 
and,  2)  the  new  NMC  global  spectral  interpolation  (SSI)  objective  analysis  scheme  (Parrish 
and  Derber,  1991),  which  was  implemented  on  June  25,  1991,  performs  the  analysis  at  the 
model  levels,  which  also  eliminates  a  source  of  error. 

The  retrieval  code  has  now  been  ported  to  the  CRAY  YMP  and  parallel 
analysis/forecast  experiments  conducted  at  T62  resolution.  A  full  resolution  (T126)  impact 
test  will  be  conducted  soon,  with  operational  implementation  anticipated  in  early  1992. 
Enhancements  to  the  pilot  system,  as  discussed  above  (e.g.,  the  hybrid  model  first  guess  - 
classification  approach,  use  of  the  NASA  cloud  clearing,  etc.),  will  be  implemented 
following  successful  parallel  analysis/forecast  experiments  on  the  NMC  CRAY  YMP. 
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500  MB  RMS  TEMP  ERROR  FOR  NORTHERN  HEMISPHERE: 

6-HOUR  FORECAST  VS  SATELLITE  RETRIEVAL 
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INTERACTIVE  ANALYSIS/FORECAST/RETRIFVAI  SYSTEM 


Fig.  1  Schematic  diagram  of  an  interactive  analysis/forecast/ 
retrieval  system  under  joint  development  by  NESDIS  and 
NIK. 


RETRIEVAL  -  RAOB  BIAS  ERROR 

WITH  AND  W/O  BIAS  CORRECTION 


SIGMA*  1000 

UNCORRECTED  BTEMPS  CORRECTED  BTEMPS 
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Concluding  Remarks 


Implementation  of  the  interactive  retrievals 
represents  a  significant  enhancement  in  our 
ability  to  fully  utilize  satellite  radiances  for 
NWP.  However,  further  work  remains  in 
several  areas: 

-  Air-mass  classification 

-  Cloud-clearing 

-  Use  of  error  covariance  matrix  at  each 
gridpoint 

-  Enhanced  water  vapor  retrievals 

-  Other  retrieved  parameters  (i.e., 
precipitation) 
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CLASSIFICATION  RETRIEVAL  APPROACHES  FOR  DMSP 


Larry  McMillin 


National  Oceanic  and  Atmospheric  Administration 
National  Environmental  Satellite,  Data,  and  Information  Service 
Satellite  Research  Laboratory 
Washington,  D.C.  20233 


Classification  methods  have  been  used  to  derive  temperature  and  water  vapor 
profiles  from  satellite  radiances  measured  by  the  TIROS  Operational  Vertical  Sounder 
(TOVS).  When  compared  to  approaches  that  start  with  a  climatological  first  guess,  the 
classification  approaches  have  demonstrated  an  ability  to  recover  more  of  the  information 
available  in  the  radiances  and  increase  the  accuracy  of  the  resulting  retrievals  by  0.2  to 
0.3  K  rms  (see  McMillin  1991).  The  empirical  classification  regression  described  by 
McMillin  has  been  successfully  applied  to  microwave  data  from  the  MSU  sensor  with 
similar  results.  At  this  meeting,  results  of  the  method  using  SSMT  data  will  be  presented. 
In  addition,  comparisons  of  this  classification  approach  with  alternative  classification 
approaches  will  be  discussed. 

McMillin,  L.M.,  1991:  Evaluation  of  a  Gassification  Method  for  Retrieving  Atmospheric 
Temperatures  from  Satellite  Measurements,  J.  Appl.  Meteor.,  30,  432-446. 
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Pressure  (hPo) 


Classification  Retrieval  Approaches  for  DMSP 


Larry  McMillin 


National  Environmental  Satellite,  Data,  and  Information  Service 
Satellite  Research  Laboratory 
Washington,  DC 


Retrieval  Statistics  (Raob  -  RTV)  for:  noaalO,  ALL 


Bias  ond  SO  errors  (K,  g/Kg) 


l 
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Overview  of  Method 


1.  Calculate  an  average  temperature  for  the  profile 

2.  Adjust  the  measured  radiances 

3.  Calculate  the  principal  component  scores  of  the  eigenvectors 

4.  Calculate  the  class  index 

5.  Get  the  class  mean  temperature  and  radiance 

6.  Apply  the  regression  coefficients 


Calculate  the  average  temperature  for  the  profile 

1.  Temperature  is  unknown  so  average  radiances 

2.  Average  channels  1,2,3,5,6,10 

3.  Denote  the  sample  average  by  "a" 


Adjust  the  Measured  Radiances 

1.  Subtract  the  profile  average  from  all  radiances 

2.  V*  =  T*  -  a  -  V* 

T*  IS  THE  VECTOR  OF  RADIANCE  TEMPERATURES 
V*  IS  THE  AVERAGE  FOR  THE  DEPENDENT  SAMPLE 
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Calculate  the  Principal  Component  Scores 

1.  Multiply  the  adjusted  radiances  by  the  eigenvector  coefficients  for 
CHANNELS  1-6  AND  10. 

2.  a  =  EV* 

a  ARE  PRINCIPAL  COMPONENT  SCORES 
E  ARE  THE  EIGENVECTOR  COEFFICIENTS 


Calculate  the  Class  Index 

1.  Use  4  EIGENVECTORS 

2.  14  DIVISIONS  FOR  THE  1ST 

3.  4  FOR  THE  SECOND 

4.  2  FOR  THE  THIRD 

5.  2  FOR  THE  FOURTH 


Get  the  Class  Mean  Temperature  and  Radiance 
1.  The  class  means  are  generated  from  a  sample  of  collocations. 

T(c)  ARE  THE  MEAN  TEMPERATURES  FOR  THE  CLASS 

T* (C)  ARE  THE  MEAN  BRIGHTNESS  TEMPERATURES  FOR  THE  CLASS 
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STANDARD  DEVIATION  (k) 


100  200  300  400  500  600  700  800  900  1000 

PRESSURE  (mb) 


OPERATIONAL  CLASS  -  REGRESSION 


Conclusions 


1.  A  CLASSIFICATION  -  hEGtffeSSION  SYSTEM  PRODUCES  A  SIGNIFICANT  (SEVERAL 
TENTHS  OF  A  OEGREEJ  INCREASE  IN  ACCURACY  OVER  THE  CURRENT  OPERATION. 
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USING  SSM/T  DATA  AS  PROOF  OF  CONCEPT  FOR  AN 
AMSU-A  RETRIEVAL  ALGORITHM 


Henry  E.  Fleming 

NOAA/NESDIS 
Satellite  Research  Laboratory 
World  Weather  Building,  Room  810 
Washington,  DC  20233 

Eugene  J.  Kratz 

SM  Systems  &  Research  Corporation 
8401  Corporate  Drive,  Suite  510 
Landover,  MD  20785 


Brightness  temperature  data  from  DMSP/SSM/T  have  been  used  as  proof  of 
concept  for  a  proposed  AMSU-A  operational  retrieval  algorithm.  The  algorithm  is  for  the 
simultaneous  retrieval  of  the  atmospheric  temperature  profile,  the  surface  temperature,  and 
the  surface  emissivity.  The  solution  is  based  on  the  minimum  variance  method,  which 
minimizes  the  expected  mean-square  error  for  a  linear  unbiased  estimate.  Since  the 
brightness  temperature  observations  sense  these  meteorological  parameters  simultaneously, 
retrieving  them  simultaneously  guarantees  that  the  solution  satisfies  the  radiance. 

Surface  temperature  and  emissivity  occur  as  a  product  in  the  boundary  term  of  the 
radiative  transfer  equation,  therefore,  this  product  is  retrieved  as  a  single  parameter,  and 
then  typically  the  surface  temperature  is  obtained  from  some  independent  source  and 
divided  into  the  retrieved  product  to  extract  the  surface  emissivity.  A  novel  feature  of  our 
method  is  that  we  retrieve  these  two  parameters  separately  using  only  the  brightness 
temperature  measurements,  i.e.,  independent  knowledge  of  the  surface  temperature  is  not 
required.  The  deri  'ation  of  this  procedure,  based  on  the  minimum  variance  approach,  will 
be  presented.  Results  will  be  given  to  demonstrate  the  validity  of  the  retrieval  algorithm. 
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USING  SSM/T  DATA  AS  PROOF  OF  CONCEPT 
FOR  AN  AMSU-A  RETRIEVAL  ALGORITHM 


Henry  E.  Fleming 
Satellite  Research  Laboratory 
NOAA/NESDIS 

Eugene  J .  Kratz 

SM  Systems  &  Research  Corporation 


1 

Ti  =  (e*ts)  Ti(ps)  +  j  t(p)  dii^e) 

rs 

=  ith  chan,  brightness  temperature 
e  =  surface  emissivity  g  50.5  GHz 

t(p)  =  atmospheric  temperature  profile 
ts  =  surface  (skin)  temperature 

J"i(P)  =  ith  chan,  standard  transmittance 
r i (P *  € )  =  effective  transmittance,  where 

^i(P#€)  =  {1  -  [1  -€] (ri(ps)/Ti(p) )l  )rj(p) 
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MANIPULATION  OF  THE  RADIATIVE  TRANSFER  EO . 

1.  Assume  the  initial  guess  r  ^  r 

2.  Subtract  the  initial  state  RTE  from  the 
exact  RTE 

3.  Let  6x  =  x  -  "x  for  any  quantity  x 

4.  Set  €  =  7  +  5e  in  dri(p,e) 

5.  Reflectivity  r  =  1  -  e,  hence 

-Se  =  Sr 

1 

6Ti  =  6{e ts)  TjJPs)  +  |  <5t(p)  dTjJp,?) 

rs 

1 

-  «5e  J  t(p)  [Tifpgj/Tifp)  ]*  dr  i  (p) 


1 

t*i  =  -J  t(p)  { [Ti(ps)/ri(p) y  dT± (p) } 

f  c: 
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SSM/T  REFLECTED  COMPONENT 
WEIGHTING  FUNCTIONS 


PRESSURE  (MB) 


<5T  =  &e  t*  +  6  ( ets)  rs  +  A  <5t 


PART  I 


Let  e  =  e ,  then 


«5T  =  6  ( ets)  rs  +  A  5t 


PART  II 


Let  ST  -  A  St  =  SS,  then 


SS  =  6e  t*  +  <5(ets)  ts 

which  is  an  overdetermined  3x2  system. 
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TEMPERATURE  AND  EMISSIVITY  RESULTS  FROM  A 
SSM/T  RETRIEVAL  ALGORITHM 

Eugene  J.  Kratz 

SM  Systems  &  Research  Corporation 
8401  Corporate  Drive,  Suite  510 
Landover,  MD  20785 

Henry  E.  Fleming 

NOAA/NESDIS 
Satellite  Research  Laboratory 
World  Weather  Building,  Room  810 
Washington,  DC  20233 


A  physical  retrieval  algorithm  has  been  developed  utilizing  DMSP  SSM/T 
brightness  temperature  data.  The  retrieved  products  include  atmospheric  temperature 
profiles,  the  surface  (skin)  temperature,  and  the  surface  emissivity  at  50  GHz.  Only  a  brief 
description  of  the  algorithm  will  be  given  because  the  emphasis  of  the  talk  will  be  on  the 
accuracy  of  the  retrieved  products.  Not  only  will  the  retrievals  be  compared  with 
coincident  radiosonde  temperature  profiles,  but  also  with  the  NESDIS  Operational  Shared 
Processing  profiles.  The  Shared  Processing  algorithm  uses  a  linear  regression  approach. 
The  physical  retrieval  algorithm  is  shown  to  be  at  least  competitive  with  the  operational 
retrieval  algorithm,  and  in  some  cases  exceeds  the  accuracy  of  the  Shared  Processing 
results.  A  comparative  analysis  of  the  results  from  the  two  retrieval  approaches  will  be 
given. 
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TEMPERATURE  AND  EMISSIVITY  RESULTS 


FROM  A  SSM/T  RETRIEVAL  ALGORITHM 


Eugene  J.  Kratz 

SM  Systems  &  Research  Corporation 


Henry  E.  Fleming 
Satellite  Research  Laboratory 
NOAA/NESDIS 


PhYSICAL  vs.  OP.  REGRESSION  RETRIEVAL  ERRORS  (596) 
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measured  sst 


RETRIEVED  VS  MEASURED  SST 


INITIAL  GUESS  AND  RETRIEVED  EMISSIVITIES 


Emiss i v i ty 


PAINLESS  EXTRACTION:  SUBTRACTIVE  TEMPERATURE  SENSING 
FROM  SATELLITE  RADIANCES  USING  THE  ZETA  TRANSFORM 

Jean  I.F.  King 

Phillips  Laboratory 
Geophysics  Directorate 
Hanscom  AFB,  MA  01731-5000 


The  radiance  integral  specifies  the  satellite-observed  flux  R(z)  as  a  complex  folding 
or  convolution  of  the  Planck  temperature  profile  B(z)  with  the  weighting  function  w  (z  -  z) 
which  represents  the  medium  interaction.  The  temperature  profile  is  then  identified  with  the 
solution  of  a  Fredholm  equation  of  the  first  kind.  In  conventional  approaches  a 
discretization  process  expresses  this  equation  as  a  finite  sum  quadrature  with  the  matrix 
then  inverted  to  yield  the  Planck  intensity. 

We  pursue  a  different  tack.  A  moment  method  is  used  to  approximate  the  profiles 
as  three  interacting  bilaterial  wave  packets.  A  zeta-function  transform  applied  to  the  three 
packets  reduces  the  complex  interaction  of  the  radiance  and  temperature  profiles  to  a 
subtractive  arithmetic  operation.  By  eliminating  the  inversion  operation  as  a  source  of 
error,  we  are  positioned  to  determine  the  limits  of  retrieval  accuracy  imposed  by  the 
inevitable  noise  in  the  radiance  observations  and  their  number  limitations.  Examples  will 
be  given  which  quantitatively  specify  this  noise  vs.  accuracy  relationship. 
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ZETA  STATE  VECTOR  ANALYSIS 


TEMPERATURE  INFERENCING  REDUCED 

FROM  MATRIX  INVERSION  TO  COLUMN  VECTOR 

ARITHMETIC 


DR.  JEAN  I.F.  KING 

PHILLIPS  LABORATORY/GEOPHYSICS  DIRECTORATE 
t#tC/NESDIS/DOD  CONFERENCE  ON 
DMSP  RETRIEVAL  PRODUCTS 
15  April  1992 
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TRICOLOR  ZETA  NUMBER  THEOREM:  THE  KEY 

ZETA  NUMBER  SET 
1  N  1 

=  n  ,  ?n  =  rn  (n  =  1.2....N;  N  =  2,3,...) 

CORRESPONDING  DENSITY  PROFILE 
Rn(2)  =  R-(z)  +  R  +  (z) 

R_(z)  =  H(-z)V-^^  (Re  a;  <  O  ) 

i=1  rn  I  ai ) 

R+( 2  )  =  -H  ( z  )  .2  r^:  (Re  a,-  >  O  ) 

i  =m  rn  l  ai  / 

FROM  DENSITY  PROFILE  TO  ZETA  SET 

R  (z)^xn  (moments  )-*£n  =  rn  Tn  =  1 

n! 

FROM  ZETA  SET  TO  DENSITY  PROFILE 
i'n  =  aj->Rn  (  z) 

SIMILARLY  FOR  gn  =  Cn  ,  Gn  AND  bn  =  Cn  ,  Bn  PAIRS 


J  zn  R  (  z)  dz 


i_^rT  ^  =  c5  (t, )  = 


Tj  =  (X, ) 


<7^5  (5  =  S5(X,)  = 


T,3  T,5  T2  T22 
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RMS  Delta  Jejnperat.ure 


PHYSICALLY  RETRIEVED  TEMPERATURE  ACCURACY 
VS.  NO.  OF  RADIANCE  MOMENTS  (ZETA  TRANSFORM  METHOD) 


RADIANCE  MOMENTS 
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Pressure  (mb) 


Weighting  Function  Dependence  on  Interred  Profile 
Elsasser  Weighting  Function  (m=0.5) 


^B  “  ^R 

-  (-1)"  Co 

Radiance 

Moment 

0 

0.00000 

0.00000 

0.00000 

1 

0.00000 

0.63518 

0.63518 

2 

1.64493 

0.61685 

1.02808 

3 

0.00000 

0.35060 

0.35060 

4 

0.54116 

0.25367 

0.28749 

5 

0.00000 

0.20090 

0.20090 

6 

0.33911 

0.16691 

0.17220 

7 

0.00000 

0.14292 

0.14292 

8 

0.25102 

0.12502 

0.12600 

9 

0.00000 

0.11112 

0.11112 

10 

0.20020 

0.10000 

0.10020 

Retrieval  Error  Dependence 
on  Weighting  Function  Specification 


Temperature  (K) 
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HMS  Prlta  IVfnprratuf 


Retrieval  Error  Dependence 
on  Weighting  Function  Specification 


ZETA  VECTOR  ANALYSIS  MAKES  POSSIBLE 

•  REDUCTION  OF  INVERSION  TO  COLUMN  ARITHMETIC 

•  SOLUTION  OF  BLACK  BOX  PROBLEM 

•  ANALYSIS  OF  RETRIEVAL  SENSITIVITY  AND  ACCURACY 

•  INTERACTIVE  INVERSION 
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SSM/I  MEASUREMENTS  AS  PREDICTORS  OF  THE  RESPONSE  OF 
AMSU-A  TO  SURFACE  AND  ATMOSPHERIC  PHENOMENA 

P.W.  Rosenkranz 

Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 


SSM/I  measures  at  four  frequencies,  19, 22,  37  and  85  GHz,  which  covers  a  range 
similar  to  the  24, 31, 50  and  89  GHz  window  channels  on  AMSU-A.  There  are 
differences  in  incidence  angle  (constant  53°  for  SSM/I;  varying  for  AMSU-A)  and 
polarization  (V  and  H,  except  V  only  at  22  GHz,  for  SSM/I;  mixed  linear,  depending  on 
scan  position  for  AMSU-A),  but  these  differences  can  be  reconciled  with  the  aid  of  some 
reasonable  assumptions  about  angle  and  polarization  dependence.  Thus,  SSM/I  data  yields 
information  about  signatures  of  surface  phenomena  such  as  rough  seas,  ice  and  snow,  and 
atmospheric  phenomena  such  as  precipitation,  that  will  be  relevant  to  the  interpretation  of 
measurements  to  be  made  by  AMSU-A. 
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SMI 


AMSU-A 


WINDOW  FREQUENCIES 

19/  22/  37/  86  GHz 

2 ,/  31/  50/  89  GHz 

POLARIZATIONS 

V/  H  (EXCEPT  22) 

ROTATES  FROM  -48*  TO  +48 

INCIDENCE  ANGLE 

53* 

VARIES  FROM  2*  TO  57° 

Window  channels  - 

Radiative  transfer  model  (non-scattering): 

SURFACE  AT  Tg  WITH  EMISSIVITY  e 
ISOTHERMAL  ATMOSPHERE  AT  Ta 
T  SCALE  DEFINED  SO  TC0SMIC  *  0 

tB  ■  \  *  t2  <'TS  -  V 

Transmittance  of  the  atmosphere: 

=  EXP(-(a(v,Ta)  +  B(v,Ta)V  +  Y(v/Ta)L)/u) 

V  =  H20  vapor  burden 
L  =  LIQUID  Ho0  BURDEN 
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O  Hollinger  (1.41  GHz) 

0  Swift  (4.0  GHz) 

+  Sasaki  ei  al.  (6  7  GHz) 
A  Swift  (7.5  GHz) 

•  Hollinger  (8.36  GHz) 
x  Sasaki  et  al.  (18.6  GHz) 
A  Hollinger  (19.34  GHz) 
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Incidence  angle  (degrees) 

Fig.  4.  Dependence  of  wind  speed  sensitivity  of  horizontally  polarized 
brightness  temperature  at  different  frequencies  on  incidence  angle. 

O  Hollinger  (1.41  GHz) 

V  Swift  (1.4  GHz) 

0  Swift  (4.0  GHz) 

+  Sasaki  el  at.  (6.7  GHz) 

▼  Swift  (7.5  GHz) 

•  Hollinger  (8.36  GHz) 
x  Sasaki  ei  al.  (18.6  GHz) 

A  Hollinger  (19.34  GHz) 
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Session  4:  Applications  to  Numerical  Weather  Prediction  (Kalnay,  chair) 


In  this  session  several  operational  and  research  atmospheric  data  assimilation  and 
forecast  systems  using  DMSP  products  were  described.  In  the  last  part  of  the  session 
simple  model  applications  of  data  to  retrieve  additional  parameters  were  presented. 

The  first  paper  was  presented  by  Pat  Phoebus  from  NOARL.  She  indicated  that  the 
Navy  is  most  advanced  in  the  operational  assimilation  of  DMSP  environmental  products: 
SSM/T  atmospheric  soundings,  SSMI  surface  wind  speeds,  and  water  are  assimilated  in 
the  Navy  system.  They  are  also  experimenting  with  the  assimilation  of  SSM/I  total 
precipitable  water  by  creating  moisture  profiles  based  on  the  forecast  profiles  modified  to 
contain  the  same  total  as  the  SSM/I  observation.  The  present  lack  of  computational 
resources  to  make  parallel  integrations  has  not  allowed  a  rigorous  testing  of  whether  these 
data  have  a  positive  impact  upon  the  forecasts,  but  indirect  evidence  suggests  that  at  least 
the  surface  wind  speeds  are  beneficial. 

Ken  Mitchell,  from  NMC,  reviewed  the  current  research  on  improving  the 
sophistication  of  the  surface  parameterizations  in  the  NMC  models,  which  in  turn  makes 
them  more  sensitive  to  the  accurate  definition  of  the  current  state  of  surface  parameters  such 
as  ice  and  snow  cover  and  soil  moisture,  as  well  as  cloud  cover.  He  showed  that  the 
manually  derived  NESDIS  weekly  snow  cover  analysis,  based  on  AVHRR  imagery,  is 
highly  reliable,  but  it  becomes  obsolete  over  a  few  days,  so  that  NMC  is  exploring  the  use 
of  a  NESDIS  experimental  snow  cover  product,  based  on  a  SSM/I  algorithm  developed  by 
N.  Grody.  Although  there  are  some  areas  where  it  underestimates  snow  cover,  it  generally 
looks  very  promising. 

Tsann  Yu,  from  NMC,  presented  results  of  experimental  assimilations  of  SSM/I 
surface  wind  speeds.  Extensive  verifications  indicated  a  small  but  positive  impact  in  the 
forecasts,  as  well  as  a  reduction  of  the  data  rejected  by  the  Quality  Control  (QC),  indicating 
better  agreement  between  the  6-hour  forecast  and  the  data.  It  was  shown  that  the  use  of 
SSM/I  wind  speed  data  produced  forecasts  which  were  clearly  in  better  agreement  with 
buoy  data,  although  the  analyses  were  slightly  worse.  This  apparent  contradiction  was 
explained  by  the  fact  that  in  the  presence  of  abundant  SSM/I  additional  data,  it  is  much 
more  difficult  for  the  analysis  to  fit  the  relatively  sparse  buoy  data.  It  was  also  pointed  out 
that  SSM/I  data  near  coastal  areas  should  not  be  included  in  the  "superobs"  (1  degree 
averaged)  SSM/I  observations.  As  a  result  of  these  encouraging  results,  NMC  plans  to 
make  the  assimilation  of  SSM/I  wind  speeds  operational  in  the  near  future. 

L.  Phalippou,  from  ECMWF,  presented  comparisons  between  the  SSM/I  total 
precipitable  water  (TPW)  and  the  ECMWF  fields,  as  a  preliminary  evaluation  of  the  current 
ECMWF  analysis  and  its  possible  use  for  data  QC.  Serious  discrepancies  were  found, 
especially  west  of  the  continents  over  stratus  regions  and  regions  of  high  precipitation.  The 
differences  could  not  be  explained  by  rain  contamination  upon  the  SSM/I  algorithm,  since  a 
correction  would  increase  the  discrepancy. 

John  Derber,  from  NMC,  presented  results  of  an  experimental  assimilation  of 
SSM/I  precipitable  water  into  the  NMC  analysis  system  (based  on  a  variational  approach 
denoted  Spectral  Statistical  Interpolation  or  SSI).  He  showed  that  in  the  tropics,  the 
assimilation  of  the  TPW  resulted  in  an  increase  in  the  moisture  immediately  above  the 
boundary  layer  (where  the  model  is  already  close  to  saturation).  As  a  result,  the  excess 
TPW  was  rained  out  during  the  short  range  forecasts.  This  discrepancy  may  be  due  to 
problems  in  the  model  cumulus  parameterization,  which  may  result  in  too  dry  a  tropical 
model  atmosphere.  It  could  also  be  due  to  excessively  high  TPW  estimates  from  the  SSM/I 
since  the  retrieved  TPW  is  quite  sensitive  to  the  details  of  the  algorithm.  For  the  moment 
NMC  is  still  experimenting  with  this  data. 

Bob  Atlas,  from  NASA/GLA,  presented  several  experiments  on  the  assimilation  of 
SSM/I  wind  speeds,  as  well  as  Seasat  and  ERS-I  scatterometer  winds,  and  over  the 
horizon  (OTH)  radars.  Although  the  OTH  has  several  promising  features,  it  also  has 
serious  drawbacks.  The  tests  with  SSM/I  winds  showed  much  better  results  when  the 
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stability  of  the  boundary  layer  was  taken  into  account  in  the  vertical  distribution  of 
observational  increments  due  to  SSM/I  data. 

Andreas  Goroch,  from  NOARL,  showed  the  results  of  a  simple  model  to  describe 
radiative  transfer  through  a  layer  of  brine  over  ice.  He  showed  that  with  the  combination  of 
visible,  infrared  and  microwave  measurements  it  should  be  possible  to  distinguish  between 
flooded  ice,  broken  ice,  and  open  water  regions. 

Chris  Grassotti,  from  AES,  Canada,  presented  a  model  to  determine  precipitation 
rates  at  both  short  (hourly)  and  long  (monthly)  time  scales,  based  on  a  multispectral 
statistical  classification  scheme.  The  results  without  using  SSM/I  data  looked  promising 
when  compared  to  VIS/IR  data,  and  SSM/I  data  will  be  incorporated  in  the  near  future. 

Irene  Rubenstein,  from  EOL,  Canada,  discussed  the  effect  of  the  influence  of 
evolving  cloud  hydrometers  upon  the  satellite  microwave  observations.  Results  using 
observations  from  the  Labrador  Ice  Margin  Experiment  suggested  agreement  between  the 
model  and  observations. 
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SESSION  4  -  APPLICATIONS  TO  NUMERICAL  WEATHER  PREDICTION 


THE  ASSIMILATION  OF  DMSP  RETRIEVAL  PRODUCTS  INTO 
THE  NAVY’S  ATMOSPHERIC  PREDICTION  SYSTEMS 

James  S.  Goerss  and  Patricia  A.  Phoebus 


Naval  Oceanographic  and  Atmospheric  Research  Laboratory 
Atmospheric  Directorate 
Monterey,  CA  93943 


The  utilization  of  various  DMSP  retrieval  products  by  the  Navy  Operational  Global 
and  Regional  Atmospheric  Prediction  Systems  (NOGAPS  and  NORAPS)  run  at  Fleet 
Numerical  Oceanographic  Center  (FNOC)  is  described  in  this  paper.  DMSP  SSM/T 
temperature  soundings  arc  incorporated  into  both  the  global  and  regional  prediction  systems 
as  layer  thicknesses  by  a  multivariate  optimum  interpolation  (MVOI)  analysis.  The 
operational  assimilation  of  SSM/I  wind  speeds  into  the  NOGAPS  MVOI  analysis  was 
initiated  in  September  1990.  Moisture  profiles  are  produced  for  inclusion  in  a  global 
moisture  analysis  using  SSM/I  total  column  precipitable  water  measurements.  The 
NOGAPS  background  or  first-guess  moisture  profiles  at  the  location  of  the  SSM/I 
observations  are  adjusted  so  that  their  total  column  precipitable  water  matches  what  was 
observed. 


161 


DMSP  SATELLITE  DATA  USE  IN 
NAVY  ATMOSPHERIC  PREDICTION  SYSTEMS 


SSM/T  temperature  soundings  from  F10  and  F11 
(NOGAPS  and  NORAPS) 


SSM/I  wind  speeds  from  F10  and  F11  (NOGAPS) 


SSM/I  total  column  precipitable  water  from 
F10  and  F11  (to  be  used  by  NOGAPS) 


DMSP  SSM/T  SOUNDINGS 
TYPICAL  6-HOUR  COVERAGE 


analysis  processing 

OE  SATELLITE  SOUNDING  DATA 


•  Mean  layer  temperatures  used  to  compute  layer  thicknesses,  if  neces¬ 
sary. 

e  Gross-error  checks  performed  on  each  individual  thickness  increment. 
—  Tolerance  T  depends  on  observation  error  and  prediction  error. 

-  Observation  error  for  DSMP  same  as  for  NOAA  clear  retrievals. 

-  Rejection  criteria  depends  on  observation  location. 


Table  6.  Gross-Error  Check  Thickness  Tolerances 


T  varies  with  level  and  data  type  j 

|  Data  Increment  Type 

|  No  Checking 

\  Further  Checking 

|  Rejection  | 

Satellite  Thickness 
Polar  Areas 
Below  50  mb 
50  to  20  mb 
20  to  10  mb 
Extra-tropics 
Below  50  mb 
50  to  20  mb 
20  to  10  mb 


<2.4T 

<  ZST 

<  5  ST 

<  ur 

<  2-4T 

<  3.6T 


None 

None 

None 

None 

None 

None 


>  1AT 
>3  .ST 
>5  ST 

>  1  ST 

>  2  AT 

>  3£T 


Tropics 
Below  50  mb 
50  to  20  mb 
20  to  10  mb 


<L2T 

<  1.9T 

<  23T 


None 

None 

None 


> 

> 

> 


1ST 

1ST 

1ST 


•  Lapse-rate  checks  performed  on  sounding. 

-  Compare  to  vertical  temperature  structure  of  the  background  field. 

•  1000-700  mb  thickness  increment  compared  to  500-300  mb 
thickness  increment. 

-  Lowest  sit  layers  of  the  sounding  may  be  rejected  (below  300  mb). 
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GLOBAL  MOISTURE  ANALYSIS 


•  DMSP  SSM/1  Preclpl  table  Water  Measurements  Are  Used  to 
Adjust  NOGAPS  Moisture  Profiles  Over  Oceans 


PRECStTABLE  WATER 


•  Synthetic  Profiles  Supplement  Radiosonde  Data  Over  Land, 
Resulting  In  Improved  Global  Moisture  Analysis 


SYNTHETIC  MOISTURE  PROFILES 


Produced  at  locations  of  DMSP  SSM/I  total  column 
precipitable  water  (PW)  measurements 


NOGAPS  background  values  of  specific  humidity 
are  multiplied  by  the  ratio  of  SSM/I  PW  to  NOGAPS  PW 
(Equivalent  to  solving  variationally  using  weights 
proportional  to  PW  in  each  layer  (O'Brien,  1970)) 


Adjusted  values  used  to  create  adjusted  profiles  of 
dew-point  depression  for  input  to  moisture  analysis 
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MOISTURE  ANALYSIS  INCREMENTS 
850  MB  DEW-POINT  DEPRESSION 


Demonstrates  all-weather  capability  and  global  coverage  of  SSM/I. 
Ships  tend  to  avoid  high  wind  speed  areas. 


WIND  SPEED  CRTEG0RY  IM/5) 


NOGAPS 

SSM/I  DATA  PROCESSING 


INPUT  processes  output 


NOGAPS 

SSM/I  DATA  PROCESSING 


INPUT  PROCESSES  OUTPUT 
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SSM/I  Quality  Control 


•  Perform  preliminary  checks  prior  to  analysis. 

—  Exclude  observations  not  in  rain-free  areas. 

-  Exclude  observations  outside  the  3-25  m/sec  range. 

•  Compare  SSM/I  wind  speeds  to  preliminary  analysis  wind 
speeds. 

—  Reject  if  background  wind  is  light  and  SSM/I  wind 
speed  is  not. 

—  Reject  if  there  is  a  large  difference  between  SSM/I  and 
background. 


Preliminary 

Rejection 

Percentage 

Analysis 

Criteria 

Rejected 

<4m  s'1 

SSM/I  >7ms-‘ 

■i.iumrrua 

~  1.5% 

|  Diff  |>  10  m  s_1 

~  0.1% 

«  Compare  SSM/I  superob  wind  component  increments  to 
surrounding  data. 

-  Reject  if  no  other  supporting  data. 

-  Very  few  SSM/I  obs  are  rejected  by  this  test. 


NOGAPS  Analysis  Differences 
with  SSM/I  Winds  minus  without  SSM/I  Winds 
12Z,  17  September  1991 


Maximum  Difference  =  £  mb  Maximum  Difference  =  10  ktt . 
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NOGAPS  48-lu-  Forecast  Differences 
with  SSM/1  Winds  minus  without  SSM/I  Winds 
valid  00Z,  15  Match  19!)  1 


1U00  mb  HEIGHT 


Afazimvm  DiJJeixncc  =  82  m 


Ship/Buoy  Vector  RMS  Errors 


Observation  -  Background  Differences 

RMS 

After  SSM/I 

N  Hem 

S  Hem 

mjhKKBH 

wSSaM 

•  Sli  p  and  buoy  observations  agree  more  closely  with  the  background 
wind  speeds  a/ter  the  SSM/1  data  have  been  assimilated. 

•  Improvement  is  most  pronounced  in  the  Southern  Hemisphere. 
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IMPACT  OF  SSM/I-BASED  SNOW/ICE  ANALYSES  IN  NMC’S  ETA 

MODEL 

Kenneth  E.  Mitchell  and  Thomas  L.  Black 

NOAA/NWS/National  Meteorological  Center 
Washington,  DC  20233 

Fedor  Messinger 

UCAR  Visiting  Scientist 
NOAA/NWS  National  Meteorological  Center 
Washington,  DC  20233 

Norman  C.  Grody 

NOAA/NESDIS,  Office  of  Research  and  Applications 
Washington,  DC  20233 


At  various  operational  numerical  weather  prediction  (NWP)  centers,  forecast 
models  on  global,  regional,  and  meso-scales  are  incorporating  increasingly  complex 
parameterizations  of  the  land/atmosphere  interface.  These  parameterizations  require 
improved  specifications  of  surface  conditions,  such  as  snow/ice  cover,  vegetation,  albedo, 
and  soil  moisture.  Such  is  the  case  in  the  experimental  mesoscale  forecast  model  known  as 
the  ETA  model  (so  named  for  its  step-mountain  eta  coordinate)  under  development  at  the 
National  Meteorological  Center  (NMC).  The  physical  parameterizations  in  the  ETA  model 
are  described  in  detail  by  Janjic  (1990). 

Following  the  current  practice  in  NMC's  operational  Regional  Analysis  and 
Forecast  System  (RAFS),  which  includes  the  Nested  Grid  Model  (NGM),  experimental 
runs  of  the  ETA  model  to  date  have  used  a  relatively  coarse  resolution  (nominally  2-degree 
grid)  snow/ice  analysis  produced  weekly  on  Mondays  by  the  Synoptic  Analysis  Branch 
(SAB)  of  NESDIS.  This  operational  NESDIS  snow/ice  analysis  is  manually  derived  on  a 
digitizing  table  by  an  analyst  who  inspects  and  interprets  the  most  recent  cloud-free  visible 
and  IR  imagery  from  the  AVHRR  sensor  of  NOAA's  polar  orbiters.  This  interpretation  is 
complicated  in  winter  by  persistent  cloud  cover  at  high  latitudes  and  the  lack  of  visible 
imagery  during  the  arctic  night. 

The  spatial  and  temporal  resolution  of  the  NESDIS  operational  snow  analysis  is 
insufficient  for  the  high  resolution  specification  of  snow/ice  cover  needed  by  NMCs 
regional  forecast  models.  The  study  by  Petersen  and  Hoke  (1989)  illustrated  the 
detrimental  impact  of  a  week-old  snow/ice  analysis  in  the  RAFS.  In  their  case  study,  the 
untimely  snow  'ice  analysis  led  to  poor  low-level  temperature  forecasts  and  poor  delineation 
of  the  rain-snow  boundary  in  an  area  where  the  snow  cover  had  recently  changed. 

As  an  alternative  to  the  weekly  analysis  of  snow/ice  cover,  NESDIS  is  producing 
an  experimental,  real-time,  automated,  high-resolution  (47  km)  snow/ice  analysis  via 
retrieval  algorithms  (Grody,  1991)  applied  to  the  Special  Sensor  Microwave/Imager 
(SSM/I)  instrument  aboard  the  polar  orbiters  of  the  Defense  Meteorological  Satellite 
Program  (DMSP).  The  SSM/I  has  a  conical  scan  with  a  swath  width  of  1400  km.  At  the 
middle  and  high  latitudes  associated  with  snow/ice  cover,  the  areas  scanned  by  consecutive 
orbits  of  the  SSM/I  give  total  hemispheric  coverage  in  less  than  24  hours  (i.e.  daily 
coverage).  In  addition,  the  SSM/1  retrieval  of  snow/ice  is  not  contaminated  by  cloud  cover 
or  precipitation,  nor  does  it  require  sunlight. 

In  this  study,  we  illustrate  the  impact  of  replacing  the  weekly,  2-degree,  operational 
NESDIS  snow/ice  analysis  with  the  latter  "daily"  47-km,  experimental  SSM/I  snow/ice 
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analysis  as  input  to  the  surface  specification  for  experimental  mesoscale  runs  of  NMCs 
ETA  model  on  a  30- km  grid.  We  focus  on  cases  of  recent  significant  changes  in  snow- 
cover  being  missed  in  the  operational  analysis  due  to  the  weekly  update  frequency.  The 
errors  in  snowcover  in  turn  produce  nontrivial  errors  in  surface  heat  fluxes  (owing  to 
snowcover  induced  errors  in  the  surface  albedo  and  skin  temperature),  and  hence  errors  in 
the  low-level  ETA  model  temperature  forecasts.  These  low-level  temperature  errors  in  turn 
lead  to  errors  in  the  inferred  rain/snow  boundary.  We  will  show  that  both  types  of  errors 
can  be  avoided  in  the  ETA  model  by  alternatively  using  the  timely  specification  of  snow/ice 
cover  afforded  by  the  SSM/I  retrieval. 

References: 

Janjic,  Z.I.,  1990:  The  step-mountain  coordinates:  physical  package.  Mon.  Wea.  Rev., 
118:  1427-1443. 

Petersen,  R.A.,  and  J.E.  Hoke,  1989:  The  effect  of  snow  cover  on  the  Regional  Analysis 
and  Forecast  System  (RAFS)  low-level  forecasts.  Wea.  and  Forecasting.  4:  253- 
257. 

Grody,  N.C.,  1991:  Classification  of  Snow  Cover  and  Precipitation  Using  :’ie  Special 
Sensor  Microwave  Imager.  J.  Geophys.  Res.,  7423-7435. 
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IMPACT  OF  SSM/I-BASED  SNOW/ICE 
ANALYSES  IN  THE  NMC  ETA  MODEL 

OR 

COMPARISONS  OF  SSM/I  SNOW/ICE  RETRIEVALS 
WITH  NESDIS  AND  USAF  OPERATIONAL 
SNOW/ICE  ANALYSES 


KEN  MITCHELL  fNMOT 
J.  PEREIRA  (USAF),  M.  PECNICK  (CENTEL), 
C.  BOETTCHER  (SMSRC),  N.  GRODY  (NESDIS) 


NMC/NESDIS/DOD  DMSP  CONFERENCE 
APRIL  14-15 ,  1992 


DMSP  PRODUCTS 

(NWP  PERSPECTIVE) 

*  SSM/T-1  (NESDIS) 

-  Temperature  profiles 

*  SSM/T-2  (NESDIS) 

-  Water  vapor  profiles 

*  SSM/I  (FNOC/NRL) 

-  Ocean  surface  wind  speed 

-  Ocean  column  water  vapor 

-  Snow/ice 

-  Rain  rates 

*  OLS  (AFGWC) 

-  RTNEPH  cloud  analysis 

--  NWP  cloud  forecast  schemes 
--  rain  rate  estimates  (AGROMET) 
--  access  via  NESDIS/SPN 
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Time  (Hours) 


Global  Mean  Water  Balance 

2/10  -’iH 

Prec  Evap  Prec  Evap 

Op  Op  SSMI pw  SSMI  pw 


Verification 

GSM/A 

GSM/N 

GSM/E 

QSM/QL* 


Global  mean  (mm/day) 


NMC  2001  VISION 


1.  A  NATIONAL  MESOSCALE  MODEL 

-  ULTIMATELY  ON  4  KM  GRID 

2.  UNLIKELY  THAT  OBSERVATIONAL  DATA  WILL 
PROVIDE  A  COMPLETE  DESCRIPTION  OF  THE 
INITIAL  CONDITIONS  ON  SUCH  A  FINE  GRID 

3.  SCIENTIFIC  BASIS  ON  WHICH  OUR  HIGH 
RESOLUTION  MESOSCALE  EXPECTATIONS  REST 
IS  THE  CONCEPT  THAT  MOST  MESOSCALE 
CIRCULATIONS  ARE  DRIVEN  BY  INTERACTIONS 
OF  LARGE  SCALE  WITH  LOWER  BOUNDARY 

-  DETAILED  TOPOGRAPHY 

-  SURFACE  HEAT  AND  MOISTURE  FLUXES 

—  LAND  (Soil  wet.,  Vegetation,  Snow,  Albedo,  ZO) 
--  SEA  (SST's,  Sea  Ice) 

-  CLOUD  COVER 


ETA  MODEL  SUMMARY 

PHYSICS:  (JANJIC,  MWR,  1990) 

-  BETTS/MILLER  CONVECTION 

-  GFDL  RADIATION 

-  MELLOR/YAMADA  2.5  PBL 

-  FORCE/RESTORE  SLAB  SFC 

9  Tg  ■  (l-A)SJ,  +LJ  +  Lt  +  H  +  MA*Ep  +  G  ♦  SNW 

9X  1  2  2  3-4  4  5 

1  -  ALBEDO 

2  -  ClOUDS 

3  -  VEGETATION 

4  -  SFC/SOIL  MOISTURE 

5  -  SNOW/ICE  COVER 
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SNOW  ANALYSIS  CHARACTERISTICS 


NESDIS 

$SM/| 

USAF 

Timeliness 

WEEKLY 

DAILY 

DAILY 

Resolution 

190  km 

23  km 

46  km 

Analyst 

YES 

NO 

YES 

Stations 

NO 

NO 

YES 

ice  Flag 

NO 

YES 

YES 

Reliability 

GOOD 

FAIR 

GOOD 

Cloud  Mask 

YES 

NO 

YES 

5 
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5tM  Lt*U.  PRESSURE  IM8I 

1000-500  «  THICKNESS  (OflM)  R8-H  ETH  FCST 

VflLIO  I2Z  9  OEC  91  80KM  E-GRID 


NESDIS/GRODY  SNOW  RETRIEVAL 
CONSIDERED  IMPROVEMENTS 


1  -  SPACECRAFT  SPECIFIC  COEFFICIENTS 

2  -  MAP  SPATIALLY  AVERAGED  FOOTPRINTS 

3  -  THRESHOLD  TUNING 

-  MERGED  F-10  AND  F-11  MASTER  MAP 

5  -  RESOLVE  ASCENDING  VS  DESCENDING 
PRODUCT  DIFFERENCES 
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CONCLUSIONS 


1)  WEEKLY  NESDIS  OPERATIONAL  SNOW/ICE  ANAL 
IS  TOO  CRUDE  FOR  ETA  MESOSCALE  FCSTS 

2)  NESDIS  SSM/I  SNOW/ICE  RETRIEVAL  NOT  RELIABLE 

—  Misses  large  areas  of  snow 

—  False  ice  ribbons  along  coasts 

—  Lacks  retrieval  over  coastal  points 

3)  USAF  SNOW/ICE  ANALYSIS  BY  ITSELF  WILL  BE 
SOURCE  OF  FIRST  PARALLEL  TESTS 

4)  FOR  THE  FUTURE,  A  BLEND  OF  USAF  AND  SSM/I 
ANALYSES  MAY  BE  SUPERIOR 

—  Requires  improved  retrieval  algorithm 

5)  URGE  AFGWC  TO  USE  SSM/I  SNOW/ICE  PRODUCT 
AS  A  SUBJECTIVE  INPUT  TO  USAF  SNOW/ICE  ANAL 
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USE  OF  SSM/I  WIND  SPEED  DATA  IN  THE  OPERATIONAL 
NUMERICAL  WEATHER  PREDICTION  SYSTEM  AT  NMC 


T.W.  Yu,  W.H.  Gemmill,  and  J.  Woollen 

Development  Division 
National  Meteorological  Center 
Washington,  DC  20233 


The  Special  Sensor  Microwave  Imager  (SSM/I)  system  on  board  the  U.S.  DMSP 
satellite  is  capable  of  measuring  global  ocean  surface  wind  speeds,  among  other  things,  and 
these  wind  speed  data  are  now  routinely  available  at  the  National  Meteorological  Center 
(NMC).  However,  the  SSM/I  wind  speed  data  do  not  include  wind  directions,  and 
therefore  additional  provisions  have  to  be  included  in  the  NMC's  operational  analysis 
system  to  handle  the  wind  data.  A  procedure  for  assigning  wind  directions  to  the  SSM/I 
wind  speed  data  for  use  in  the  operational  global  data  assimilation  system  has  been 
designed  and  tested.  Several  data  quality  control  procedures  specially  designed  for  the 
operational  use  of  the  SSM/I  wind  speed  data  have  been  developed,  which  arc  being 
discussed  in  a  separate  paper.  Results  of  analyses  and  forecasts  based  on  the  use  of  the 
SSM/I  wind  data  for  an  extended  period  in  the  NMCs  operational  assimilation  and  forecast 
system  will  be  discussed.  These  results  are  compared  with  those  of  the  analyses  and 
forecasts  which  do  not  include  the  SSM/I  wind  data  to  assess  the  impact  of  the  data  on  the 
operational  numerical  weather  prediction. 
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COMPARISON  BETWEEN  SSM/I  AND  ECMWF  TOTAL  PRECIPITABLE 

WATER 


L.  Phalippou 


ECMWF 
Shinfield  Park 

Reading  RG2  9AX  United  Kingdom 


The  atmospheric  parameters  which  can  be  retrieved  from  data  from  the  Special 
Sensor  Microwave/Imager  (SSM/I)  on  the  DMSP  satellites  include  the  total  precipitable 
water  (TPW),  cloud  liquid  water  (CLW)  and  surface  marine  wind  speed.  These  three 
parameters  are  of  importance  for  numerical  weather  prediction  (NWP).  This  paper  presents 
the  first  results  of  a  comparison  between  the  SSM/I-TPW  and  the  ECMWF  model  humidity 
field.  SSM/I  data  have  been  acquired  in  the  WetNet/McIDAS  context  for  the  July  to 
September  1987  period.  The  goal  of  this  comparison  is  to  identify  possible  weaknesses  of 
the  ECMWF  humidity  analysis  and  to  control  the  quality  of  SSM/1  products. 
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WATER  VAPOUR  FROM  SSM/I  DATA 


COMPARISON  WITH  ECMWF  ANALYZED  HUMIDITY  FIELDS 


Laurent  PHALIPPOU 

ECMWF 

Shiniield  Parle 

READING 

ENGLAND 

Tel  44  734  499  656 

Fax  44  734  869  450 

Email  LaurentPhalippou@ccmwf.co.uk 


OPERATIONAL  DATA  ASSIMILATION  -  FORECAST  CYCLE 
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SSM/I  at  ECMWF 


THE  WETNET  PROGRAM 

5  years  project  run  by  NASA/MSFC 

Based  on  SSM/I  data  +  McIDAS  software 

Goal:  cooperative  research 
+  EOSDIS  preparation 

WETNET  DATA 

SSM/I  Brightness  temperatures 

SSM/I  Products 

GOES  Imagery  (Infra-red) 

ECMWF  receives  WetNet  data  since  April  91 
Available  now:  July  87  — >  November  87 
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Total  Precipitable  Water  from  SSM/I 


Background 

TPW  =Jgqdp 

g:  gravity  constant 

q:  water  vapour  mixing  ratio 

TPW  Units:  kg/ms  =  0.1  g/cm1  =  1.  mm 

WetNet  TPW 

AlishOLse  et  al.  algorithm  (NESDIS) 

Regression  radiosonde  TPW  vs  SSM/I  Tb 
(for  19  small  islands) 

Period  July  87  to  April  88 

Uses  3  channels: 

TPW=  232.  -  0.14  Tbl9v  -  1.8  Tb22v 

-  0.36  Tb37v  +  0.006  (Tb22v)**2 


QUALITY  CONTOOt  OF  SSM/1-TPW  AT  ECMWF 
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■RtrioJ  ZSStpt  to  IS  6<£  +•  63i)wt«  ZSW.v  U«t=  4-S  <ivj» 
model  -  SSM/I  at  SSM/I  location  (closa  to  a  radiosonde) 


Lat 

Sampf 

Mean 

Ana 

Mean 

SSM/I 

S.D. 

Ana 

S.D. 

SSM/I 

Bias 

90. 

60. 

1891 

10.3 

10.8 

4.2 

3.8 

-0.5 

60. 

55. 

847 

11.0 

11.1 

4.9 

3.8 

-0.1 

55. 

25. 

6226 

25.3 

25.2 

13.0 

12.9 

0.1 

25. 

20. 

914 

44.8 

44.0 

9.8 

10.1 

0.8 

20. 

0. 

3604 

47.8 

48.2 

8.6 

7.5 

-0.4 

0. 

-20. 

353 

50.2 

48.9 

6.4 

7.3 

1.3 

-20. 

-  25. 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

-25. 

-55. 

252 

20.9 

18.8 

7.4 

7.8 

2.1 

-55. 

-60. 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

-60. 

-90. 

26 

6.9 

11.2 

2.2 

3.3 

-4.3 

Total 

90.  -90. 

14113 

30.0 

30.0 

17.2 

17.0 

0.0 

-  ^Soo  r«.obs 
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COMPARISON  BETWEEN  ECMWF  ANALYSEDHUM1DITY 
ANb  SSM/I-TPW  FIELDS 


ANALYSEDHUMIDITY  FIELDS 

Period:  July  id  December  87 
T106  Analysis  ( resotutton-125  km  ) 

4  analysed  Helds  per  day  (  00,  06,  12,  18  UTC ) 
Average  over  14  days  (  56  Helds  ) 


SSM/1  TPW  FIELDS 

AvaiaWe  over  sea  only 
65  km  resolution 

Average  over  14  days  (  ascend.  ♦  descend,  orbits  ) 
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COMPARISON  BETWEEN  ECMWF  HUMIDITY 
AND  SMMR-TPW 

Scanning  Multichannel  Microwave  Radiometer 
Aboard  Nimbus  7  (1978)  and  Seasat 
Frequencies:  6.,  10,  18.,  21.,  37.  GHz  (VH) 
Products:  as  for  SSM/I  in  particular  TPW 

ECMWF  vs.  SMMR-TPW 

Study  done  by  Eymard  et  al.  (JAOT,  1989) 
Comparison  for  20  days  average  (  Jan.-Feb  79) 
Analysis  run  for  FGGE  purposes 
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IMPACT  OF  CLOUDS  ON 
SSM/I-TPW  RETRIEVAL 

Sensitivity  of  TPW  to  clouds 

Cloud  described  by: 

Liquid  water  density  (d  in  g/m’) 

Base  end  Top  of  the  cloud  (Z^jz,  in  in) 

->.  CLW  =  d»(Z,-ZJ 

Sensitivity: 

S=  6TPW/  8CLW=  2  (8TPW/  6TJ*(6TJ  6  CLW) 
6T J  8CLW  is  evaluated  with  a  RT  model 


Results 

S  <  0  ->  TPW  is  under-estimated  in  presence  of 
clouds 

ECMWF-SSM/I<0  cannot  be  explained  by 
clouds  (at  least  not  those  used  in  this  study) 

Stratus  clouds  do  not  explained  large  positive 
differences 
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Table  3:  Sensitivity  of  SSM/I  tuul  predpilablo  water 
10  clou<l  liquid  water 


Cloud  Layer 
mOir 

Liquid  Density 
g/tn’ 

Mld-Lat  ScnslL 

Subarctic  Scnait. 

Tropical  Scusii. 

950-900 

0.1 

■ 

-3.77- 

-5.07 

900-850 

0.1 

I 

-3.75 

-5.05 

850-800 

0.1 

-7.77 

-3.57 

-5.73 

700  650 

0  1 

-8.69 

-2.68 

-7.25 

950900 

1. 

-4.43 

- 

-3.85 

900  850 

1. 

-4.60 

- 

-3.97 

850  800 

1. 

-4.32 

- 

-4.10 

700650 

-3.06 

-4.64 

n*c  sensitivity  la  unitlcss  (l.c.  kg.m  *Ag  nt  * ) 

Mow  to  read  litis  table: 

for  a  cloud  layer  of  thickness  950-900  mb.  l.c.  about  500  m,  with  a  density  of  .1  g/tn’  lltc  cloud  liquid 
content  (CLW)  Is  : 

CLW-  500  m  *  0.1E-3  kg/m*  -  0.05  kg/m* 

Tlie  error  In  1TW  due  to  cloud  Is  Uten  for  example: 

S*CLW«-6.88*0.05“0.344  kg/m* 


SUMMARY 


Good  agreement  between  SSM/I-TPW  and  ECMWF 
analyses  near  radiosondes. 


For  some  areas,  large  differences  between  ECMWF  and 
SSM/I-TPW  have  been  found  and  they  cannot  be 
explained  by  clouds  only.  Further  investigations  are  in 
progress. 


It  is  expected  that  SSM/I-TPW  assimilation  experiment 
will  start  soon  at  ECMWF. 
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METEOROLOGICAL  IMPACT  OF  SURFACE  WIND  DIRECTIONS 
MEASURED  BY  OVER-THE-HORIZON  RADAR  (OTHR)  AND  WIND 
SPEEDS  MEASURED  BY  THE  SPACEBORNE  SSM/I  MICROWAVE 

RADIOMETER 


R.  Atlas,  J.A.  Weinman 

NASA  Goddard  Space  Flight  Center 
Greenbelt,  MD  20771 

T.M.  Georges 

NOAA/Environmental  Research  Laboratory 
325  Broadway 
Boulder,  CO  80303 


The  lack  of  adequate  observational  data  has  been  recognized  as  a  major  factor 
limiting  the  validity  of  weather  forecasts.  Numerical  Weather  Prediction,  NWP,  models 
require  the  complete  specification  of  initial  state  variables,  yet  only  a  fraction  of  such  data 
are  available  at  any  given  time.  Remotely  sensed  measurements  offer  an  effective  way  to 
supplement  those  from  the  conventional  synoptic  network.  The  present  tests  were 
conducted  to  demonstrate  the  impact  on  NWP  of  sea  surface  wind  speeds  obtained  from  the 
SSM/I  spacebome  microwave  sensor  on  DMSP  and  the  sea  surface  wind  directions 
obtained  from  Air  Force  Over  the  Horizon  Radars,  OTH-B.  Georges  et  al.  (1990)  showed 
that  those  DoD  surveillance  OTHRs  provide  a  unique  opportunity  to  obtain  meteorological 
data  from  distances  within  2000  km  of  the  east  and  west  coasts  of  the  United  States. 

Quantitative  surface  winds  were  measured  from  the  west  coast  OTH-B  during  one 
week  in  May  1991  and  from  the  east  coast  OTH-B  during  a  week  in  June  199’  Those 
data  were  incorporated  into  the  Goddard  Laboratory  for  Atmospheres  NWP  model  along 
with  wind  speeds  measured  by  SSM/I.  Results  of  this  analysis  will  be  presented  and  the 
impact  of  the  sea  surface  wind  velocity  measurements  will  be  discussed. 

Georges,  T.M.,  G.D.  Thome,  1990:  An  opportunity  for  Long-Distance  Oceanographic 
and  Meteorological  Monitoring  Using  Over-The-Horizon  Defense  Radars,  Bull. 
Am.  Met.  Soc.,71.  1739-1745. 
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TEOROLCGICAL  IMPACT  OF  SURFACE  WIND  DIRECTIONS  MEASURED  BY 
CVER-THE-HORIZON  RADAR  ( OTHR)  AND  WIND  SPEEDS  MEASURED  BY 
THE  SPACEBORNE  SSM/I  MICROWAVE  RADIOMETER. 


R. 

ATLAS  (1) 

T.M. 

GEORGES  (2) 

J  .  A . 

WEINMAN  (1) 

(1)  NASA  Goddard  Space  Flight  Center 
Greenbelt ,  MD  20771 


(2)  NOAA/Environmental  Research  Laboratory 
325  Broadway 
Boulder,  CO  80303 


NEED  FOR  REMOTELY  SENSED  SURFACE  WIND  DATA 


SURFACE  WIND  VELOCITY  MEASUREMENTS  ARE  REQUIRED  TO: 

1 .  DRIVE  OCEAN  MODELS  AND  SURFACE  WAVE  MODELS. 

2.  CALCULATE  SURFACE  HEAT.  MOISTURE  AND  MOMENTUM  FLUXES 

3.  CONSTRUCT  CLIMATOLOGIES  OF  SURFACE  FIELDS. 

4.  STUDY  A  WIDE  RANGE  OF  ATMOSPHERIC  AND  OCEANIC  PHENOMENA 
-SUCH  AS  EXPLOSIVE  CYCLOGENESIS,  FORMATION  OF  SST  ANOMALIES 
AND  OTHER  ASPECTS  OF  AIR-SEA  INTERACTION. 

5  IMPROVE  NUMERICAL  WEATHER  PREDICTION  -  BY  PROVtCXNG  BOTH 
INITIAL  DATA  AND  VERIFICATION  DATA  FOR  MODEL  PREDICTIONS 
AND  PARAMETERIZATIONS. 

CONVENTIONAL  DATA  OVER  THE  OCEANS  IS  NOT  ADEQUATE  FOR  THESE 
APPLICATIONS  -OATA  IS  USUALLY  SPARSE  AND  NOT  OF  SUFFICIENT 
ACCURACY 
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Table  1-1.  Conventional  Ocean  Surface  Wind  Stress 


Type 

Time 

Period 

Spatial 

Coverage 

Temporal 

Sampl ing 

Motes 

Ship 

reports 

18S0  to 
present 

G1  obal ,  but 
sporadic 

Highly  variable 
with  location 
and  tine 

Typically  of 
poor  accuracy 

Ocean 

Weather 

Ships 

1950  to 

late 

1970's 

Fewer  than  a 
dozen  ships, 
mid-latitudes 

Several  times 
a  day  at  fixed 
stations 

Discontinued 

Moored 

platforms 

I960  to 
present 

Predominantly 
near  coasts 

Several  times 
a  day  at  fixed 
stations 

Inferred 

from 

surface 

pressure 

maps 

1960  to 
present 

Extra-tropical 

Typicilly  ?-« 
per  day 

Requires 
assumed  mode) 
relating 
pressure  to 
wind 

Islands 

1950's  to 
present 

Genera l ly 
poor 

Several  times 
per  day 

May  be 
affected  by 
micro- 

meteorological 

factors 

OTHER  POSSIBILITIES  FOR  SURFACE  WINDS 


NEW  SCATTEROMETER  LAUNCHES 
ERS-1  (1991) 

NSCAT  (1996) 

THESE  SATELLITES  WILL  YIELD  GREATLY  IMPROVED  WINDS; 
MORE  ACCURATE  AND  WITH  FEWER  AMBIGUITIES. 


SEVERAL  SOURCES  EXIST  NOW  FOR  SURFACE  WIND  SPEEDS 
GEOSAT  (ALTIMETER) 

NIMBUS-7  (SMMR) 

DMSP  (SSMI) 

HOWEVER  DIRECTIONS  MUST  BE  INTERRED  FOR  THESE 
WIND  SPEED  DATA  IN  ORDER  TO  USE  THEM  IN  CURRENT 
ANALYSIS/FORECAST  SYSTEMS 
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REAL  DATA  RESULTS 


METHODS  TO  ASSIGN  DIRECTIONS 
TO  WIND  SPEED  DATA 

1.  Use  first  guess'  wind  directions  from  a  GCM 

>  Wind  component  Interpolated  to  obs  locations; 
resulting  directions  assigned  to  speed  data. 

>  Places  complete  reliance  on  qualify  ol  GCM  (6  H) 
forecast. 

2.  Use  analyzed  wind  directions  (STD  ANALYSIS ) 

>  Analyze  all  conventional  data  In  a  first  step; 
then  Interpolate  analyzed  winds  to  observation 
locations. 

>  less  of  a  reliance  on  the  quality  of  the  GCM 
first  guess. 

>  Similar  in  nature  to  the  GLA  objeclive  dealiasing 
o(  SASS  winds. 


3.  Use  analyzed  surface  pressures  (EKMAN  BALANCE) 

>  Very  similar  to  the  approach  described  In  Yu 
(1987). 

>  This  method  involves  estimating  the  surface  drag 
from  the  wind  speeds,  first  guess  pressure 
gradients,  and  an  assumed  balance  (Ekman 
Balance)  relation.  Wind  vectors  are  then 
generated  from  these  quantifies.  There  Is  no  Input 
from  the  model  first  guess  winds. 

4  Use  analyzed  surface  pressures  and  winds  (HYBRID) 

>  Combines  previous  two  methods;  uses  model 
winds  everywhere  but  in  those  areas  (e.g.  South 
Polar  regions)  where  the  Ekman  Balance 
assumption  yields  better  results  than  the  model; 
a  result  of  very  sparse  conventional  data. 
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5.  Use  analyzed  surtoce  pressures  and  winds  (GEN'!.  BAl) 


>  A  generalization  ot  the  Ekman  Balance; 
more  terms  from  the  momentum  equation  are 
used  to  estimate  the  directions. 

>  This  approach  should  extend  the  range  of 
validity  (compared  to  Ekman  approach)  in 
relation  to  the  model  first  guess. 

6  Use  surface  pressures  and  winds  (VARIATIONAL) 

>  Basically  like  std  analysis  method;  a  high 
resolution  variational  analysis  combining 
first  guess,  conventional  data  and  surface 
wind  speed  data  is  performed.  Direction  of 
analyzed  surface  vector  wind  assigned  to 
the  speed  data 

>  Variational  analysis  step  can  contain 
weights  relating  quality  of  the  first  guess 

to  the  other  data  sources.  Explicit  dynamical 
balances,  as  well  as  constraints  on  smoothness 
can  be  put  Info  the  functional. 


GLOBAL  DIRECTIONAL  ERROR  (IN  DEGREES!  AFTER  ANALYSIS  OF  SIMULATED 

DATA 

_ 0-5m/s  5-10  m/s  10-15  m/s  15-20  m/s  ALL  SPEEDS 


METHOD 


1 

41.8 

18.8 

18.0 

24.9 

27.2 

2 

35.0 

18.1 

17.3 

23.1 

23.3 

3 

62.3 

30.8 

17.4 

9.0 

40.1 

4 

35.0 

14.9 

12.2 

8.3 

21.8 

6 

33.8 

15.0 

14.4 

20.4 

2119 

number:  532 

681 

218 

15 

1446 

Method  1:  Model  first  guess 
Method  2:  Analyzed  wind  directions 
Method  3:  Ekman  balance 

Method  4-  Hybrid  (combination  ot  methods  2  and  3) 
Method  5:  General  balance 
Method  6;  Variational  analysis 
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GENERAL  BALANCE  DIRECTION  ASSIGNMENT  METHOD 


•  Basis  for  method:  SPEED  EQUATION 


ds  v;  „„  dS  1  V  „  c2  ( 1 

—  +  V  •  VS  ♦  o» -  +- - ■  VP  +  Cn  S  -  0 

dt  dp  p  s  D 


•  If  V  -  S  COS'f  i  +  5  S/JVy  /  then  the  direction  equation  is: 


c  Bracket  terms  multiplying  the  SIN, COS  terms  are 
affected  by  Pressure  data  and  the  SSM/I  speed  data 

>  A  simple  scalar  analysis  of  SSM/I  speed  performed, 
then  gradients  are  calculated 

>  A  diagnostic  relation,  also  using  the  SSM/I  speed  data, 
is  used  to  estimate  the  friction  term  C0S2 

>  All  other  terms  are  estimated  from  fhe  GCM  first  guess 
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COMPARISON  OF  METHODS  1  AND  5  USING  SIMULATED  DATA 


NATURE  RUN 
(PERFECT)  SSM/1 
SURFACE  WINDS 


GENERAL  BALANCE 
SSM/I  SURFACE  WIND 
VECTORS 


METHOD  1 

SSM/I  SURFACE  WIND 
VECTORS 


Fig.  13 


211 


VARIATIONAL  DIRECTION  ASSIGNMENT 
FOR  SSM/I  WIND  SPEEDS 


PROBLEM: 

GIVEN: 

METHOD: 


CHECKING: 


Want  surface  vector  winds  from  SSM/I  data 

(1)  Sparse  Conventional  data  (ships,  buoys) 

(2)  First  guess  (ECMWF  analysis) 

(3)  SSM/I  surface  wind  speeds 

(4)  Kinematic,  dynamic  constraints 
(smoothness,  divergence,  vorticity  tendency) 

Obtain  analyzed  wind  Va  by  performing  a  variational 

analysis  using  the  given  information;  assign  the 

resulting  direction  to  the  SSM/I  wind  speeds: 

>  Construct  a  functional  F  consisting  of  terms 
measuring  the  departure  of  a  solution  (wind 
vectors)  from  the  above  terms. 

>  Find  the  solution  which  minimizes  F  using  a 
direct  minimization  algorithm.  Follows  work 
at  GLA  by  R.  Hoffman  on  dealiasing  SEASAT 
winds. 

>  Interpolate  the  directions  of  V„  to  the  SSM/I 
locations  -  these  are  the  assigned  directions  of 
Method  6 


Resulting  wind  voctors  are  checked  against  r 
set  of  withheld  buoys. 
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For  our  analysis,  we  minimize  the  following  functional  F: 

F  =  A,,SC  +  X^S,  +  XjS,  +  X4SveJ  +  +  X6Svo,  + 

Here,  the  ^.-weights  control  the  amount  of  influence  each  constraint 
has  on  the  final  analysis. 

Penalty  terms  used  in  the  functional  measuring  departures: 

.  sc  =  i  <vc  -  vy  from  conventional  data 

.  s,  =  e  (V,  -  VJ2  from  first  guess 

•  S,  =  E  ( I I  -  I  Va  l  )2  from  SSM/I  wind  speeds 

Three  smoothness  constraints  (vs  the  first  guess) 

.  Sve,  =  E  (V2(u,  -  UQ))2  +  E  (V2(v,  -  va))2 

•  Sd(V  =  E  (V2(X,  - 

.  SVOf  =  E  (V2(')/|  -  v,/a))2 

Finally,  a  dynamical  constraint  on  the  analyzed  vorticitv  tendency: 

.  sdyo  =  E  (  )2 
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QUALITY  CONTROL  USED  IN  VARIATIONAL  ANALYSIS 


Three  checks  are  made  to  check  a  wind  observation 

voe$  =  <Uo«  ,  v^)  against  a  reference  wind 

VMf  =  (u,*,:  ,  vg(f),  wn(~h  in  this  case  is  taken  from  the  output  of 

the  first  minimization. 

1.  Vector  velocity  differences 


((Uobs  -  uBEf)’  +  <Voes  -  vRW)i)'/1  >  y,  IVAVI 

where 

VAV  =  (V^  +  VKf)/2 

2.  Speed  differences 

I  OV^I  -  IVKfl)  l  >  Yj  IVAVI 

3.  Direction  differences 

cos(  0o&s  -  0„l(  )  >  Yj  cos(  0CU1  ) 


here  we  use  0CUI  =  60  degrees. 


F  is  minimized  by  means  o f  two  direct  minimization  steps: 

I .  First  step  uses  gB  conventional  data,  using  a  loose 
convergence  criterion.  (The  conventional  data  were 
subject  to  a  mild  quality  control  In  a  pre-processing 
step) 

Z  Repeat  the  minimization  process,  but  now  perform  gross 
checks  on  the  conventional  data  (vs  the  first  minimization 
field).  Use  a  tighter  convergence  criterion,  since  the  first 
step  should  get  a  solution  near  the  global  minimum  of  F. 

The  cunent  version  of  the  variational  analysis  uses  the  same  X's  for 
both  steps;  the  greatest  weights- are  given  fo  the  conventional 
data  and  the  SSM/I  speeds,  the  least  weights  are  given  to  the 
smoothness  constraints  and  the  first  guess. 


COMPARISON  OF  DIRECTIONAL  ASSIGNMENT  VS.  OEAUASING  OF 
REAL  SEASAT  SCATTEROMETEn  DATA  USING  METHOD  6 


JULY-OCTOSER  COLLOCATION  STATISTICS  USING  INDEPENDENT  DATA 


MEAN  ABSOLUTE  ERROR 

BIAS 

WIND  SPEED:  1.7  m/S 

0.24  m/s 

DIRECTION:  21.7  deg 

0.10  deg 

NUMBER  OF  COLLOCATIONS:  3870 

SSM/I  (METHOD  61  VS  PMEL  BUOYS 

BIAS 

WIND  SPEED:  1.7  m/s 

*1.6  m/s 

DIRECTION:  27.5  deg 

•6.1  deg 

NUMBER  OF  COLLOCATIONS:  124 


IMPACT  OP  KM  OAT  A  OH  GCA  ANALYSIS  PON  «  AUG.  1MT 


MO  m  croro't*-' IN.  •<  -  »*<  fm>«  in-jrjC* 

i  •»  w  -  i  x  o  - 
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Contioi  and  Control  «  S’.M/i  <uo  to  8  j()  <t\o> 


•iOC  •*  «OPOfC«*HN  *A  IO*»  •  V*.  c«l*«  TAOAICS 
*•»  ms  »os  i  tx  0  nv 


fiouip  ?t>  *L?'T«soT'0'e  SK  ,r,i 

ne*gM  <.y  omtXv  c  a"  *: *a  1  *x>i  to*  *  cxa  c  ou'i 
Conitcn  iC)  C  •  SSM/i  <st^cx:eonv>  C  *  SSM'i 
(«X)  »o  8f£  n-*)i  and  <  •  cloud  nocked  wi nos 
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ECMUF  vs  SSMI  *4  50AY  JUL87  WIND  VEC.M/S 


•>fc 


COMPONENTS  OF  THE  GLA  GCM  PBL  SCHEME 
(Helfand  and  Labraga) 


•  A  level  2.5,  second-order  turbulence  closure 
model  (Mellor  and  Yamada,  1974;  Yamada,  1977; 
Helfand  and  Labraga,  1988) 

•  A  Monin-Obukhov  similarity  scheme  for  the 
"Extended  Surface  Layer" 

•  A  parameterization  for  the  viscous  sublayer 

•  Parameterizations  to  predict  the  surface 
roughness  parameter  z0  over  land  and  sea 
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o-TSi-i  w,M6i  UM  ypsi 

-*A-  5  Z*XArC’L-.iA,’.>’  tot*  tClODO  wi*?S 


LFLX  (55M1-EC1000)  (PA5C)  AUGUST  1987 


••Ml'C.IMtlCGI  •  0  «tOO*t  -01  COMTOUI  IMT(*Y«L>  t.00 


Impact  of  SSMI  Surface  Winds  on  Air-Sea  U-Wind  Stress  (N/m2) 

August  1987 


a  <o  -a  n  -o  « »  -a  n  -a  ej  -a.u  too  102  1.01  •  «  0. or  •  at  >•. it 


Impact  of  SSMI  Surface  Winds  on  Air-Sea  V-Wind  Stress  (N/m2) 

August  1987 
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How  does  the  ionosphere  distort  spectrum  of  the  sea  echo? 
Azwrwth  •  1S09deg  157.1  deg  160.9 deg 


Tins  range  Doppler  pioi  shows  how  the  Bragg  hoes  ol  the  sea  echo  are 
shifted  and  spread  in  frequency  by  reflection  from  an  irregular  ionosphere 
Tick  marks  al  the  bottom  indicate  the  theoretical  Doppler  shills  of  ttie  sea 
echo  Bragg  lines,  and  the  Doppler  shift  scale  goes  from  *1.25  Hz  lo  1  2b 
Hz.  The  relative  strengths  of  the  advancing  and  receding  Bragg  lines  in  the 
three  beams  reveal  the  circulation  in  hurricane  Claudette 


How  does  OTH  Radar  detect  sea  surface  conditions7 


Ionosphere 


^ocean  =  y  “  7-5  m 

Doppler  Shill  fo  =  ±0.1  =  iO-45  Hz 


Sea-Echo  Spectrum: 


Advancing 

Waves 


I  J  I 

♦  1  Hz  I 


kv. 


heceding 
.  Waves 


uv'n — , 

0  fn  -1Hz 


n  ”  ’  n 

Doppler  Shift 


I  lie  limn  ol  ‘r-'fl.H  i>  tUiltiM*  Imni  moving  ommii  w.iv  nor*". illy 

discarded  »n  d'-lmiM*  .ipplt*  .itinriv  <  onlam’.  mfoMii.ilino  ,iU*il  <»«  »mh  mm 
face  condi ti'  >ns 
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RADAR  DIRECTIONAL  ERROR 


DIRECTIONAL  ERROR 


--X--GLAFG  RADAR  -  GLA  -**-  GLA  AN  ■«-  RADAR -GLA 
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CONCLUSIONS 


1.  SSM/I  WIND  SPEEDS  CAN  BE  EFFECTIVELY  ASSIMILATED  USING  DIRECTIONAL  ASSIGNMENT 
METHODS. 

-  PRELIMINARY  RESULTS  SHOW  A  GENERALLY  POSITIVE  IMPACT  OF  SSM/I  WINDS  IN  DATA 
ASSIMILATION,  BUT  POTENTIAL  FOP.  NEGATIVE  IMPACTS  EXISTS,  PARTICULARLY  WITH  SIMPLE 
DIRECTION  ASSIGNMENT  METHODS  '  YHERE  THE  FG  DIRECTIONS  HAVE  LARGE  ERRORS. 

-  CAUTION  MUST  BE  USED  IN  VERTICALLY  EXTENDING  THE  INFLUENCE  OF  SURFACE  WINDS  ABOVE 
THE  PBL.  THIS  EXTENSION  SHOULD  DEPEND  UPON  THE  SYNOPTIC  SITUATION. 

-  LARGE  POTENTIAL  EXISTS  FOR  SSM/I  TO  IMPROVE  AIR-SEA  FLUXES  AS  PART  OF  A  DATA 
ASSIMILATION  CYCLE;  THIS  MAY  HOLD  THE  GREATEST  POTENTIAL  FOR  IMPROVING  NWP  AS 
WELL. 

2.  OTHR  OBSERVATIONS  CONTAIN  USEFUL  INFORMATION  FOR  A  VARIETY  OF  METEOROLOGICAL 
AND  OCEANOGRAPHIC  APPLICATIONS,  BUT  THE  UTILITY  OF  THESE  DATA  ARE  LIMITED  BY  THE 
DIRECTIONAL  AMBIGUITY  AND  THE  OVERLY  SIMPLISTIC  METHODOLOGY  BY  V/HICH  SURFACE 
WINDS  ARE  INFERRED. 

-  THE  UTILITY  OF  OTHR  WINDS  COULD  BE  SUBSTANTIALLY  ENHANCED  BY  AN  IMPROVED 
ALGORITHM  FOR  SURFACE  WIND  DETERMINATION  AND  BY  THE  APPLICATION  OF  OBJECTIVE 
AMBIGUITY  REMOVAL  METHODS,  NOT  REQUIRING  A  MODEL  FG  (eg.  AI/PATTERN 
RECOGNITION/NEURAL  NET  TECHNIQUES). 
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OBSERVATIONS  OF  FLOODED  ICE  IN  ARCTIC  REGIONS 


Andreas  K.  Goroch  and  Robert  W.  Fett 


Naval  Oceanographic  and  Atmospheric  Research  Laboratory 
Atmospheric  Directorate 
Monterey,  CA  93943-5001 


The  remote  detection  of  flooding  of  the  Arctic  ice  pack  has  been  difficult  to 
accomplish  using  conventional  methodologies.  We  have  examined  a  case  of  flooded  Arctic 
ice  with  NOAA  and  DMSP  visible  and  infrared  measurements,  supplemented  by  DMSP 
microwave  imager  (SSM/I)  data.  Analysis  of  visible  and  infrared  data  for  a  sunglint  region 
was  used  to  show  the  distribution  of  flooding  at  27  N,  135  W  on  8  June,  1989. 

A  simple  model  is  described  for  the  radiative  transfer  at  microwave  wavelengths 
through  a  brine  layer  over  ice.  The  predictions  of  the  model  were  found  to  be  consistent 
with  the  spatial  behavior  of  the  SSM/I  measurements.  This  case  study  demonstrated  the 
utility  of  using  combined  visible,  infrared,  and  microwave  measurements  to  differentiate 
among  flooded  ice,  broken  ice,  and  open  water  regions,  even  in  the  presence  of  cloud 
cover  and  surface  fog. 
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Observations  of  Flooded  Ice  in  Arctic 

Regions 


Andreas  K.  Goroch  and  Robert  W.  Fett 


Purpose 

□  Identify  flooded  regions  over  ice  using  microwave  imagery 


Approach 


□  Develop  physical  model 


□  Identify  flooded  regions  with  optical  imagery  (Vis,  NIR, MIR, Thermal) 


n  Identify  microwave  characteristics 
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Flooded  Ice  In  Arctic  Regions 

Frequency  Dependence  of  Surface  Emission 


qi 


Hooded  Ice  in  Arctic  Regions 


SSM/I  Thover  water,  ice,  and  flooded  ice 

19v  vs  37  v  8rLghtness  LeraperaLure 

200 
270 
250 
2S0 
240 
230 
220 
210 
200 

o  ooaoooooo 
cr*o  —  r*  r*  ■»*  u->  uj 

—  OJ  <N  04  CM  <N  <N  OJ 

I9v  Tb 

Flooded  Ice  in  Arctic  Regions 


SSM/I  Brightness  Temperature  Variations 


CLASSIFICATION-BASED  RAINFALL  ESTIMATES  USING  SATELLITE 
DATA  AND  NUMERICAL  MODEL  OUTPUT 

Christopher  Grassotti  and  Louis  Garand 

Atmospheric  Environmental  Service 
2121  Trans-Canada  Highway 
Dorval,  H9P  1J3 
Canada 


As  part  of  ongoing  efforts  at  AES  to  utilize  satellite  data  within  the  Canadian 
Meteorological  Center  forecast  and  analysis  system  we  have  been  exploring  methods  of 
retrieving  rainfall  at  both  short  (hourly)  and  longer  (monthly)  timescales.  It  is  anticipated 
that  retrieved  rain  rates  at  hourly  timescales  may  be  incorporated  within  a  diabatic 
initialization  procedure  at  the  start  of  each  forecast  run  -  hopefully  reducing  the  "spin-up" 
problem. 

The  retrieval  method  to  be  presented  uses  an  objective  multifeature  classification 
approach  in  which  each  analysis  box  (1.25°  lat/lon)  is  classified  based  on  a  variety  of 
features  extracted  from  both  visible  and  infrared  imagery.  Notably,  we  also  use  predicted 
fields  of  vertical  motion  from  the  operational  forecast  model  as  an  additional  feature.  A 
class-dependent  rain  rate  is  then  assigned  to  the  analysis  box.  Preliminary  results  (GPCP) 
AIP-1  (Japan)  have  indicated  that  the  addition  of  the  model  fields  do  improve  estimates  of 
rainfall  compared  with  those  obtained  using  a  Vli/iR  only  approach.  Currently, 
microwave  observations  are  not  used  in  the  retrievals,  however  we  do  intend  to  include 
SSM/I  measurements  in  application  of  this  retrieval  method  to  the  GPCP  A1P-2  (Great 
Britain)  data  set. 


CLASSIFICATION-BASED  RAINFALL  ESTIMATES 
USING  SATELLITE  DATA  AND  NUMERICAL 
FORECAST  MODEL  OUTPUT 


Christopher  Grassotti  and  Louis  Garand 


Atmospheric  Environment  Service 
Dorval,  Quebec 
CANADA 


15  April  1992 


OHJECTIVE:  Improvement  of  rainfall  estimation  at 

synoptic  and  climate  scales.  Ultimately 
to  be  used  within  a  diabatic 
initialization  procedure  at  start  of  CMC 
operational  forecast  runs. 

APPROACH:  Combine  information  from  both  satellite 
(VIS/IR)  and  numerical  model  forecast 
fields,  (presently  no  MW  data  used.) 


TECHNIQUE:  Multifeature  classification  using 

minimum  distance  criteria  in  feature 
space.  Mean  rain  rate  is  assigned  to  each 
class. 


DATA  SET:  The  GPCP  first  AIP  observing  period 
(June,  July/Aug.  ’89). 


1.  CLASSIFICATION 


•Initial  Clustering:  establish  most  useful 
features/classes.  Use  the  minimum  squared 
error  criterion  (Duda  and  Hart.  1973). 

h=t J. 

1=1 


Jt=  Z  I*-™-)2 

xe  X, 


Move  sample  from  class  i  to  j  if  Jj  decreases 
more  than  Jj  increases. 


•Classification:  place  sample  in  class  for  which 
metric  distance  to  class  mean  is  minimized. 


2.  DATA  SETS 

June-August  1989,  GPCP  AIP  1  over  Japan.  Data 
available  hourly  at  1.25°  Iat/Ion. 


(1)  AMeDAS  Network:  Rain  gauge/Radar  composites 


(2)  VIS/IR  imagery  from  GMS.  Analyzed  using 

technique  of  Garand  (1988)  to  obtain  AL,TT,CF 
at  1.25°. 


(3)  Model  forecast  fields  of  rain  rate  interpolated 
to  same  grid.  Model  is  CMC  global  spectral 
(T106,  Kuo  convection). Sequence  of  model  6  hr 
forecasts,  initialized  every  6  hours  were 
interpolated  to  intermediate  hours  permitting 
estimation  at  every  hour. 
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3.  RETR1KVALS  (June) 

•Daytime:  3  Peat.,  30  Classes  (TT.AL.RRf). 
•Nighttime:  3  Feat.,  30  Classes  (TT,CF,RRf). 
•Decision  tree: 


(1)A11  CF  <  96  %,  RR=0. 


(2)IR  temp  screening: 

TT<210  K,  RR=5.7  mm/h 
TT<220  K,  RR=4.3  mm/h 


(3)Classification 

Classify  1.25°x!.25°  box. 
Assign  mean  RR  of  class. 


•Application: 

Dependent  set=first  and  last  7.5  days  of  June. 
Semi-independent  scl-cn  tire  month. 


RAIN  RATE  RETRIEVAL  ERRORS 


Monthly 

(mm) 

Hourly  (mm/h) 

M  o  n  t  h 

Method 

R 

RMS 

BIAS 

R 

RMS 

BIAS 

VIS+IR+Model,  30  cl. 

.846 

53.0 

1.8 

.621 

.604 

-.001 

VIS+IR+Model,  20  cl. 

.839 

54.3 

0. 1 

.61  1 

.61  1 

-.004 

June 

Model 

.806 

56.9 

- 1  5.4 

.541 

.649 

-.027 

VIS+IR 

.740 

67.2 

11.2 

4  69 

.68  1 

.012 

IR 

.628 

73.5 

1  .6 

.4  10 

.703 

-.002 

VIS+IR+Model  (June) 

.730 

95.1 

-46.0 

490 

.615 

-.07  1 

July/ 

VIS+IR+Model 

.680 

84.  1 

-  1  6.6 

.5  75 

.592 

-.028 

Aug 

Model 

.481 

92.7 

-11.1 

188 

.740 

-.0  18 

VIS+IR 

.610 

88.1 

-18.9 

.538 

.59  1 

-.03  1 
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4.  CONCLUSIONS 


•By  combining  model  forecast  estimates  with 
geostationary  satellite  data,  we  improve 
estimates  of  rainfall  at  hourly  and  monthly 
time  scales  over  estimates  obtained  from  either 
satellite  or  model  alone. 

•Method  needs  more  testing  and  validation  to 
resolve  issues  of  dependence  on  regional 
climate  and  rainfall  regime,  e.g.  validity  over 
midlatitude  stratiform  precipitation  unknown. 
(GPCP-2,Gr.  Britain). 

•Approach  can  easily  accomodate  additional 
features/predictors  (SSM/I,  cloud  texture,  etc.) 
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THE  INFLUENCE  OF  CLOUD  MICROPHYSICAL  STRUCTURE  ON 
SPACEBORNE  MULTISPECTRAL  MICROWAVE  OBSERVATIONS 

Irene  G.  Rubenstein 


Earth  Observation  Laboratory 
Institute  for  Space  and  Terrestrial  Science 
4850  Keele  Street,  North  York 
Ontario,  Canada  M3J  3K1 


The  evolution  of  microwave  radiative  properties  of  precipitating  clouds  is  explored 
in  this  study.  Numerical  simulation  of  hydrometeor  effects,  utilizing  realistic  atmospheric 
profiles  and  cloud  microphysical  structure,  on  the  upwelling  microwave  radiation  is 
evaluated.  The  effects  of  ice  hydrometeors  and  melting  layer  presence  have  been 
considered.  The  top  of  the  atmosphere  brightness  temperatures  at  several  frequencies  are 
computed  by  means  of  radiative  transfer  model  for  the  surface  backgrounds  of  ice  free, 
rough  ocean  and  low  sea  ice  concentrations.  The  sensitivity  of  the  calculated  brightness 
temperature  to  uncertainties  of  the  input  environmental  parameters  is  evaluated  for  test 
scenarios  from  Labrador  Ice  Margin  Experiment  (LIMEX)  and  other  documented 
oceanographic  research  campaigns.  Results  of  a  retrieval  procedure  developed  utilizing  the 
simulation  results  appear  to  be  in  qualitative  agreement  with  the  available  surface 
observations. 
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fuli  iict*  in  Accuracy  of  the  retrieval  algorithm  performace  will  increase  once  these  tasks  arc 
completed.  Remote  sounding  of  atmospheric  parameters  from  meteorological  satellites  has 
made  a  significant  impact  on  weather  forecasting  and  our  understanding  of  large-scale  cir¬ 
culation  patterns.  The  derivation  procedures  for  quantitative  estimates  of  vertical  cloud 
composition  and  structure  can  now  be  tested  with  actual  data. 
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ATMOSPHERIC  AND  ENVIRONMENTAL  RESEARCH,  INC. 


Subject;  NMCINESDIS/DoD  Workshop  on  DMSP  Retrieval  Products 

Dear  Colleague: 

A  recurring  theme  at  recent  satellite  meteorology  and  numerical  weather  prediction  symposia  has  been  the  availability 
and  relevance  of  products  derived  from  the  Defense  Meteorological  Satellite  Program  (DMSP)  sensors,  particularly  the  SSM/T 
and  SSM/1  products.  As  real  time  access  to  these  data  has  been  improved  through  the  Shared  Processing  efforts  of  both  the 
Department  of  Defense  (DoD)  and  NOAA's  National  Environmental  Satellite,  Data,  and  Information  Service  (NESDIS),  the 
scientific  community's  interest  in  the  use  of  these  data  has  increased. 

In  response  to  this  interest,  a  Workshop  on  DMSP  Retrieval  Products  has  been  planned  by  scientists  at  the  NOAA/ 
National  Meteorological  Center  (NMC),  NOAA/NESDIS,  and  the  DoD.  The  objective  of  the  workshop  is  to  bring  together 
individuals  from  both  the  defense  and  civilian  communities  with  a  common  goal  of  understanding  the  use  of  DMSP  retrieval 
products.  The  format  will  consist  of  contributed  papers  and  sufficient  time  for  discussions.  It  is  hoped  that  this  interaction 
will  benefit  all  concerned.  The  workshop  will  be  organized  by  Atmospheric  and  Environmental  Research,  Inc.  (AER). 

The  Workshop  on  DMSP  Retrieval  Products  will  be  held  14-15  April  1992  at  the  NOAA  Science  Center  in 
Washington,  DC.  A  tentative  agenda  is  attached.  You  are  invited  to  contribute  a  1 5  minute  presentation  describing  your 
recent  DMSP  related  work.  Please  submit  a  short  abstract  for  inclusion  in  the  proceedings  of  the  meeting  by  15  November. 

If  you  do  not  wish  to  present  a  paper  but  are  interested  in  attending,  please  call  or  write  to  the  workshop  chairman  (617-547- 
6207  or  rgi@aer.com). 

We  look  forward  to  your  active  participation. 


Sincerely, 


Ron  Isaacs 
Workshop  Chairman 


P.  Krishna  Rao,  NOAA/NESDIS 


Robert  McClatchey,  AF  Phillips  Laboratory 
Workshop  Organizing  Committee 


840  MEMORIAL  DRIVE.  CAMBRIDGE,  MA  02139  TEL;  (617)  547-6207 
FAX:  (617)661-6479  FAX:  (617)  876-7158  TELEX:  95-1417 
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Final  Agenda 


NMC/NESDIS/DOD  CONFERENCE  ON  DMSP  RETRIEVAL  PRODUCTS 

14-15  APRIL  1992 

Ron  Isaacs,  AER,  Inc.  Eugenia  Kalnay,  NOAA/NMC 

Workshop  Chairman  P.  Krishna  Rao,  NOAA/NESDIS 

Robert  McClatchey,  AF  PL/GPA 

Workshop  Organizing  Committee 

Tuesday,  14  April 
Morning: 

8:15  Opening  Remarks:  Kalnay/Rao/McClatchey 

8:30  Session  1  -  Overview  of  Current  and  Future  DMSP  Sensors  and  Products 

Session  Chairman:  Ron  Isaacs,  AER,  Inc. 

8:30-8:50  "SSM/I:  Status  Report" 

J.  Hollinger  (Naval  Research  Laboratory) 

8:50-9:10  "Calibration/ Validation  of  the  DMSP  Microwave  Water  Vapor 

Sounder  (SSM/T-2)" 

V.  Falcone,  J.  Morrissey,  M.  Griffin  (Phillips  Laboratory),  D.  Boucher  and 
B.  Thomas  (The  Aerospace  Corporation),  R.  Isaacs  and  J.  Pickle 
(Atmospheric  and  Environmental  Research,  Inc.) 

9:10-9:30  "Performance  of  DMSP  Special  Sensor  Microwave  Humidity 
Profiler  (SSM/T-2):  Preliminary  Results” 

D.  Boucher  and  B.  Thomas  (The  Aerospace  Corporation) 

9:30-9:50  "DMSP  Block  6:  A  Data  Driven  Approach" 

M.W.  Borden,  F.  Kelly,  G.  Mandt  (HQ,  Air  Force  Systems  Division), 

D. L.  Glackin  and  J.  Bohlson  (The  Aerospace  Corporation) 

9:50  Session  2  -  Operational  and  Quasi-Operational  Retrieval  Approaches 
and  Products 

Session  Chairman:  R.  McClatchey,  Air  Force  Phillips  Lab, 

Geophysics  Directorate  (GPAS) 

9:50-10:10  "NOAA/NESDIS  Operational  SSM/T  Products  and  Distribution" 

E.  Burdsall  (NOAA/NESDIS) 

10:10-10:25  Break 

10:25-10:45  "NOAA/NESDIS  SSM/T  Algorithms  and  Evaluation" 

T.  Reale  and  D.  Donahue  (NOAA/NESDIS) 

10:45-11:05  "Defense  Meteorological  Satellite  Program  Microwave 

Radiometer  Processing  at  Air  Force  Global  Weather  Central" 

T.J.  Neu  (HQ  AFGWC) 
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Tuesday,  14  April 
Morning: 


Session  2  -  Operational  and  Quasi-Operational  Retrieval  Approaches 
and  Products  (continued) 

11:05-11:25  "DMSP  Special  Sensor  Microwave/Imager  Processing  at  the 
Fleet  Numerical  Oceanography  Center:  Operational  Status, 
Applications  and  Plans" 

M. C.  Colton  and  C.J.  Cornelius  (Fleet  Numerical  Oceanography  Center) 

11:25-11:45  "An  Overview  of  the  Experimental  SSM/I  Orbit-by-Orbit 
Products  System  (SSMIPROD)  at  NOAA/NESDIS" 

J.  Fiore  (S  M  Systems  and  Research  Corporation) 

Tuesday,  14  April 
Afternoon: 

1:00  Session  2  -  Operational  and  Quasi-Operational  Retrieval  Approaches 
and  Products  (continued) 

1:00-1:20  "SSM/I  Mapped  Experimental  Environmental  Image  Products" 

C.  Boettcher  (S  M  Systems  and  Research  Corporation) 

1:20-1:40  "SSM/I  Total  Precipitable  Water  Vapor  Algorithms:  A  Reprise 
and  Update" 

J.  Alishouse  and  R.  Ferraro  (NOAA/NESDIS) 

1:40-2:00  "The  Use  of  the  DMSP  SSM/I  for  the  Generation  of 

Precipitation  Products  at  NOAA/NESDIS:  Part  I:  A  Status 
Report" 

R.  Ferraro  and  N.  Grody  (NOAA/NESDIS) 

2:00-2:20  "The  Use  of  the  DMSP  SSM/I  for  the  Generation  of  Precipitation 
Products  at  NOAA/NESDIS:  Part  II:  Scientific  Results" 

N.  Grody  and  R.  Ferraro  (NOAA/NESDIS) 

2:20-2:40  "Production  and  Evaluation  of  Experimental  SSM/I  Ice  and 
Wind  Products  at  NOAA/NESDIS" 

W.  Pichel  (NOAA/NESDIS) 

2:40-3:00  "Real  Time  Quality  Control  of  SSM/I  Wind  Speed  Data" 

M.  Waters,  W.S.  Richardson,  W.H.  Gemmill,  and  C.M.  Caruso  (NOAA) 

3:00  Break 

3:15-3:35  "Status  of  NESDIS  Distribution  of  SSM/I  Products" 

P.  Taylor  (NOAA/NESDIS) 

3:35-3:55  "Processing  and  Distributing  SSM/I  Data  via  the  Wet  Net  Project" 

H.M.  Goodman  (NASA  Marshall  Space  Flight  Center) 
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Tuesday,  14  April 


Afternoon: 

3:55  Session  3  -  New  Retrieval  Approaches  and  Products;  Relevance  to  AMSU 

Session  Chairman:  George  Ohring,  NOAAINESDIS 

3:55-4:15  "Joint  NESD1S/NASA/NMC  Effort  to  Develop  an  Advanced 
Satellite  Retrieval  System" 

W.  Baker  (NOAA) 

4:15-4:35  "Classification  Retrieval  Approaches  for  DMSP" 

L.  McMillin  (NOAA/NESDIS) 

4:35-4:55  "Using  SSM/T  Data  as  Proof  of  Concept  for  an  AMSU- A 
Retrieval  Algorithm” 

H.  Fleming  (NOAA/NESDIS)  and  E.  Kratz  (S  M  Systems  &  Research  Corp.) 

4:55-5:15  "Temperature  and  Emissivity  Results  From  a  SSM/T  Retrieval 
Algorithm" 

E.  Kratz  (S  M  Systems  and  Research  Corporation)  and 
H.  Fleming  (NOAA/NESDIS) 

Wednesday,  15  April 

Morning: 

Session  3  -  New  Retrieval  Approaches  and  Products;  Relevance  to  AMSU 
(continued) 

8:30-8:50  "Painless  Extraction:  Subtractive  Temperature  Sensing  from 
Satellite  Radiances  Using  the  Zeta  Transform" 

J.  King  (Phillips  Laboratory,  Geophysics  Directorate) 

8:50-9:10  "SSM/I  Measurements  as  Predictors  of  the  Response  of 
AMSU-A  to  Surface  and  Atmospheric  Phenomena" 

P.  Rosenkranz  (Massachusetts  Institute  of  Technology) 

9:10  Session  4  -  Applications  to  Numerical  Weather  Prediction 

Session  Chairman:  E.  Kalnay,  NOAAINMC 

9:10-9:30  "The  Assimilation  of  DMSP  Retrieval  Products  Into  the  Navy's 
Atmospheric  Prediction  Systems" 

J.  Goerss  and  P.  Pheobus  (NOARL) 

9:30-9:50  "Impact  of  SSM/I-based  Snow/ice  Analyses  in  NMC's  Eta  Model" 

K.  Mitchell  andT.  Black  (NMC);  F.  Messinger  (NOAA/NWS);  N.  Grody 
(NOAA/NESDIS) 

9:50-10:10  "Use  of  SSM/I  Wind  Speed  Data  in  the  Operational  Numerical 
Weather  Prediction  System  at  NMC" 

T.  W.  Yu,  W.  H.  Gemmill  and  J.  Woollen  (NMC) 

10:10  Break 
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Wednesday,  IS  April 
Morning: 

Session  4  -  Applications  to  Numerical  Weather  Prediction  (continued) 

10:25-10:45  "Comparison  Between  SSM/I  and  ECMWF  Total  Precipitable 
Water" 

L.  Phalippou  (ECMWF) 

10:45-11:05  "Meteorological  Impact  of  Surface  Wind  Directions  Measured 
by  Over-the-Horizon  Radar  (OTHR)  and  Wind  Speeds 
Measured  by  the  SSM/I  Microwave  Radiometer" 

R.  Atlas  (NASA  Goddard  Space  Flight  Center),  T.M.  Georges 
(NOAA/Environmental  Research  Laboratory),  and  J.A.  Weinman  (NASA 
Goddard  Space  Flight  Center) 

11:05-11:25  "Observations  of  Flooded  Ice  in  Arctic  Regions" 

A.  Goroch  and  R.  Fett  (NOARL) 

11:25-11:45  "Classification-Based  Rainfall  Estimates  Using  Satellite  Data 
and  Numerical  Forecast  Model  Output" 

C.  Grassotti  and  L.  Garand  (Atmospheric  Environmental  Service) 

11:45  12:05  "The  Influence  of  Cloud  Microphysical  Structure  on 
Spaceborne  Multispectral  Microwave  Observations" 

I.  Rubenstein  (Institute  for  Space  and  Terrestrial  Science) 

Afternoon: 

1:15  Session  5  -  Open  Discussion  (issues,  concerns,  etc.) 


2:00  Session  6  -  Workshop  Issues  ( recommendations ,  workshop  report) 
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DMSP  Conference  Attendee  Address  List 


Rao  Achutuni 

NOAA/NESDIS/ORH/PS  B 
NOAA  Science  Center 
Washington,  DC  20233 
(301)763-8251 

Jill  Aldridge 
NOAA/NESDIS 
Room  3045/FB#4 
Washington,  DC  20233 
(301)763-4310 

Dr.  John  Alishouse 
NOAA/NESDIS 

Room  810,  NOAA  Science  Center 
Washington,  DC  20233 
(301)763-8136 

Dr.  Chuck  Arnold 
Int’l  GEWEX  Project  Office 
409  3rd  Street  S.W.,  Suite  203 
Washington,  DC  20024 
(202)  863-0012 

LtCol.  Edwin  Arrance 

Defense  Meteorological  Satellite  Program 

Space  Systems  Division 

SSD/MIA 

PO  Box  92960 

Los  Angeles  AFB,  CA  90009-2960 
(310)  336-4090 

Dr.  Robert  Atlas 

NASA  Goddard  Space  Flight  Center 
Greenbelt,  MD  20771 
(301)774-6042 

Genia  Bain 

Goddard  Space  Flight  Center 
Greenbelt,  MD  20771 
(301)285-8787 

Dr.  Wayman  Laker 
NOAA 

National  Meteorological  Center 
W/NMC2.  WWB,  Room  204 
Washington,  DC  20233 
(301)763-8005 

Capt.  Leslie  Belsma 

AFSD/AWS-OL-F 

Los  Angeles  AFB,  CA  90009 

(310)416-7720 


Hal  Bloom 

Hughes  STX  Corporation 
4400  Forbes  Road 
Lanham,  MD  20706 
(301)763-4482 

Corey  M.  Boettcher 

S  M  Systems  and  Research  Corporation 
8401  Corporate  Drive,  Suite  510 
Landover,  MD  20785 
(301)763-4384 

Mark  W.  Borden 

HQ/Air  Force  Space  Systems  Division 
(SSD/MWB) 

PO  Box  92960 

Los  Angeles  AFB,  CA  90009-92960 
(310)  336-4494 

Dr.  Donald  J.  Boucher 
The  Aerospace  Corporation 
PO  Box  92957 
Los  Angeles,  CA  90009 
(310)  336-8837 

Dennis  Bungato 
Goddard  Space  Flight  Center 
Greenbelt,  MD  20771 
(301)286-2705 

Ellen  L.  Burdsall 
NOAA/NESDIS 
Room  3045/FB#4 
Washington,  DC  20233 
(301)763-4310 

Chris  M.  Caruso 

NOAA  Ocean  Products  Center 

5200  Auth  Road 

Camp  Springs,  MD  20746 

(301)763-8030 

David  T.  Chang 
AER,  Inc. 

840  Memorial  Drive 
Cambridge,  MA  02139 
(617)547-6207 

Marie  Colton 

Commanding  Officer  Code  70 
FNOC 
Airport  Road 
Monterey,  CA  93943 
(408)  647-4623 


260 


DMSP  Conference  Attendee  Address  List 


Dawn  Conway 

BFEC/NavaJ  Research  Laboratory 

(301)767-2663 

(call  and  check  address) 

Jim  Cornelius 

Commanding  Officer  Code  70 
FNOC 
Airport  Road 
Monterey,  CA  93943 
(408)  647-4398 

Dr.  James  Dodge 
NASA  Headquarters 
600  Independence  Avenue,  S.W. 
Washington,  DC  20546 
(202)  453-1680 

David  R.  Donahue 
SMSRC/DMSD 
Room  3204,  FB#4 
Washington,  DC  20233 
(301)763-4383 

Paul  DouPont 

Naval  Research  Laboratory  (NRL) 
Code  8311 

Washington,  DC  20375-5000 
(202)  767-5417 

Henry  F.  Drahos,  Jr. 

NOAA/NESDIS 

Room  3212/Fed  Building  #4 

Washington,  DC  20233 

(301)763-4380 

Lt.  Comm.  Chris  Durbin 
British  Navy 
NOM  (OP) 

Lacon  House 
Theobalds  Road 
London  WC1  X8RY 
ENGLAND,  U  K. 

071-430-6739 

Vincent  Falcone 

Phillip  Laboratory  Geophysics 

Directorate  (PL/GP/LY) 

Hanscom  AFB,  MA  01731 
(617)  377-4029 


Ralph  R.  Ferraro 

NOAA/NESDIS 

Satellite  Research  Laboratory 

Washington,  DC  20233 

(301)763-8042 

Joseph  V.  Fiore 
SMSRC/DMSD 
Room  3204,  FB#4 
Washington,  DC  20233 
(301)763-8211 

Dr.  Henry  E.  Fleming 
NOAA/NESDIS 

Room  810,  NOAA  Science  Center 
Washington,  DC  20233 
(301)763-8136 

William  H.  Gemmill 

NOAA/NMC 

5200  Auth  Road 

Camp  Springs,  MD  20746 

(301)763-8133 

Dr.  Morton  Glass 
PL/GPAS 

Hanscom  AFB,  MA  01731 
(617)  277-2946 

Mitchell  Goldberg 
NOAA/NESDIS/SRL 
NOAA  Science  Center 
Washington,  DC  20233 
(301)763-8136 

Dr.  Andreas  Goroch 
Naval  Research  Laboratory 
Monterey,  CA  93940 
(408)  647-4725 

Christopher  Grassotti 
Atmospheric  Environment  Service 
2121  Trans-Canada  Highway 
Dorval,  H9P  1J3 
Canada 

(514)  421-4749 

Shelton  Green 
Code  SSD/MIN 
Headquarters  Space  Division 
LA  AFB,  CA  90009-2960 
(310)336-4408 


261 


DMSP  Conference  Attendee  Address  List 


Dr.  Michael  Griffin 
PL/GPAS 

Hanscom  AFB,  MA  01731 
(617)  547-6207 

Dr.  Norman  C.  Grody 
NOAA/NESDIS 

Room  712/NOAA  Science  Center 
Washington,  DC  20233 
(301)763-8136 

Arnold  Gruher 
NOAA/NESDIS 
NOAA  Science  Center 
Washington,  DC  20233 
(301)  763-8053 

Dr.  Kenneth  R.  Hardy 

Lockheed  Missile  &  Space  Co.,  Inc. 

0/91-01  B/255 

325 1  Hanover  Street 

Palo  Alto,  CA  94304-1191 

(415)  424-2825 

David  Hogan 
GE  Astrospace  Division 
PO  Box  800 

Princeton,  NJ  08543-0800 
(609)  490-3355 

Dr.  James  P.  Hollinger 

Naval  Research  Laboratory  (NRL) 

Code  8311 

Washington,  DC  20375-5000 
(202)  767-3398 

Nicholas  Horvath 
General  Electric 
PO  Box  8048 

Attn:  Bldg-A,  Room  33A30 
Philadelphia,  PA  19101 
(215)531-1385 

Dr.  Keith  0  Hutchison 
Ixickheed  Missiles  &  Space  Co. 
6800  Bourleson  Road 
T2-30,  Bldg.  310 
Austin,  TX  78744 
(512)  386-1678 


Ronald  G.  Isaacs 
Atmospheric  and  Environmental 
Research,  Inc. 

840  Memorial  Drive 
Cambridge,  MA  02139 
(617)  547-6207 

John  Jandik 
ITT-ACD  -  Code  646 
1919  West  Cook  Road 
PO  Box  3700 
Ft.  Wayne,  IN  46801 
(219)  487-6357 

Dr.  Eugenia  Kalnay 
NOAA/NMC 

National  Meteorological  Center 
Washington,  DC  20233 
(301)763-8005 

Vaishali  Kapoor 
NOAA/NESDIS 
FB#4,  Room  3212 
Landover,  MD 
(301)763-4384 

Dr.  John  P.  Kerekes 
MIT  Lincoln  Laboratory 
L-143 

244  Wood  Street 
Lexington,  M A  02173-9108 
(617)  981-0805 

Dr.  Jean  I.F.  King 
Geophysics  Directorate 
Phillips  Laboratory 
Hanscom  AFB,  MA  01731-5000 
(617)  377-2976 

Eugene  J.  Kratz 

SM  Systems  &  Research  Corporation 
8401  Corporate  Drive,  Suite  510 
Landover,  MD  20785 
(301)763-4383 

Dr.  Lowell  Krawitz 
GE  Astrospace  Division 
PO  Box  800 

Princeton,  NJ  08543-0800 
(609)  490-3386 


262 


DMSP  Conference  Attendee  Address  List 


Herb  Kroehl 

NOAA/NGDC 

(303)497-6121 

Frank  LaFontaine 
USRA 

ES43  Bldg.  4481 
Huntsville,  AL  35812 
(205)  544-4348 

Dr.  John  A.  Leese 

International  GEWEX  Project  Office 

Suite  203 

409  Third  Street,  S.W. 

Washington,  DC  20024 
(301)405-4851 

Dr.  Robert  McClatchey 
Phillip  Laboratory  Geophysics 
Directorate  (PL/GPA) 

Hanscom  AFB,  MA  01731 
(617)  377-2975 

Dr.  Larry  McMillin 
NOAA/NESDIS 
Satellite  Research  Laboratory 
Washington,  DC  20233 
(301)763-8136 

Dr.  Kenneth  E.  Mitchell 
NOAA/NMC 

Room  204,  NOAA  Science  Center 
Washington,  DC  20233 
(301)763-8161 

Douglas  Moore 
SMSRC/DMSD 
Room  3204,  FB#4 
Washington,  DC  20233 
(301)763-8211 

James  Morrissey 
PL/GPAS 

Hanscom  AFB,  MA  01731 
(617)  377-4405 

Lt.  Cdr.  Cynthia  Nelson 
Office  of  the  Federal  Coordinator 
for  Meteorologiccl  Science  and 
Supporting  Service 
Suite  900 

6010  Executive  Bldg. 

Rockville,  MD  20852 
(301)443-8704 


Tom  J.  Neu,  Capt.,  USAF 
HQ  AFGWC,  SYS/Satellite 
106  Peacekeeper  Drive,  STE  2N3 
Offutt  AFB,  NE  68113-4039 
(402)  294-7444 

Donald  Norquist 
Phillips  Laboratory  (GPAP) 
Hanscom  AFB,  MA  01731-5000 
(617) 377-2962 

Charles  Novak 
NOAA/NESDIS 
Rm  3212/FB4 
Washington,  DC  20233 
(301)763-4380 

Dr.  George  Ohring 
AFGWC/NOAA/NESDIS 
Room  712/NOAA  Science  Center 
Washington,  DC  20233 
(301)763-8078 

Capt.  John  Pereira 
AFGWC/NESDIS  Liaison 
Room  0308,  FB#4 
Washington,  DC  20233 
(301)763-4384 

Patricia  A.  Phoebus 
Naval  Research  Laboratory 
Atmospheric  Directorate 
Monterey,  CA  93943 
(408)  647-4700 

Laurent  Phalippou 
ECMWF 
Shinfield  Park 
Reading  RG2  9AX 
United  Kingdom 
+44  734  496  656 

William  Pichel 
NOAA/NESDIS 

Room  102/NOAA  Science  Center 
Washington,  DC  20233 
(301)763-8231 


263 


DMSP  Conference  Attendee  Address  List 


John  Pickle 

Irene  Rubenstein 

Atmospheric  and  Environmental 

Earth  Observation  Laboratory/ 

Research,  Inc. 

Institute  for  Space  and  Terrestrial 

840  Memorial  Drive 

Science 

Cambridge,  MA  02139 

4850  Keele  Street 

(617)  547-6207 

North  York,  Ontario 

CANADA  M3J  3K1 

Dr.  Gene  Poe 

Aerojet  Electrosystems 

(416)  665-5410 

1 100  W.  Hollyville 

John  Ruggiero 

PO  Box  296 

Code  SSD/MIN 

Azuza,  CA  91702 
(818)812-1154 

Headquarters  Space  Division 

LA  AFB,  CA  90009-2960 

(310)  336-4410 

John  Pritchard 

NOSA/NESDIS 

Glenn  Sandlin 

Room  102/NOAA  Science  Center 

Naval  Research  Laboratory  (NRL) 

Washington,  DC  20233 

Code  8311 

(301)763-4310 

Washington,  DC  20375-5000 

(check  address) 

(202)  767-3443 

Dr.  P.  Krishna  Rao 

Noelle  A.  Scott 

NOAA/NESDIS 

ARA/LMB 

Room  701/NOAA  Science  Center 

Centre  National  de  le 

Washington,  DC  20233 

Recherche  Scientifique 

(301)763-8127 

Ecole  Polytechnique 

Palaiseau  Cedex  F-91128 

Dr.  C.R.Nagaraja  Rao 

FRANCE 

NOAA/NESDIS/SRL 

NOAA  Science  Center 

(33)  1  69334532 

Washington,  DC  20233 

Greg  Sharfen 

(301)763-8136 

National  Snow  and  Ice  Data  Center 
University  of  Colorado 

Anthony  Reale 

CIRES 

NOAA/NESDIS 

Box  449 

Room  3212/Fed  Building  #4 

Boulder,  CO  80309 

Washington,  DC  20233 

(303)492-6197 

(301)763-4380 

Mike  Shields 

W.S.  Richardson 

MIT/Lincoln  Labs 

NOAA  Ocean  Products  Center 

244  Wood  Street 

5200  Auth  Road 

Lexington,  MA  02173-9108 

Camp  Springs.  MD  20746 

(617)  981-4919 

(301)763  8030 

David  Spangler 

Dr.  Phillip  Rosenknmz 

Code  4211 

MIT 

Naval  Research  Laboratory 

77  Massachusetts  Avenue 

4555  Overlook  Ave. 

Cambridge,  MA  02 1 39 

Washington,  DC  20375-5 000 

(617)  253-3073 

(202)  767-5417 

264 


DMSP  Conference  Attendee  Address  List 


Steve  Swadley 

Bruce  Whistler 

Aerojet  Electrosystems 

SMSRC  (DMSP) 

1 100  W.  Hollyville 

NOAA/NESDIS 

PO  Box  296 

FB  #4,  Room  3212 

Azuza,  CA  91702 

Landover,  MD 

(818)812-1426 

(301)763-3875 

Pamela  M.  Taylor 

Harold  Woolf 

NOAA/NESDIS 

U.  Wisconsin 

Room  3212/Fed  Building  #4 

NOAA/NESDIS 

Washington,  DC  20233 

Room  203 

(301)763-4380 

1225  W.  Dayton  Street 

Madison,  WI  53706 

Bruce  Thomas 

The  Aerospace  Corporation 

(608)  264-5325 

PO  Box  819 

Man  Li  Wu 

Bellevue,  NE  68005 

Goddard  Space  Flight  Center 

(402)  291-3140 

Sher  Wagoner 

Greenbelt,  MD  20771 
(301)286-4087 

NOAA/ERL/FSL 

Tsann  W.  Yu 

R/E  S4  3100  Marine 

Development  Division 

Boulder,  CO  80303 

NOAA/National  Meteorological  Center 

(303)  938-2085 

5200  Auth  Road 

Camp  Springs,  MD  20746 

Bob  Wagoner 

NOAA/NESDIS 

Room  3045/FB#4 

Washington,  DC  20233 
(301)763-4310 

Marshall  Waters 

NOAA  Ocean  Products  Center 

5200  Auth  Road 

Camp  Springs,  MD  20746 
(301)763-8030 

Captain  Steve  Way 

SSD/MIB 

PO  Box  92960 

Los  Angeles,  CA  90009-2960 
(316)  336-4071 

Allan  Weiner 

SAF/SS 

Pentagon  -  Room  4C-1052 

Washington,  DC  20330 
(703)614-1772 

(301)763-8133 

265 


AUTHOR  INDEX 


Alishouse,  John  C . 77 

Atlas,  Robert . 204 

Baker,  Wayman  . 123 

Black,  Thomas  L . 169 

Boettcher,  Corey  M . 72 

Bohlson,  J . 25 

Borden,  Mark  W . 25 

Boucher,  Donald  J . 11,  17 

Burdsall,  Ellen  L . 33 

Caruso,  Chris  M . 96 

Colton,  Marie  C . 58 

Cornelius,  C.  James  . . 58 

Donahue,  David  . 38,  72 

Falcone,  Vincent  . 11 

Ferraro,  Raloh  R . 77,  82 

Fett,  Robert  W  . 227 

Fiore,  Joseph  V.,  Jr . 66 

Fleming,  Henry  E . 137,  141 

Garand,  Louis  . 232 

Gemmill,  William  H . 96,  183 

Georges,  T.M . 204 

Glackin,  D.L . 25 

Goerss,  James  S . 161 

Goodman,  H.  Michael . 108 

Goroch,  Andreas  K . 227 

Grassotti,  Christopher . 232 

Griffin,  Michael  . 11 


Grody,  Norman  C . 82,  169 

Hollinger,  James  P . 2 

Isaacs,  Ronald  G . 11 

Kelly,  Frank  . 25 

King,  Jean  I.F . 145 

Kratz,  Eugene  J . 137,  141 

Mandt,  G . 25 

McMillin,  Larry  . 132 

Messinger,  Fedor  . 169 

Mitchell,  Kenneth  E . 169 

Morrissey,  James  . 11 

Neu,  Thomas  J . 45 

Phalippou,  Laurent  . 195 

Phoebus,  Patricia  A . 161 

Pichel,  William  . 89 

Pickle,  John  . 11 

Reale,  Tony  L . 38 

Richardson,  W.S . 89 

Rosenkranz,  Phillip  W . 152 

Rubenstein,  Irene  G . 245 

Taylor,  Pamela  M . 106 

Thomas,  Bruce  H . 1 1,17 

Waters,  Marshall  . 96 

Weinman,  James  . 204 

Woollen,  J . 183 

Yu,  Tsann  W . 183 


V 


266 


